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NOTICE TO AUTHORS 


1. Communications.—Papers must be communicated to the Society bya Fellow. They 
should be accompanied by a summary at the beginning of the paper conveying briefly the 
content of the paper, and drawing attention to important new information and to the main 
conclusions. ‘The summary should be intelligible in itself, without reference to the paper, 
te a reader with some knowledge of the subject; it should not normally exceed 200 words 
in length. Authors are requested to submit MSS. in duplicate. These should be 
typed using double spacing and leaving a margin of not less than one inch on the 
left-hand side. Corrections to the MSS. should be made in the text and not in the 
margin. Unless a paper reaches the Secretaries more than seven days before a Council 
meeting it will not normally be considered at that meeting. By Council decision, MSS. of 
accepted papers are retained by the Society for one year after publication; unless their 
return is then requested by the author, they are destroyed. 


2. Presentation.—Authors are allowed considerable latitude, but they are requested to 
follow the general style and arrangement of Monthly Notices. References to literature 
should be given in the standard form, including a date, for printing either as footnotes or in 
a numbered list at the end of the paper. Each reference should give the name and 
initials of the author cited, irrespective of the occurrence of the name in the text (some 
latitude being permissible, however, in the case of an author referring to his own work). 
The following examples indicate the style of reference appropriate for a paper and a book, 
respectively :— 

A. Corlin, Zs. f. Astrophys., 15, 239, 1938. 

H. Jeffreys, Theory of Probability, 2nd edn., section 5.45, p. 258, Oxford, 1948. 

Alternatively, the Harvard reference system may be used. 


3. Notation.—For technical astronomical terms, authors should conform closely to the 
recommendations of Commission 3 of the International Astronomical Union (Trans. 
1.A.U., Vol. V1, p. 345, 1938). Council has decided to adopt the 1.A.U. 3-letter abbrevi- 


ations for constellations where contraction is desirable (Vol. IV, p. 221, 1932). In general 
matters, authors should follow the recommendations in Symbols, Signs and Abbreviations 
(London: Royal Society, 1951) except where these conflict with I.A.U. practice. 


4. Diagrams,—These should be designed to appear upright on the page, drawn 
about twice the size required in print and prepared for direct photographic 
reproduction except for the lettering, which should be inserted in pencil. 
Legends should be given in the manuscript indicating where in the text the 
figure should appear. Blocks are retained by the Society for 10 years; unless the author 
requires them before the end of this period they are then destroyed. 


5. J ables.—These should be arranged so that they can be printed upright on 
the page. 


6. Proofs.—Costs of alteration exceeding 5 per cent of composition must be borne by 
the author. Fellows are warned that such costs have risen sharply in recent years, and it 
is in their own and the Society’s interests to seek the maximum conciseness and simplifi- 
cation of symbols and equations consistent with clarity. 


7. Revised Manuscripts.-When papers are submitted in revised form it is especially 
requested that they be accompanied by the original MSS. 


Reading of Papers at Meetings 
8. When submitting oes authors are requested to indicate whether they will be 


willing and able to read the paper at the next or some subsequent meeting, and approxi- 
mately how long they would like to be allotted for speaking. 


9. Postcards giving the programme of each meeting are issued some days before the 
meeting concerned. Fellows wishing to receive such cards whether for Ordinary 
Meetings or for the Geophysical Discussions or both should notify the Assistant Secretary. 
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ROYAL ASTRONOMICAL SOCIETY 


MEEK TING OF 1956 MAY 
Professor Sir Harold Jettreys, President, in the Char 
Che election by the Council of the following Fellows was duly confirmed : 
(seorge Allison Kiby, M.S« 2 Shortland Street, Khandallah, Wellington, 
New Zealand (proposed by Kk. D. Adam 
*Michael Godfrev Fkishet B.Sc.. West Wood, Carleton, Pontefract. Yor! 
shire (proposed by ¢ W Allen) 
Walter J. Miller, S. J., Fordham University Astronomical | 
Ne \ York | >) \ (proposed by R \ d kk Woolk Vv) 
Walter Orr Robert Director of the High Altitude Observator Boulder, 


i Ty 


iboratory 


Colorado, U.S.A (proposed DV ». 4 hapman) ind 
Ilenry Peter krancis Swinnerton-Dyer, M.A lrinity ollege, Cambridge 
(propo ed Dy B | | Pavel) 


lhe election by the Coune 


onfirmed 


John Keith Welwyn Davi 
(propo ed by Z Kopal) 


seventy-foul 


last 1 ecting imclue 


(y rar \ micourlcul Oo} the Planet Wa 





THRE DAMPING OF THE VARIATION OF LATITUDI 


Harold Jeffreys 


Summary 


pointed out a mistake in an estimate of the viscosity 
damping of the 4-! thly variation of 
riticism 1s diminished when account 1 
and correction of this leads to the con- 
just be due to dissipation he core, but 
innot be due to viscosit there Ihe l 
ysity in the shell at depths greater th: 


would be about 10*° ¢.g.s. unit or to elastic after- 


In the latter case a wide range of parameters in the law is permis- 


In the third edition of The Earth 1 used three data relevant to the admissible 
ity of the karth’s core. 
(a) ‘The transmission of sharp longitudinal pulses through the core forbids 
a kinematic iscosity greater than 10° cm*/sec.* 
(4b) An extrapolation based on Bridgman’s work sugyeste d that 1 cm?/sec is 
probably of the right order of magnitude. 
(c) On the assumption that damping of the 14-monthly variation of latitude 


is due to viscosity in the core, I found a viscosity of about 10° cm?/sec. | 
(d) | found also that within a wide range of possibilities elastic afterworking 

in the shell could account for the damping 
Bondi and Gold in an interesting paper§ have argued that no type of 
dissipation in the core can account for the damping of the variation of latitude 
hey say, however, in the first paragraph on p. 42 that I considered it 
probable, or at any rate possible”’ that viscosity in 1@ core accounts for 


in view of the above comparison thi tatement ems exaggerated In 


particulas it we accept (>), and | see no reason to doubt it, their result © tal 


it concerns Vi itv Is an Immediate con equen 

\s they point out there is a serious mistake 1n (¢). ‘hey are correct in aying 
that my discussion confuses two meanings of (/, m), which at first denote direction 
cosines of the instantaneous axis with respect to the geographic axis and late: 
denote those of the geograph axis with re pect to the mean frame of referenc 
llowever, the argument contains a second mistake, which they mention but do 
not explain in detail his is worth attention, particularly ill be 
relevant to further discussions. In the simple case of a rigi hell and a 
homogeneot core | took a general rotation of the hell, referred 
rotating axe expressed by direction cosines (/,m), and put l 


additional motion of the core 1s re pres nted by a compl juantity 


*H ettre' / rf i P| 244 t/ 
MN., 838, 41-46 H. Jeftre 
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not give a inple rotation, being atlected by the ellipticity but the resulting 
displacements differ trom those given by a rotation within a tactor of order e, 
where e€ is the ellipticity of the core boundary Lhe ratio « (1 —€) + Of) 
he slip depends only on (hus with an oblate boundary ellipticity decreases 
the slip, but prolatenes would increase it In any case the ellipticity Is quit 
unimportant (i his is not true tos precession and nutation there the ellipticity 
has a dominating etlect on account of the proximity of two tree periods.) 

For an clastic shell the elastic displacements of the boundary produce similat 
normal displacements to those produced by a rotation of the shell, and it was 
convenient to introduce a fictitious rotation that would produce the same 
normal displace ments as the actual rotation and the elastic displacements together 
Unfortunately | supposed in The Earth that it would also produce the same 
tangential displacement , and this ts false lhe motion of the core is the 
resultant of those produc d by and ¢, with now about (1 €)¢ In the 
case | solved tor an elasti hell was about o 1 4¢ With a factor e the slip 1S 


equal to that produced by a rotation ¢ (i Ct) and hence to the same accuracy 


to that produced by ¢ For an elastic shell as for a rigid one, therefore, the slip 


is practically the same as if the core did not move. Consequently the factor o-14 
should be suppressed in the ratio of the slip to ¢ Incidentally Bondi and Gold 
(second paragraph on p. 45) are wrong in supposing that the rotation of the shell 
varies trom C’ to with radius; the mean rotation of a level surface 1 
throughout 

In equation (13) of The Earth, p. 246, | suppress the factor 0-02 (= 0-14") 
and apply the factor 1/420 to /, m [his divides the expression, for given v, by 
about 3600. ‘hus to fit my estimate of the damping would make the thicknes 
of the boundary layer about 1400 km and the kinematic viscosity about 10!” cm?/se 
‘his viscosity is inconsistent not only with estimate (+), mentioned above, but 
with the purely geophysical estimate (a). Further, at such thickness, ordinary 
boundary layer theory would need a considerable correction, and the consideration 
of Bondi and Gold becomes relevant, that too high a cosity would drag the 
core around with the she I] and the dissipation would tall off avain. It does howe vel 
seem, apart from other considerations, that suitable viscosity could give damping 
of the right order of magnitude 

| think that some points of Bond: and Gold discussion need further 
clarificatior particularly im \ f remark made vhen the papel! Wa read,* 


thei n plified model g rr tl naximum damp vy a time tactor 


rotation peed } the dynamic il ellipticity of the shell 
Ince elliptre ity of the core is neglected) and & is the 
of inertia of the core and shell since the free period 
1, this amounts in a period to about exp ( r), and 
the amplitude to e! of its original value is about 
maximum damping at first sight, would be consistent 
vith Walker and Young estimate 


are whether the uppo ed maximum 
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mum, and if so, how far the actual value 1 
¢ theorems of minimum rates of dissipation for 
a COSITY 


likely to fall 


steady motion 
and electric currents, and since the period in th 
are presumably applicable in 


ir maximum 1s an 
| motion would be 


cad 


is Ca 
this case I therefore se 


over-estimat« As Bondi and Gold say 
substantially different from theirs and the 
to be nearly the maximum consistent with the typ 
Whether thi instance, electro 
iy be left to personal judgment. Without 
it seen to be po ible to 


no reason 


Is an argument against, for 


I 


more detailed 
say that the damping cannot 
in the core but might just possibly be due to electromagnet 


been 


uggested that the upper bound given by (a) might be irrelevant 
iriation of latitude, for 


two reasons. 
on. the period. 


First, the 


viscosity might depend 


Hlowever, for normal liquids the 


letermined either from steady motion in <¢ 


plate ind the two methods agree 
cocthcient of viscosity 


cosity can bi 
1 tube or from torsional oscillations of 
Secondly, the discovery of the 

for liquids introduces a doubt as to whether an estimate 


second 
of viscosity from longitudinal waves would be applicable to a mainly distortional 


movement. ‘The second coethcient, however, expresses a resistance to expansion 
ind will give extra damping of longitudinal waves without 
di tortion thu . es 


longitudinal waves 


a corresponding 
iddition for important, the estimate of the 
I 


viscosity given by will be too high 


ordinary 
trenythened 


and the argument 


| did not give in The Earth an explicit estimate of the viscosity in tl 


ceded to account for the damping 


hell 
lor elastico-vi ity 


247 


ot re | ixation 


‘his is easily found from (1g) on p 


is infinite and the equation leads to a time 
of rigidity 


$0 
2600 days, 
etlective coefficient of viscosity would be of the 
ot out of the question I he 
ol 


order of 10 


YIn/Cm Sc 
gravity anomalies would be consistent 
trength at de pths over 700 km at the cost of a str neth of about 
cm* above that level. An infinite set of 


solutions attributing the 
ifterworking would also be possible 





OBSERVATIONS FROM AN AIRCRAFT OF THI 
ZODIACAL LIGHT AT SMALL ELONGATIONS 


D. E. Blackwell 
(Received 1956 April 30) 


Summary 


1 for polarization data for the zodiacal light between « 
liscussed It is shown that sucl 
high altitude, and « pecially from 
urements of brightne in absolut 
clongation a and 31 alony the 
vht of gooo ft over the South Pa 
aX! was then vertical 
of tl ixis Of symmetry on Ut 
interpreted in terms of electron den 
id the effect on the computed electron densiti 
erpretation are discussed Che angle between the 
clectron component and the ¢ liptre for thi ranyve of 
le than 1 


/ 


roduction During the last few years there has been a renewed interest 
mtiacal light phenomena, perhaps partly stimulated by a paper of Whippk 
h | discusses for the first time the po sibilitv that the phenomena 


ntirely to the scattering of solar radiation by free electron Since then 


nerally assumed that at least part of the zodiacal light is due t 


the remainder being due to scattering bi lhe reason 

that the observed polarization for all elongation matles 
too great to hye ittributable to dust scattering alone It i 
inite decision about this, however, because our knowledge 
ist component is still very scanty, and will not be increased until 1 


are available Briectly we may reason j ollow Lhe 
iwht differs slightly from that of the crag jar d 
of the light 1s due to electron scatte1 
itferent from that of the Sun, Rayl 
proc ind the effective particl 
il that the high polarization a 
Van de Hulst (4) ha 
ittering by larg partick 
30. a Value of per c 
is not known, but 
dust component is not gre 
other hand, the observed pola 
f 9] lee ‘ 


ne ord ol 1) 
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Lhe important problem has therefore emerged of deducing from zodiacal 


light data the electron density as a function of distance from the Sun. ‘The 
general principles involved are well known Che contributions from dust 
and electrons have first to be separated Ihe simplest kind of observation that 
can be used tor this purpose is that of polarization. Denoting the brightnes 


of the dust and electron components as F and K respectively, and thei polariza 


tions a p’ ind p’ we have for the total observed polarization 


hk 
Po= Fy’? 


h 


Ra 5 (1) 
If p” is known, or can be assumed zero this equation can be solved by successive 
ipproximation and both A and F calculated lor the solar corona, within the 
range of elongation o S p’ can without much error be assumed zero (7). Sut 
reliable observations of the polarization of the zodiacal light have hitherto existed 
only tor elongation yreatel than a6 and tor these clongation ve have no 
nowledge whatsoever of the value of p’ It is true that plausible prediction 
in be made using idealised models, but for a completely reliable prediction it 
ould be necessary to know accurately the distribution of particles in size, thei 
distribution in space, their shapes and orientations, their surface properties and 
composition Ihe upposition of Whipple (6) that th partick ire highly 
porous and of low density complicates the problem even further here ar 
reasons, even, for supposing that the polarization of the dust component is negative 
(5) kaced with thi problem there ar only two course open to us, if the 
polarization data are to be used to derive electron densiti We must either 
measure the polarization of the dust component, or confine work to those elonga 
tion vhere it can be ately neglected In thi papel ce consider the econd 


COLT S¢ 


Need for observations at small elongations In Fig. 1 howna curve which 
vive the maximum likely polarization of the dust component It has bee 
constructed by taking the kind of value suyy ted by Whipp (1) and van de Hulst 
(4) for clongation 6« ind the value expected for Rayleigh seattering for small 


longation nd making a plau ible inte rpolation lhe reliability of an electror 
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/ 


density calculation which uses me 


asured values of the polarization at a particular 
elongation, depends upon the ratio R p* E af p* now being the likely value of 


the dust polarization [he quantity R—1 has been plotted in Fig. 2, using the 
results for p' 


deduced from observations described later in this paper, and 
those of Behr and Siedentopf (2), together with the plausible values for p* just 
suggested. No reliance can be placed on electron density deductions at 
elongations for which R—1< 0. ‘This figure demonstrates that, for elongations 
in the range 36 >« 22 , this ratio is never less than 2-7, which means that the 
error of calculation of electron de nsity, if p’ is assumed zero, 1s not likely to be 
more than 37 percent. In particular, observations in the region «<7 and around 
€=30 can be interpreted with a maximum likely error of 30 per cent, while beyond 
elongation 60 interpretation in terms of electron densities becomes practically 
worthless (unless there is prior knowledge of p*) fhe minimum clongation 
for which reliable data for the polarization of the zodiacal light have hitherto 
been obtained ts « 1. (he corresponding maximum elongation in the sola 
corona IS 5 ‘hese clongations correspond to minimum distances trom the Sun 
of 20R., and 122R and between these distances there are no electron density 
values available. It is evidently very important that measures of the zodiacal 
light should be brought inwards at least to an elongation of about 20 ( ‘ 


this will lead to reliable electron density values being obtained within the distances 


range 123k... to74Rk~. In this paper we shall show that reliable measures can 
74 pat 





Aind of ol lions needed I he prime purposs ot the c inne 


to obtain reasonably accurat 


| 
biti 
ent ot both total brightn $ and polar Zation are 
measurements of the hbrightne of the zodiacal 


e of Roach etal (3) how that ther isa relat on oltthetorn 


and elongation which repre ov 
between elongation 


vhich link the solar corona 
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how that this relation reproduces the observations with an accuracy of the order 
of § per cent at least over the range of elongation 6 —70 \s the error in the 
interpretation of the observations may be as high as 25 per cent there 1s therefore 
no need to measure the brightness of the zodiacal light in the very difficult region 


20 ‘ 20 , but it is sufficient to measure it over the range 30 e<so and 


extrapolate to smaller elongation It will become apparent later that this simpli- 


hes the problem considerably he polarization, however, presents a completely 
ditterent problem \s the values vhich are plotted in hig. Q show, it 1s not 
possible to make a reliable interpolation between the solar corona and the zodiacal 
light in the region 20 ~—e- 235 , and it is therefore essential to obtain accurate 
measurements. 

4. Previous observations.—Vhere have been many investigations of zodiacal 
light brightne 3 (mong the earlier work the most extensive and reliable 1s 
that don photoc lectrically by Elvey and Roach (8) between elongations $O and 
yo. In work of this kind the great difficulty is to decide upon the brightness 
of the background, and in a later paper by Roach et al. (6) this problem 1s discussed 
in detail Here, values of zodiacal light brightness between elongations 30 and 
100 are vive! The difthculties of measurement are shown by a comparison with 
the resu Schr and Siedentopt (2): at the elongation of 35 the brightnes 
viveln by Roach 1s twice that piven by he hy and Si de ntopl [his may be a real 
difference in brightne but Regene (g) in some very careful measurement 
does not detect any systematic variation with time over a period of 15 month 

here have, untortunately, been few quantitative investigations of polarization 
lhe first was by Dutay (10) who used a photographic method between elongation 

ind g¢ lhe second was by Behr and St dentopt (2) who useda photoc lectri 
method between elongations 35 and go Ihe two sets of data agree as to order 
of magnitude, but whereas at 40 clongation Dutay obtains the value 12°> per cent 

and Siedentopt derive 23 per cent lake-uchi (11) has made measurement 


‘ 


myation 6« .o observations exist for the Important range of elongation 


(C'‘onditions for observations in this region lo observe the zodiacal light 

mall clongations it is clearly necessary to observe as soon after sunset and a 

to the horizon as possible Suppose that unde exceptional conditions 

is possible to ob ve within 4 of the horizon. Fig. 3 shows the brightnes 

of the zodiacal lig ind of the sky that may be expected at altitude 4 for different 

lepressior the Sun below the horizon he zodiacal light curve show 

the observed brightne uncorrected tor atmo pheri extinction Lhe data tor 

the sky brightne in this figure, which are approximat only, have been obtained 

photograph illy from an aircraft by methods described late Some ground 

observations at smaller de pre ions of the Sun below the horizon have also been 

ised in constructing the diagram, all of which refers to an ettective wave-length 
of 6200 A 

lhe diagran lemon it under these conditior the zodiacal light 


equal in. brightne » th \ ackground at an elongation of 22 If the 


1 


itmospheri conditiotr r atistactory it should thus be possible to make 
observations to at least this clongation Llowever, the sky brightne increase 
very rapidly with decreasing depression of the Sun below the horizon and it ts 
clear, first, that good tin vy of the observations is necessary, and econdly, that 


the mille | clonyation al vhich ob ervation can be mad is harply defined. 
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Che second intersection of the two curves at elongation 20) 


is of no significances 
because there is no need to observe at this clongation so close to the horizon 


wriation of sky brightne and zodiacal 


helo the horizon 1//] measurements + 


It is clearly essential to have exceptionally good observing conditions tot 


quantitative observations of this nature. In particular, it is not casy to make 


uch observations from the ground because the light path passes very close to th 
l-arth’s surface for a great distance, and is therefore much influenced by low lying 
dust and haze. ‘lo avoid this difficulty it is desirable to observe ove 
cither from an isolated mountain peak or from an aircraft Lhe 


f an aircraft is that it enables observation to be made 


a sca horizon 
great advantage 


over open ocean, fat from 


large land masses, where the air is virtually free from dust and smoke Under 


these conditions the atmospheric extinction is small (as is demonstrated later) 


and its variation with altitude regular. From an elevated position usetul 


observations can be made down to an altitude ¢ above the depre 


where ¢ 1s the depression of the horizon 


sed horizon, 
In this proble m there is little to be 
gained from observing at a height* greater than 1ooooft; at this height the 
observer 1s already above most haze and dust, and observations can be made te 


vithin 2. of the horizon before the haze layer is penetrated by the light path 


lo reduce as much as possible the effect of atmo pheru extinction the 
ecliptic should be nearly vertical. 


6. Observations he observations of the zodiacal light de 


1) ide 
¢ lip 


cribed ere VCC 
as part of the programme of an ¢ xpedition to Fiji to observe the total sola 


e of IQ ss June 20 hese observation were planned to be complementary 


to the e lips« obs« rvations, all being made at about the same tims Lhey were 


made 
over sea in an open Sunderland aircraft of the Royal New Zealand Air borce at 
i he ight of gooo ft and in the approximate po tion 


Latitude 16 ‘> 


Longitude 17 \W 


In thas } iper the 


ide reters to ang ot line 
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i 


‘Two flights were made for this purpose, both in excellent weather conditions. 


No photographs were taken near sunset because of the confusion at that season 


between the Milky Way and the zodiacal light. 


Viean depr 
of Sun below 


rizor 


5900 6100 6300 
wavelength 


imilar in design to that used to observe the total solar eclip 


~ 


a description of this camera has already been published (7) 


1954 June 3 
I here is no need in work of thi 
measuring the brightnes of a featureless surtace, and Im any 
the irregular motion of the aircraft during the relatively 


kind to use high quality lenses; we are simply 


Cast vithout 


claborate apparatu } 


long exposure oon vitiates good definition. Three simple len es of white 
gla were used, each of diameter 7-1cm and focal length 12°gcim, vith relative 


aperture f/i-d, giving a pl ite cale of o $4 permm and a field of 60 ach len 


used in conjunction with a Chance OY1 filter, giving the spectral sensitivity 


va 
} for a source at Hooo dev. K lhe effective wave-length*® is clos 


hown in Fig 
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to 6300 A l'wo of the three lenses were also backed by suitably orientated sheets 
of HN22 Polaroid filters [he three lenses formed three separate images on to 
one Ilford HP3 plat Under these conditions an exposure of 5 min gave an 
excellent direct photograph, but an underexposed polarization photograph. 
Photographs | and II were given exposures of 15 min and 5 min respectively, 
and IIl and IV were given 8 min and 4min. Nos. | and II were taken while the 
zodiacal light was still almost uncontaminated by sunlit sky near the horizon, 
[1] and IV were taken later to obtain the polarization of the dawn sky. 

Several hours later on the ground an absolute calibration from a sub-standard 
lamp was put on to another plate of the same batch through a quartz step wedge, 
and the zodiacal light plate and calibration plate developed together Ihe 
substandard lamp was later compared with a second standard which had been 
calibrated by the National Physical Laboratory ‘he brightness of the zodiacal 
light could thus be obtained in absolute units, and by using the value for the 
mean brightness of the Sun given by Minnaert (12) the zodiacal light could bs 
related directly to the average brightness of the solar dis Precautions mentioned 
In a previous paper (7) were taken to ensure photometric accuracy Particular 
attention has b paid to the variation of Camera sensitivity with distance from 
the optical axi 

‘The camera was mounted on a table behind an open hate h in the side of the 
aircraft \ crude atte mipt to guide was made, but much reliance wa pl iced on 
the stability of the aircraft It was found that in reasonably good flying condition 
this results in nearly round star images rather more than 1 in diameter for an 
exposure time of 5 min Such definition is amply good enough for this problem 

Veasurement of atmospheric transmission he atmospheric transmission 
Was measured on one occasion only from the aircraft at a height of gooo ft, but 
as the sky was of excellent quality on this occasion, and on all other occasion 
when observations were made from the aircraft, the results can be assumed 
to be representative [he extinction was measured by taking calibrated photo 
graphs of the rising moon In the reduction of the observations the extinction 
formula 


[= Ig" (2) 
Was used, where / 


! the observed intensity, & an absorption coefficient and m, the total mass of ai 


q 1 the intensity of radiation entering the Earth atmosphe re 


along the line of sight Lhe unit of m, is the ma or at along a vertical line of 
ight starting from the Earth’s surtac« Ihe quantity m, has been « valuated by 
numerical integration for an aircraft height of goooft and for various altitud 

above the horizon, using a standard atmosphere (13) \ plot of log / against m 
now vives the value of A \ good hit toa traight line was obtained tor value ol 
m, ranging trom g°S5 air masses to 20°8 air masse the gradient of which gives the 


value k = 0-o8e | air ma 


\ comparison of thi sbi ith tho valu given by Allen (14) for the 
various atmospheric constituents demon | most complete absenc 
of atmo pheric du t during the ( eri rvatior \tthi geographi 
yosition the contribution ’ np I. i ud Ist 


tron 
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i total k=0°077. ‘This means that the contribution from dust 1s 


For “fairly clear conditions” Allen give 
ontributior i ( c vhich 


»< 


10 2 the dust 
is about seven times greater than the 


valu 
found in the pre ent york 


We are torced by the nature of the problem to work as close to the 
po ible 0 that we should now ask the question 
extinction be measured ? 


horizon 
how accurately should the 
It is seen from equation (2) that the error AJ, in 
deriving the true intensity J, from a measured intensity / is related to the error 


\k in measurement of k, by the relation 


Ade \k | 
/. j . Og 


he] 


| plott d avainst altitude in big tor / 
| 


lhe probable error in & is of the order of 
deduce that at all altitudes greater than 

rizon the error in derivation of J, from this cause 1 
\t 10 above the horizon it is only 3 per cent lhese error 
than those due to uncertainty in the | 
itfcient to kno k to within 


' 


‘y background, and it 1 
10 pel cent Wi emphasize avain at 
ible to treat atmospheric absorption quantitatively at such 
only because of the absence of dust or 
thy to the 


hat its po 


smoke haze along the line of 
geographical po itions and partly to the use of an aircraft. 


it (Dsorptior 
It 


coefhicients refer to a point source outside the karth 
the usual practice In Inve 


xtinetion which 1 


itmosphere 
r 

tigations of the night sky light to 

one half of this (15, 16) hi procedur 

ily oOrptior | chietly duc to 


¢ a Value tor 


used becaus« 

cattering and in any direction the obse1 

cattered trom all parts of the night sky In thi 
olid angle ubtended by the zodiacal light 


ree extinction 1 used 


mae —— Se 


case it is not istified 


compat itively 





, F , 
the zodiacal light at small elongations 


position on lite photograph tin photoyraph how 


ivnitude, are suthcient data for a plate solution by a 
modification of Turner's method ( 17) in which the p! ite distortion of the optical 
ystem 1s first measured The positions of the star image m the better plates 
can be measured with a ruler to within 0-2 mm o1 1, and the first plate solution 
reproduces their positions with this accuracy he average position of thi 
horizon is reasonably well defined and the altitude of any part of the sky above 
the horizon is taken directly from the photog iph L her thus no need to have 


any KnNO le duc ol the ysition or orientation of thre j I rofl thre preel c tine 


it which the photograph was taken 


g. Axis of symmetry of zodiacal light lhe first two plates are well suited to 
this measurement because when they were taken the eclipty as vertical at the 
horizon to within a fraction of a deyvree [hese photographs were taken with 

Sun at the average distance of 1g -7 and 18 -1 below the horizon respectively 

[1 gives the deviation of the axis of symmetry of the zodiacal light from the 
er liptic averaged trom these two plat ~ All deviation ire to the north of the 


ecliptr Lhe i¢ ! y \ ot measurement 1 oft the ordet ot l I hye iverave 


viation on this occasion was only 0-17, Which in circular measure about thi 
I hese results support the view of Regener (g) that the zodiacal light 

vimmetrical with respect to the ecliptu Pre umably the axy 

tends towards the olat equator! ith decreasing distance trom the 

he change trom the ecliptic to t olar equat must occu 


to the Sun because on this o¢ isi ne ) een thre 


lhe brightne of the zodiacal light and of th hackvround ur task 
ure the brightness of the zodiacal light a function of distance trom 
i clongation ( bor thi purpose | i been used because 
is the earliest photograph of the series and hen the Sun wa 
so far below the horizon that, exce pt tora band very e to the horizon, none of 
kground brightnes is due to the cattering of sunlight in the karth 
examination of the appearance th daw! 
of sunlight is at first confine 
pat thy K\ uddenly DCcyvi rig ! i ipo 
never tim Lhe brightne Na ( mn) ure cannin 
he plat ith the Cambridge Observatories microphoton Phe procedure 
outlined by Roach ef a/ (3) ha been used to obtain ) round down te 
* Lhe re ulting data obtained trot 
miphasised at this stage that aircraft 
than ground observations for 


the brightne 
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below the horizon this is not appreciably diminished by 


from an aircraft hese observations are not In good agreement with 


those of Roach ef al. (3) lor « sample at elongation 35 Wwe obtain a brightness 
of 3-%g x10 'B vhere Bis the average brightness of the solar disk at A6300 A. 
Roach* obtained the value, on 1952 November 18/19, of 5:15 x10 *B 


l 
, 


Regener, in his careful work, obtained 7-29 x 10 '*B, using a photoelectric 
cell without a colour filter Behr and Siedentopf (2), on the other hand, at 
vave-length 5400 A, obtained the value of 3-46 10 '&-. It is important to 
decide whether these discrepancies represent a real variation or are the result of 
errors of observation lhe following points are relevant to this matter. First, 
t is most undesirable to use stars of known stellar magnitude as an intermediate 
tandard of surface brightness lhe surface brightness of the zodiacal light 
has eventually to be related to the surface brightness of the Sun, and the extinction 
produced by the Earth's atmosphere enters once into all reductions. A more 


direct and satistactory procedure is therefore to measure the brightness of the 


\/ 


zodiacal light in absolute units and to make a good estimate of the atmosphe ri 
extinction at the time of observation lhe atmospheric absorption trom a height 
of gooo ft and at an altitude of 15) above the horizon (corresponding to an elonga- 
tion of 35 ) of the order of 20 per cent at wave-length 6300 A under good 
condition it should be possible to determine the extinction coefficient with 
an error of | than 10 per cent, and hence the error in the extra-terrestrial flux 
from this cause should be less than 2 per cent Second, the error in measurement 
of brightne it elongation 35 due to an error in the estimation of the background 

mall \t thi elongation, according to these measurements the contribution 
of the zodiacal light is about 67 per cent of the total, and even if the estimate of th 
background light 1s 50 per cent in error, the error in the zodiacal light measurement 
vill he only 24 pel cent Vhe discrepancies noted above ire hue hy larger than 
this; in the present work even the sum of zodiacal light and background is still 
not equal to the zodiacal light brightness observed by Roach: the extreme 
variation in the values given by Roach is 49 per cent. In spite of the comments 


I} slue tor ¢ considerably lower than the ur i I ive-length 





No. 4, 1956 of the zodiacal light at small elongations 


of Regener, it is dithicult to believe that the observed variation 


are other than real 
‘Thirdly, as the zodiacal light is due to the scattering of sunlight, its brightness 


should be referred to that of the average Sun and not that of the centre of the 


disk. In particular, the colour of the zodiacal light should be compared with 
the colour on the average solar disk. 


~ 
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Phe data tor the brightness of the zodiacal light obtained from plate I can 
now be applied to plate I], taken a few minutes afterwards, to find the increased 
contribution of the sky brightness from scattered sunlight. ‘lhe results are 
hown in Figs. 7 and & for two depressions of the Sun below the horizon. Each 
of these diagrams gives data for various altitudes above the horizon for one 
position of the Sun below the horizon; in Fig. 3 the data relate to one altitude of 
observation and a range of positions of the Sun. All the diagrams show that 
observations of this nature near the horizon cannot be made with the Sun less 
than about 18 below the horizon ‘his is an absolute limit to observation by 
this method, which will not be appreciably improved upon by observing from 
a greate! height. 

11. Measurement of polarization. —Vhe polarization photographs on plate I] 
ire of good quality and better exposed than those on other plates; only thes« 
photographs have therefore been used in the reductions. 

he observed polarization is the sum of contributions from the zodiacal 
light, seattered sunlight and night sky background. If P, is the observed 
polarization, we have 


(Iy, + Ig + 1y)Po = InP 2+ IgPs tly Prx (4) 


here 1,, 1, and 1, are the intensities of the three components, and P,, P., Px 
their polarizations l'o obtain P, from the observed polarization it is necessary 
to know I,, be, = Pe, Px. 

Ihe separation of the observed intensity into the two parts /, and /,+ /, 
has already been discussed. A careful separation of /, and /y is not required 
because under these conditions the contribution from scattered sunlight, /¢, 
both in intensity and polarization, is far greater than that from the night sky. 
\n approximate value for /, has been taken from plate I. 

Scattered sunlight predominates in the photographs on plate IV, and this 
plate has been used for the measurement of its polarization, At an altitude of 
above the depressed horizon the observed polarization is 10-0 per cent. 
Putting in a correction for the zodiacal light contribution, this gives the valu 
of g:2 per cent for the polarization of the sunlit background. ‘This is close to the 
value of g°6 per cent, which would be predicted by a simple application of Rayleigh 
cattering theory at this elongation. ‘lhe values of P, for various altitudes 
ibove the horizon have therefore been computed from this theory. ‘This is 
prete rable to a direct measurement because at higher altitudes there is an increasing 
confusion with the zodiacal light and at lower altitudes the photographic density 
becomes too high for accurate measurement. 

Ihe value chosen for Py is not critical. Van de Hulst (18) has pointed out 
that the night sky background should be polarized very near the horizon becaus« 
here it is chiefly scattered light that is seen. Following the computations of 

in de Hulst we have taken the polarization of the sky emission to be 6 per cent 
at the horizon, diminishing to 1 per cent 10 above the horizon. 

lhe observed polarizations on plate Il vary from 7 per cent at elongation 
10 to 15 percent at elongation 30. ‘The data resulting from the reduction of the 
smoothed observations are shown in Fig. 9 \lso on Fig. g are plotted the data 
from the 1954 eclipse (§, 7) and those taken from the work of Behr and Siedentopf 


(2) he present results appear to confirm the high polarizations obtained 


yy Behr and Siedentopt, but do not agree with the smaller polarizations found 
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I44 


by Dufay* (10). It is unfortunate that insufficient exposure prevents the present 
observations being continued beyond elongation 30 , because a combination of 


these data with those of Behr and Siedentopf makes it likely that there is a 


naximum of polarization near elongation 35 


iB 
\ potent source of error in all measurements near the horizon under thes« 
nditions can result from near specular reflection trom the sea \t an angle of 
lence of 54 the reflection coefficient is about 52 per cent and the degree 


arization of the reflected light is 18 per cent. Serious error ill result if 


Iii 


any of this light intrudes upon the zodiacal light. ‘his will occur if the guidiny 
(which results in sharp star images and a blurred horizon) or the lens aberration 
ause some mixing. ‘l’o avoid the effects of diurnal motion it 1s not sufficient to 
guide well; it is also necessary to use a short exposure to give a sharp horizon 
lhe data presented here are probably not seriously affected by errors of this kind 


* Dufay makes no mention of h rrection for the sky backgrouns his observation 


undercorrected 
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‘The error of these measurements is unlikely to be greater than +3 per cent 


in absolute measure, i.e., the maximum relative error varies from +30 per cent 
at 22° elongation to +4 15 per cent at 30 elongation. 

It is now possible to interpolate the polarization measurements between 
the solar corona and the zodiacal light. ‘lable III contains polarization data 
for a smoothed model between elongations 5° and 50’, which is based partly on 
sehr and Siedentopf’s data, and is partly an interpolation between the outer 
corona and the inner zodiacal light. 

‘The polarization data of Fig. 9 and the intensity data of Fig. 6 have been used 
to calculate electron densities in the zodiacal light region along the ecliptic 
assuming that the polarization of the F component is zero for the range of 
elongation 5 to 50. ‘The results are shown plotted in Fig. 10 (middle line) 
together with data previously obtained on the outer corona. ‘These values ar 
lower than those given in a previous publication (2) because the smaller values 


for the brightness of the zodiacal light found in the present paper have been used, 





Fic. 1 Electron densities in the outer corona and inner zodiacal light regions Viiddle curve 
observed electron density Outer curves—maximum errors resulting from uncertainties in obser- 


vational data and in polarization of F component, 


At large elongations («50 ) considerable errors in electron density values 
can result from the assumption p* =o. Maximum and minimum values of 
electron density have been computed using the data of Fig. 1 for the maximum 
likely polarization of the F component (positive or negative). No account has 
been taken of the possibility of a sudden cut-off of electron density. ‘The results 
are shown in Fig. 10 as the upper and lower curves. ‘These curves also include 
the effects of an uncertainty in the brightness of the zodiacal light of 20 per cent, 
and an uncertainty in the polarization at elongation 30 of +2 per cent absolute. 

[It is important to know how symmetrical about the ecliptic is the electron 
component. The likely deviation from the ecliptic can be determined from these 





No. 4, 1956 of the zodiacal light at small elongations 379 


observations only if the shape of the electron cloud is known. ‘lo form an ordet 
of magnitude of the quantities involved we assume that the electron cloud is an 


ellipsoid of revolution with major and minor axes in the ratio of 5:1—a not 


improbable value judging from the electron density data obtained at the 1954 
eclipse (5, 7). Knowing, now, the contribution of the electrons to the total 
zodiacal light, and also the fact that the axis of the total light is within 0°-3 of the 
ecliptic (in circular measure about the Sun), we can calculate the maximum 
likely deviation of the axis of the electron component from the ecliptic. ‘The 
result is that the axis of symmetry of the electron component is probably not more 
than about 1 (in circular measure) from the ecliptic 
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RADIO OBSERVATIONS 
OF A LUNAR OCCULTATION OF ‘THE CRAB NEBULA 


H. Costain, B. Elsmore and G. R. Whitfield 


Summary 


Observations were made at wave-lengths of 3 7m and 7°9 m of the lunar 
occultation of the Crab Nebula on 1956 January 24 At 3°7 m, observations 
ere made throughout the oc ultation, and the distribution of radio brightness 


vas derived "The position of the radio source agrees closely 


icTo the nebula v 
with that of the optical nebula but its size 1s somewhat larger. 

\t 7°9m, severe man-made interference occurred during the emersion 
but observation of the immersion suggested that at this wave-length the 
ource vas considerably more extended Recent optical evidence has 

vested that the emission from the Crab Nebula ts due to the acceleration 
of high energy electrons 1n weak magnetic field a variation of the spectral 

tribution of the radio emission from centre to limb may be associated with 
riation in the energy, pectrum of the high energy’ particle 3 or ot 
field 

comparing the time of obscuration at 3°7 m with the computed value, 

an estimate w ade of the refraction occurring in the lunar atmosphere 
in electron den vy at the Moon’s surface of 10* cm was derived, a figure 

hich correspond » a surface density of the lunar atmosphere about 10 

f that of the | 


1. Introduction. During the period 1955-1956, two occultations’ of the 
Crab Nebula by the Moon were visible from Cambridge. Observations of 
these were made to determine in greate detail the distribution of radio brightness 
cross the source, previously derived by interferometric methods (Baldwin, 1954; 
Viills, 1952), and to determine with more certainty the density of the lunar 
itmeao ph re deduced from the occultation of the radio source in Gemini ( Elsmore 
ind Whitfield, 1955). 

Ihe first occultation occurred low in the sky on 1955 November 30, and the 
observations were seriously affected by man-made interference. Only the second 
occultation, that of 1956 January 24, will be discussed 

I pparatus. Observations were made at wave-lengths of 3°7 m and 7-9m. 
At both wave-lengths, phase-switching receivers (Ryle, 1952) were used in 
conjunction with fixed aerial systems. ‘The positions of the source, as seen from 
Cambridge at times of immersion and emersion, are indicated in Fig. 1. 

(a) The 3-7 m wave-length aerial system.—'Vhe aerial system was designed to 
fulfil two conflicting requirements : first to discriminate between the radia- 
tion from the source and that from the general galactic structure, and secondly to 
allow continuous observation of the source for several minutes on either side of 


the calculated times of immersion and emersion. ‘These conditions were satisfied 


by using two widely-spaced aerials on an axis which produced interference 


maxima running tangentially to the path of the source, as indicated in Fig. 1. 
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‘The two aerials were separated by a distance of 213 m ona horizontal axi 
of azimuth 120 East of North, giving an effective aperture of 37 wave -length 
Each of the two aerials consisted of a broadside array of eight half-wave dipok 

(6) The 7-gm wave-length aerial system.—At 7-9m wave-length, an existin 
aerial system was used which consisted of two groups of aerials mounted 
56 wave-lengths apart on an East-West axis. One group contained 4, the othe: 
16 Yagi aerials. ‘The maxima of the reception pattern in the direction of the 
occultation were approximately horizontal and 2 apart. ‘lwo different phasing 
arrangements of the aerials in the larger group were necessary to observe bot! 


immersion and emersion. 








Lhe path of the “ure cen from Cambridge on 1956 January 24 The axis of th 
tem has an a nt ‘ - / and the position of the maxima and minima ar 


and espectivel 


2. Observations (a) 3°7m wave-length he record of the occultation | 
shown in Fig. 2 together with a mean curve determined from records obtained 
on days before and after the occultation Besides the deflection produced by 
the Crab Nebula, the record show minor contributions from other sources 

the prin ipal one being the nebulosity IC $43 In Gemint (maximum at 
16" cgm-e U.’[,) The effect « sources may be completely removed b 
ubtracting the two curves from one another lhe difference, at any tim 
represents the radiation from the area occulted by the Moon In order to 
expr this ditference as a traction of the total intensity, it is necessary to knov 
the detlection that would have been produced at this time, by the whole sources 
lhis wa computed from a knowledge of the position of the source and the 
geometry ot the receiving system lhe occultation curves obtained at imme ion 


and emersion are shown in Fig. 3 (a) 


Ihe path of the Moon’s centre, as viewed from Cambridye pa ed very neal 


to the position of the centre of the Crab Nebula thus at immersion and 
emersion two nearly identical scans acro the source were obtained. lhe two 
curves were superposed for comparison in Fig. 3(a); the time interval for the 
best fit was found to be §9°6 minutes. ‘The calculated time interval between 


immersion and emersion was 59:2 minutes. ‘The discrepancy of 0-4 minute 
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1 
~ 





Fic. 2.—-The record of the occultation of 1956 January 24, obtained at 3-7 m The smooth curve 
is the mean of records obtained on days before and after the occultation The dashed lines show the 


calculated times of immersion and emersion of the optical centre of the Crab Nebula 


may be due to refraction in the lunar atmosphere and the 
result will be discussed in Section 4(c). 
(b) 7-gm wave-length.— 'Vhe observations at 7 


by man-made interference. ‘This limited the 


significance of this 


gm were very seriously affected 
accuracy of the measurements at 
immersion to 15-20 per cent and made observation of the emersion impossible. 
I'he immersion curve is shown in Fig. 3 (6). 
Interpretation: (a) Distribution of brightness 


During immersion, succes- 
trips of the source 


are covered by the passage of the Moon from West to 


2 60... 
* 
O50 














— 


T_T 
fey 
5 


6 


biG. 3 


1956 January 24 at (a) A= 3:7 mand (b) A 


and comersion a ine a ind ae respectit ely The time 


of wmmmersion ptica of 


7°gm. Immersion 
are relative to the calculated time 


the Crab Nebula (emersion displaced 59°6 minutes) 
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Kast. ‘lherefore the curve shown in Fig. 3(a) represents the integral of the 
distribution of brightness projected on an East-West axis. ‘This curve was 
differentiated to yield the distribution shown in Fig. 4. ‘The shape of the outer 
limbs cannot be determined accurately, but the area beneath the curve for radial 
distances greater than 2-5 minutes of arc cannot be appreciably different from 
that shown. 

‘The accuracy of the points in Fig. 3 (4) isnotsufficient to give a reliable measure 
of the strip distribution of brightness at 7-9 m; it is, however, possible to find 
the approximate width of the source by determining the occultation curves for 
Gaussian models of brightness distribution. ‘The curves for two such models, 
which correspond to probable error limits, are shown in Fig. 3(4), and the 
corresponding brightness distributions in Fig. 4. 

he most important feature brought out in Fig. 4 is the difference in apparent 
size of the radio source at the two wave-lengths. Whilst the radiation at 3-7 m 
is confined to a region just slightly larger than the visible nebula, a considerable 
portion of the radiation at 7-9 m comes from beyond this region. ‘This result 

uggests that the object extends to considerably greater radial distances than is 


revealed in photographs. Occultation measurements at a wave-length of 75 cm 


in Holland* have shown that, at this wave-length, the source is smaller than 


at 37m. Possible explanations for this progressive increase in size with 
increasing wave-length will be discussed in Section ¢ 

(6) Comparison with the optical Crab Nebula. ‘Vhe visible radiation from the 
Crab Nebula comes from a nearly elliptical region approximately centred on a 
double star lhe coordinates of the North-following component of this central 
double star are 


R.A. os 31m ayhe. Dec. 21° 58’ 54” (1950).1 


Using Bessel’s method, the times of immersion and emersion computed for this 
position were 16"27™7 and 17"26™-9 U.'T., respectively. ‘These times were 
confirmed by interpolation of the Moon’s topocentric coordinates provided by 
Mr Sadler of the Nautical Almanac Office. 

Since the position of the Moon at any time is known with great precision, 
an occultation provides a more accurate measure of the position of a radio stat 
than 1s possible with conventional methods. In this case, where the source is 
extended, the time of immersion of the centre of gravity of the East West 
strip distribution determines the effective right ascension (his time was 
16" 28™-0 UT. and corresponds to 

R.A. 05" 31™ 32%2 (1950). } 
his coincides quite closely with the position of the centre of the optical nebula 
and with the most accurate position determined from interferometric observations 
by Baldwin (1954), R.A. 05" 31™ 30%5 + 185 (1g50). A sketch of a photograph 
of the Crab Nebula, taken in red light (6300 A-6700 A), is included for comparison 
in Fig. 4 
Because the track of the Moon’s centre passed close to the position of the 


radio source, it was not possible to measure its declination accurately 


* We are idebted to Dr ¢ Seeger for nding us these resul I pul i 
t We are indebted to Dr D. W. Dewhirst of the Cambridge Observatory for makin 
letermination of this posituor 


Chis result is based on the mption that refraction occur 
onl If this is not the ca the position is R.A h 
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The high resolution of an observation of this kind makes it well suited to an 
investigation of fine structure in the brightness distribution. Unfortunately, 
the ultimate resolution, determined by Fresnel diffraction at the Moon’s limb, 
was not realized in this case due to the smoothing of the original records. ‘Uhis was 


necessary because of interference which reduced the eftective resolution to 














A. A. 
05" 31"40° 30° 
RIGHT ASCENSION (1950) 








MINUTES OF AR¢ 


hiness acy the Crab Neb la projected onto an East-West ax 
y m (dashed line The dotted line shows the position of the centre 


tion (tevo solution } 


he observed intensity, averaged in this way, showed no variations from a 
mooth curve greater than about 4 per cent of the total intensity. ‘Uhis result 
uyvyest that there are no localized region large! than having a surtace 
brightne more than double the mean brightness near the centre of the source. 

(c) The lunar atmosphere lhe ditference between the calculated and 
observed time of obscuration of the source, o”4 o"-26, may be attributed 
to refraction at the Moon’s limb. If it is assumed that the refraction occurs 
only at the sunlit limb of the Moon, an angle of refraction of 13” + 9” 1s obtained 


Although this value is hardly significant, it 1s very probable that the angle of 


retraction 1s less than 22 this figure provides an upper limit for the electron 


density at the Moon’s surface of 10% to 104 el. cm I his derivation 1s discussed 
in a separate paper (Elsmore, in preparation), where it 1s shown that the density 
of the lunar atmosphere must be less than ro '* of the density of the terrestrial 
atmosphere at sea-level, 

5. Discussion. —Earlier observations at Cambridge by Adgie have provided 
information on the spectral variations of the total emission of several sources, 


including the Crab Nebula, over the range 3°7-7-9 m. ‘These, together with the 
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present results, indicate that the emission from the central region of the nebula 


decreases with increasing wave-length, whilst in the outer regions the intensity 


Ihe latter variation ts in the same sense but 


increases rapidly with wave-length. 
~ 


appears considerably more marked than in the majority of radio sources 
a variation might be due to: (a) re-absorption in the outer shell, (6) scattering by 


in the electron density in the outet regions; since scattering increa 


Li 


irregularities 
as the square of the wave-length, an appreciable broadening of the source might 


occur at longer wave-lengths, (c) a real variation of the emitted spectrum i 


ditterent parts of the source 

Early models for the Crab Nebula (Greenstein and Minkowski 1953) proposed 
non-thermal sources embedded in ionized hydrogen, in order to account for the 
Hat spectrum in the range 30-1200 Me/s. ‘The electron densities assumed wer 


based on the supposition that the optical radiation was due to “free tre 


transitions 

‘The work of Oort and Walraven (1956) has shown clearly that both the optical 
and radio emission may be explained in terms of the acceleration of fast electron 
in weak magnetic fields, as originally proposed by Shklovsky (1953). 

The electron density of any re-absorbing shell is now quite inadequate t 


explain the observed spectrum. If, on the other hand, the observed distribution 


were due to scattering, it would be necessary to postulate a scale for the irregulas 
structure of a few km only 

It therefore seems more probable that the variation in the spectrum acro 
the source is due to a real change in the emitted spectrum, such as might aris 


from radial variations in the magnetic field or in the energy spectrum of the fast 


electrons. 
We are indebted to Mr M. Ryle and Dr A. Hewish for their helpful criticism 
‘This work was carried out as part of a programme of research at the Cavendish 


Laboratory, supported by a grant from the Department of Scientific and 


Industrial Research. One of us (C. H. C.) is indebted to the National Research 


Council of Canada for a Special Scholarship, and another (G. RK. W.) to ‘Trinity 
College, Cambridge, for a Coutt lrotter Students! ip 


C'avendish Laboratory, 
Cambridge 
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THE DISTRIBUTION OF BRIGHTNESS 
AT METRE WAVE-LENGTHS ACROSS THE SUN’S DISK 


R. G. Conway and P. A. O' Brien 
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(Received 1956 June 5) 


Summary 

A set of interferometer observations has been made near sunspot minimum 
to determine the distribution of brightness across the Sun at wave-lengths 
of 60cm, 1°4m, 3°7m and 7'9m At A6bocm, the distribution is 
completely determined, but at the longer wave-lengths the interferometer 
ibservations do not determine the distribution in the polar regions. Measure- 
ments of the eclipse of 1954 June 30 confirm the distribution at A 60 cm, and 
lead to estimates of the poleward extension at the longer wave-lengths. An 
ippendix considers the effect of the limited resolving power of the observations, 


ind shows that no error should result except in any small-scale features of 
distribution 


1. Introduction. — Several attempts have been made to obtain radio observations 
of the Sun which may be used to derive a model of the Sun’s atmosphere. Owing 
to the effects of retraction at the longer wave-lengths which originate in the 
corona and upper chromosphere, these regions may only be investigated if 
measurements are made of the distribution of brightness across the disk, as well 
as of the total emission. Observations at metre and decimetre wave-lengths 
have been made of solar eclipses by Christiansen, Yabsley and Mills (1949), 
Laffineur et al. (1950), Blum, Denisse and Steinberg (1952), Lafhineur, Michard 
et al. (1954), but none of these observations has led to a unique solution to the 
brightness distribution. Interferometer observations have been made by Stanier 
(1950) and O’Brien (1953). ‘The latter’s observations were used by O’Brien 
and Bell (1954) to derive a preliminary model of the corona. 

In this paper a new set of interferometer observations are described which 
vere made during a long period of exceptionally low sunspot activity. During 
the same period, observations were also made of a solar eclipse, and the two 
ets have been combined to give distributions of brightness which are believed 
to be more reliable than any previously published at these wave-lengths. It is 
hoped that these distributions may be used to derive a model of the Sun’s 


atmosphe re. 


2. The interferometer observations.— During the periods 1953 November 1-30 


and 1954 April-June (in neither of which was any appreciable solar activity 


reported), a series of observations was made with interferometers of variable 
pacing connected to phase switching receivers, as described by O’ Brien (1953). 
lhe wave-lengths used were 60 cm, 1-4 m, 3°7 m and 7-9 m. 

\t A6o cm the amplitude of the record was measured as a function both of 
the aerial spacing and of the interferometer position angle.* In addition, the 


* The interferometer position angle is defined as the angle between the interferometer plane of 
symmetry and the rotation axis of the Sun 
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phase of the records was measured for some of the larger spacings and was found 
to be negative ‘These results, which have been reported by O’Brien and 
‘Vandberg-Hanssen (1955), are shown in Fig. 1a. 

\t Ar4m, 3°7m and 7-gm, the amplitudes but not the phases of the 
records were measured, and only one interferometer azimuth was available, in 
an approximately kast-West direction. ‘The results of these observations are 
shown in Figs. 1b, ¢, d. ‘The phase of the subsidiary maximum at A1-4 m was 
not measured, but is assumed to be negative, as was found at A6ocm. ‘This 
agrees with measurements by Firor (1956), at A1-45 m. 

The absolute scale of brightness temperature is derived from the total 
emission of the Sun at each wave-length. At A6ocm and 1-4 m the aerial 
temperatures were calibrated against a noise signal generator, and the flux 
densities were derived from the calculated forward gain of the aerials. At 
A3°7 m and 7-9m the record amplitudes were compared with those of the 
radio star in ‘Laurus (the Crab Nebula). At each of the four wave-lengths a 
primary calibration has been made against a thermal noise source (Adgie and 
Smith, 1956). ‘he values of flux density so found are given in ‘Table I. The 
accuracy of each figure 1s estimated at 20 per cent. 


‘TABLE | 
Wave-length 60 cm 14m 
Flux Density 


in units 10 ~*° mks 200 








lhe observations were converted to brightness distributions by the method 
of O’Brien (1953). ‘The observations at A6ocm lead unambiguously to a 
brightness distribution (O’Brien and ‘’andberg-Hanssen, 1955) which is shown 
in Fig. 2. At At-g4m, 3:7 m and 7-9m the complete distribution cannot be 
obtained without a further assumption of symmetry, since the observations 
were made with a single value of the interferometer position angle, and hence 
lead only to a projected distribution, 


Previously it has been usual to assume that the distribution ts circularly 


symmetrical; in this case, the radial profile of the distribution is given by the 
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Fourier- Bessel transform of the observations (Machin, 1951). Alternatively, it 


is possible to interpret the observations in terms of a family of elliptical 


distributions. Ihe contours of brightness of each member of the family are 
issumed to be concentric ellipses with major axes in the Sun’s equatorial plane, 


, minor axis 
and with the same axial ratio B . Ditterent members of the same 
major axis 


family are distinguished by different values of 8. ‘The equatorial section of 
the brightness distribution of each member of the -same family is made the 
Fourier—Bessel transform of the observations; the circular distribution is 
therefore included in the family as the member with B= 1. 

Since the projection of a circle on to an arbitrary plane is an ellipse, all 
members of the same family have the same projection on to their major axis 
(apart from a scale factor of brightness equal to 1/8), and are therefore equivalent 
to the same set of interferometer observations. 

Brightness distributions at A1-4m, 3°7m and 7-gm have been derived 
upon this assumption of elliptical symmetry, and the equatorial sections at 
each wave-length are shown in Fig. 3. ‘The values of 8, which determine both the 


cale of brightness and the polar section, are discussed below in Section { 





Pic Relative brightness in the equator 


It is important to consider the effect on the derived distributions of the limited 
resolving power of the observations. A detailed analysis has been made, which 
follows the same methods as those of Bracewell and Roberts (1954), but which 
considers in particular the case of an isolated source, such as the Sun. ‘Thi 
analysis, which is contained in the \ppendix, shows that small-scale structure 
is lost, but that no error should occur in the average value of brightness over 
a finite interval of radius across the disk. Minimum values of this interval are 
given in the Appendix. 

3. The eclipse observations. \ partial eclipse of the Sun was visible at 
Cambridge on 1954 June 30. First contact was at 11! rom | I’. and last contact 
at 13"50" U.T. At the maximum phase, at 12" 34™ U.'T., 78 per cent of the 
Sun’s diameter was obscured. 


Ihe eclipse was observed at wave-lengths of 60 cm, 1:4 m and 3:7 m, using 
I } 3 


- 
phase-switching interferometers of small spacing, which recorded the Sun’s 
radiation throughout the period of the eclipse. A measurement of total power 
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—Eclipse curves for Norway The full lines show the observed intensities and the dashed 
lines show the expected intensities for different values of B 





No. 4, 1956 of brightness at metre wave-lengths across the Sun’s disk 391 


was also made at A3-7 m using a single large fixed aerial which observed the 


Sun in transit. ‘The time of transit was 12"17™ U.'T., close to the maximum 


phase of the ec lipse. 

The intensity of the Sun at each wave-length was measured on the day of 
eclipse and on contrel days before and after. ‘The variation of intensity during 
the eclipse, expressed as a fraction of the average intensity of the controls, 1s 
shown in Figs. 4a, 4, « \t A 3°7 m it was necessary to correct for the radiation 
from the radio-star in ‘Taurus (the Crab Nebula). ‘lhe correction was calculated 
from the known intensity and position of this source (Shakeshaft, Ryle, Baldwin, 
Elsmore and ‘Thomson, 1955) and is included in Fig. 4¢ 

Observations were also made at A6o cm and 1:5 m by Eriksen, Hauge and 
‘Tandberg-Hanssen (1955) in Norway, from a point where the eclipse was total 
The results of these observations are shown in Figs. sa, 6 

4. The combination of the eclipse and interferometer results..-At A60 cm, 
where the distribution was completely derived from the interferometer results, 
the eclipse measurements were used to check the interferometer measurements 
The expected eclipse curves were derived for the partial eclipse at Cambridg: 
and the total eclipse in Norway. ‘They are shown in Figs. 4a and §a, and are 


in good agreement with the observed curve 





Srizhine distributions (a) at the equator 


At A1-4m and 3:7 m the eclipse observations may be used to distinguish 
between different elliptical models derived interferometrically. At each wave 
length expected eclipse curves have been constructed for different values of 2 
a hese are shown in Figs. 4h, « and sb. ‘| he observed curves do not fit the 
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curves for B=1, showing that the distributions are not circularly symmetricab. 


lhe best fit is obtained at A1-4m with B=o-8, and at A3-7m with BP=o-6. 
[hese values agree with previous estimates: O’Brien (1953) obtained P=0-8 
it A1-4.m, and Blum, Denisse and Steinberg (1952) obtained 8 =0°6 at A 1-75 m. 
lhe value of 6 at A7‘gm has not been determine d, but the values at other wave 
lengths suggest that an upper limit is 0-6. On the other hand, a value less than 0-6 
vould imply very severe limb-darkening at the pole. In order to fix the scale 
of brightness at A7-g m the value of 6 will be taken as 0-6. 
he brightness distributions, in the form of polar and equatorial profiles, 
derived using these values of 8 are shown in Fig. 6. ‘These curves are suggested 
the basis for an investigation of the structure of the corona. ‘There is little 
loubt that they refer to the ‘“‘ quiet Sun” corona, since they were obtained in 
period when for many months there was no appreciable activity. Reliabl 
lues of absolute brightness have been obtained by calibrating against thermal 
tandard lhe limited resolving power will cause no error except in any small 
cale features of the distributions. ‘The assumption of elliptical symmetry means 
that the polar profiles at the longer wave-lengths are less certainly known than 
the equatorial profiles lhe equatorial profiles are almost independent of this 
umption, since the interferometer observations lead to the equatorial pro 
ection of each distribution, which virtually determines the equatorial profil 
thi cason, the most reliable deductions about the variation of the electron 
y and temperature with height in the corona may be made for the equatorial 
nons, based on the profile s for the equatorial plane. 
Ippendix his \ ppe ndix considers what distortion of the brightne 
distribution may result from the limited resolving power determined by th 
kimum aerial spacing used lhe treatment is based on that of Bracewell 
d Roberts (1954), and extends the analysis for the special case of a single isolated 
ource, such as the Sun 
Let P(r) be the principal solution for the brightness distribution, 1.e. th 
ourtier transform of the interferometer observations A(s), which are made tin 
ranye s§ Sy (kor impli ity, the Fourier transform 1s considered, so that 
P(r) 1 trictly the projected distribution. the arguments are not materially 
altered if P(r) is a radial profile, obtained by a Fourier- Bessel transform). Let 
the unknown difference between P(r) and the true distribution 7(r) be an 
invisible’”’ distribution /(r) 


P(r) A(s) cos (2msr) ds 


I(r) 2 A(s) cos (2msr) ds 


I(r)=2 ) i(s) cos (27sr) ds 


“0 


where i(s) is unknown for s > s9, identically zero for s ~~ 5». 
In the distribution considered, that of a single isolated sourc e, T(r) is limited 
to small values of r, and is of course never negative. It will now be shown that 
these restrictions increase the effective resolving power. 
It may be seen by inspection of Fig. 6, that the P(r) is negligible for r greater 
than a certain value r, at each wave-length. Since /(r) contains no Fourier 


components with s< so, it 1s either zero, or an oscillating function, changing 





t 
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ign twice in an interval not greater than 1/s,. 7'(r) is 
I(r) must be zero when r>»? It 


i 


always positive hence 


a function f(x) is limited to a range 
extension of the Convolution ‘Theorem 


‘smooth’’, 1.e 


ee in 
shows that its Fourier transform is 
does not change appreciably in an interval less than 1/x. /(r) 


limited to the range t To, and hence 1(S) 18 smooth within an 


interval of 1 Vo 
Since 1(S) 1s by definition zero for s s,, It Cannot have an appreciable value for 


§ Shy where Sy Sy + 1/T%o. 


[his result indicates that /(r), if it is not zero, is oscillatory 


changing sign 
on the average at least twice 1n an interval of 1/s 


lhe effective resolving power 
» an improvement of about twice \lthough the limited 
resolution may Cause an error in the value of P(r) at a single point, there 


no error in the average value of P(r) in an interval greater than 1/s, 


corresponds to s,, nots 


vill be 


since this is 


P 
he same as that of I(r). In other words 
han 1/S, 1 lo t but all tructurse whose 


VnhOs cale is greatet 
lhe quantities sy, 7) and 1 s, at each of the four wave-lengt!] 


tructure 
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and in Fig. 7 the principal distributions P(r) are shown with possible erro 
distributions /(r) having periods 1/s,. 


TABLI 
Wave-length 60 cm 1'4m 
168 Wave-lengths 
2'0 . 2 Sol. Radii 
06 o'8 7 Sol. Radii 


Maximum values of i(s) and /(r) may be estimated if the assumption is made 
that 7(r) is not an oscillating function with many subsidiary maxima. ‘This 
condition obtains if the maximum gradient of /(r) is less than the maximum 
gradient of P(r). ‘The functions /(r) shown in Fig. 7 are the maximum subject 
to this condition. 

Even with the somewhat implausible “ ripples’? which the distributions /(r) 
introduce, it can be seen that the large-scale features of the distributions remain 


unaltered. ‘‘ Limb-brightening”’ is certainly present at A60 cm, and is certainly 


absent at A 3°7 m. Furthermore, the region of maximum brightness at A6o cm 
certainly occurs within the limb, most probably at r<o-8 solar radii. 
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Summary 


New observations of the centre-limb variation « i I D lines 
compared with the earlier observations of Shanc 
and Shcherbina-Samoilova Central intensitic 
observers agree within the experimental error 
provide no evidence for any centre-limb variation ey 
points whose re sidual intensity is vreater th 
agrees with Allen’s conclusions but di 
other workers 
The line cores are formed by 
being formally accounted for by 


nechanisms for producing non-ze 


Introduction. —Vhe recent observations of the centre-limb variation of the 
Na D lines published by Priester (25) and Shcherbina-Samoilova (32) have 
failed to resolve the discrepancies in the earlier observations of Shane (31), 
\llen (2) and Houtgast (10) pointed out by Spitzer (33) \ new investigation 
of these lines was undertaken at Oxford in 1952 and is described in Section 1 
of this paper. ‘The results of all six investigations are compared in Section 2, 
and an attempt is made to determine the origin of the discrepancies. Finally, 
a brief theoretical discussion of the line core s, with partie ular reference to the 


central intensities, is given in Section 2. 


1. Observational method and results 
1.1. Observations.—Spectra of the D lines were obtained at seven points of 


the solar disk using the Oxford solar equipment (21, 22, 23). Each plate 


contains five spectra of equal exposure time; the « xposure time for different plate: 
varies from four to nine seconds. Of the five spectra on each plate, three are 
either exposures of the East or West limb (taken with a radial slit) or of the centre 
plus two intermediate points, and the remaining two are calibration exposure 

lhe calibration was performed by a Hilger rhodium-on-glass step wedge placed 
over the slit; the construction and calibration of the wedge have been described 
(4). Calibration exposures were made at the centre of the disk; an additional 
exposure at the centre of the disk, with the wedge slightly displaced, allows one 
to correct for non-uniformity of slit illumination (4) \ slit-width of 0-036 mm 
was used for all but the calibration exposures; this is the value for which th 


spectroscope apparatus function has been determined An orange filter (Ilford 


* 
25 
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Micro 5) was placed before the slit to reduce scattered light in the spectroscope. 


‘The plate vere brush developed in Parkhurst two-solution developer, which 
produces very little background fog. An interval of twenty-four hours or more 


as allowed to elapse between exposure and development in order to permit 


any initial fading of the latent image to occur. No experiments were performed 
on Eberhard effect; Redman’s experience (28), however, suggests that the effect 
is negligible when brush development is used. 

Of the twenty-two plates obtained, ten were selected for microphotometry 
five limb plates and five centre-plus-intermediate plates. Only one of the three 
limb spectra on eat h plate was mk rophotometered ; each of the final profiles for 
each position on the Sun is therefore the mean of five independent determinations 
from five different plat , each possessing its own characteristic curve. 

1.2. Microphotometry and reduction.—Microphotometry was performed with 
the Observatory microphotometer (20) using an effective slit width of 0-033 mm. 
‘This corre sponds to about one eight the half-width of D, - tests made on a number 
of narrow lines using both this slit and one having an effective width of o-o1ro mm 
showed that no resolution was lost by using the wider of the two. ‘The “ profile” 
tracings were made with a magnification of 48-2, both D, and D, being recorded 
on a single tracing. With this magnification the tracing is insufficiently long to 
record more than 4A on either side of the line pair, so in order to make sure of 
reaching the true continuum (in terms of which the profiles must be expressed) 
an additional tracing with a magnification of 6-89 was made for every “ profile’’ 
tracing hese “‘continuum’”’ tracings were measured at eighteen wave-lengths 
extending 22A to the violet of D, and a similar distance to the red of D,. ‘The 
wings of the D lines can be detected out to 13 \ on these tracings. ‘The estimated 
error in the level of the continuum derived from each tracing is 0°5 per cent. 
Thirty-five “ profile” tracings (corresponding to five profiles for each of seven 
positions on the di k) and an equal number of “continuum”’ tracings were 
measured lhe “ profile’ tracings were measured at eighty-two wave-lengths. 

Lhe limb pectra were microphotometered at values of cos@ (@=angle of 
emergence) of 0-199, 0°239, 0-299 and 0-380 with a scatter of 0-002. Prior to 
microphotometry the plates were marked at the corresponding calculated distances 
from the apparent limb by a dividing machine attached to a Cambridge micro- 
mete! \n independent check of these positions was made by microphotometry 
across the marks and the limb (i.e. perpendicular to the dispersion), the “true” 
limb being defined as the point of inflexion of the resulting curve. No systematic 
difference between these two methods of determining distances from the limb 
could be detected. ‘The centre and intermediate spectra were microphoto- 
metered at their centres, which corre spond to values of cos @ of 1-000, 0-800 and 
¢ 54 respectively 

Microphotometry of the calibration and uniformity spectra was performed 
in a direction perpendicular to the dispersion in two regions free from all but 
very faint lines. ‘These regions are A5877°8-A587g:2 and A5913:0-A5g14°1; 
according to the “‘ Rowland Revision” (30) they each only contain three lines of 
Rowland intensity ; and one of 2. ‘The resulting characteristic curves, 
corrected for non-uniformity of slit illumination, were combined to give a mean 
curve for each plate. 

‘lo form a mean profile the five individual profiles for each of the seven posi- 


tions on the disk were averaged logarithmically (at the same wave-lengths); this 
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has the effect of averaging errors in the continuum as wel errors at points in 
the profile. ‘he r.m.s. error of each of the mean curves is o-8 per cent 
(continuum = 100 per cent); this figure remains moré less constant over the 
profile, though there is usually an increase in the neighbourhood of blends 

1.3%. Correction for spectroscopu distortion.—Correction of the observed 
profiles for spectroscopic distortion involves the solution of the well-known 
equation of obliteration 


O(A) = | T(r’) A(A—X’)dn’, 


where O(A) is the observed profile, 7(A) the true profile, and A(A) the spectroscops 
apparatus-function. ‘lhe experimental determination of A(A) has been fully 
described elsewhere (4); because of th importance of the far wings of A(A) in 
fixing the observed central intensiti 1 fact originally pointed out by Allen and 
Redman (27,29), the determination was carried down to g x 10” of the peal 
intensity ‘The wings of the apparatus-function re smoothed and the oute 
portions were approximated by an inverse square la hich, in agreement with 
Redman’s results (27), fits the observed wings fairly well lhe advantage of 
obtaining an analytical formula for the outer parts is that A(A) can be extrapolated 
to infinity, which in principle is necessary for the solution of the above equation 

It is also necessary to know O(A) to infinity in order to obtain an exact solution 
though errors in the far distant parts naturally have little effect on the solution 
Account must therefore be taken of the numerous weak disturbing lines (mostly 
of atmospheric origin) whose effect 1 Cquliy ilent to a general reduction of the 
continuum (cf. Redman (28)). Beyond 6A from tl entres of ‘each of thi 
D lines the average level of the continuum was taken to be 96 ) per cent, a figui 


arrived at by considering the total amount of light re ( 1 by disturbing lin 


in a certain region in the neighbourhood of the J lin Within 6A, ho 


the effect was allowed for by superimposing on the measured wv rs of the D line 


rectangular profiles to represent the disturbing line he equivalent widths of 
these lines were obtained from the ‘‘ Utrecht Atla ‘ ilues wer 
multiplied by 2-1 to take account of the greater water vapour content at Oxford 
compared with that at Mt. Wilson Phe multiply is derived by 
measuring the equivalent widths of one or two selected atmosphi lines on each 
of the ten Oxford plates and comparing with t \tla hu he profiles of 
disturbing lines near the cores of the D lines were measured th the D line 
profiles, since these can have a large ettect on the solutio 

"The solution of the integral equation was carried out by urger and van Cittert 
method (§, 6), which has the advantage over other methods in that exact account 
can be taken of sharp irregularitic sin O(A) and A(A) IL he ore of the solution 
gives automatically the residuals O — { 7 1, where 7 nth approximati 
to 7’; this not only provide s a useful check on the occurrence of erro! 
putation but also gives a fairly precise indication of when the iteration 
be stopped. ‘The calculations reduce essentially to the obliteration of a su 
of curves, the first of which is the original observed cu he solutior 
first obtained for the profiles at the centre of the disk ; cot 
rapid, four iterations being necessary to reduce th 
random photometric error 0-8 per cent he final 
of the disk, together with the observed profile 


are observed points. 
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\n examination of the wings of the profiles obtained at the different values 
of cos@ shows that for r (residual intensity) greater than about 75 per cent any 
variation with cos @ is less than the random photometric error 0-8 per cent. ‘This 

urprising result is discussed in connection with the observations of the othe 
workers in Section 2. ‘The work in solving the integral equation is thereby 
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ABLE IA 
Truc profiles of dD, 


intensity in per cent of continuum 


cos 0 : 0°30 0°35 
AA (A) 
0°400 79°. 79 
350 


300 


150 
75 
200 
250 
300 
35° 
400 
450 


SN SN 
aw sJ 


“sJIsJ 


For all values of cos 6 
AA (A) I AA (A) AA (A) 


1°700 97° 300 0°600 
‘600 96° 750 0'650 
500 5 700 
400 650 


300 5 600 


oO 7Oo 
750 
100 

200 7 550 07 7oo 

100 500 ‘ $00 


000 45° goo 
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‘TABLE IB 


True profiles of D, 


intensity in per cent of continuum 


200 
300 
S00 
goo 
000 


100 
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shortened, since having solved the equation for the profiles at the centre of the 
disk, we can then apply it to the difference between the profile at the centre and 
the profiles at the other six position (his procedure shortens th labour 
since the range of integration is then finite lhe true profiles of the D lines for 
each of the seven positions are given in lable | 

1.4. Experimental error Following Redman (28) it is convenient to divide 
the sources of error into two classe s, namely those present in the observed profile 
before correction for spectroscopic distortion and those introduced by the proces 
of correction. In the first class we must consider random photometric error, 
error due to inaccuracy in the wedge calibration, and error due to scattered light 
in the spectroscope, t¢ le SCOP and atmosphe re 

the random photometri error 1s o’'d per cent (continuum 10 pel cent) 
except in the wings, where it 1s reduced to o-8 V7 3 pel ent by averaging the 
results obtained at the seven positions (see above) he rat errors of the 
wedge calibration have been given elsewhere (4) Lhe imum error in the 
line profiles due to this cause is 0-3 per cent in the region of the line flank’ 
(r~50 per cent), decreasing to 0-2 per centin the wing at the line centre the erro 
iS negligible. 

Scattered light in the spectroscope can be an important source of error in 
profile determination However, Plaskett (22) finds that for a slit height of 
19mm the scattered light in the Oxford spectroscope | only 1 per cent ol the 
total flux at a wave-length of A4320. At longer wave-lengths the amount 1 
smaller, and can be further reduced by using a smaller slit height and a colour 
filte: In the present work an orange filter (Ilford Micr ) 
exposures, Lhe slit height for the * profile - exposure wa pout mm 
since this is nearly the same as the height of the wedge, the small remaining amount 
of scattered light is nearly the same for the calibration as for the profile exposure 
and cannot therefore be a source of error. 

Scatter in the telescope and atmosphere has the effect that at ea h point on the 
image a certain percentage of the light comes from points where the profile may 
differ from that appropriate to the point under consideration he possible 
error reaches its greatest magnitude at the limb. Short-range scatter can be 
ignored owing to the slowly-varving nature of the profiles and to the compensatory 
effect of light scattered from both sides of the point in question Lhe extreme 
point in the present work (cos#=o0-2) corr spond to a distance of 1°560mm 
i.e. about 20”, from the limb. If wea ume complete compensation for a region 
of radius 2c urrounding the point then we need only consider scatter from 
distances greater than 2¢ vhich according to Kinman (11) 1s o-0o504 (surface 
brightness = 1-00) under normal conditions At the limb we must take one-half 


this value, namely o-o2 Consequently at a point for which ( 


97°5 per cent of the light originates at this point while the remaining 2°5 per cent 


comes from distances yreate! than 2¢ where the intensity at a given point in the 


profile may differ by, say, 10 per cent. Hence if r is the true intensity at some 
point in the limb profile (continuum = 100 per cent) the measured intensity ! 
rX97°5+(r—10) ¥ 2° 
10C 
the error therefore being o-2¢ per cent (Limb darkening increase this figure 
slightly). \t points in the profile where the centre-limb variation only a few 
per cent the error due to this cause can be ignored, 
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Errors arising from the process of correcting for spectroscopic distortion 
are of two kinds, namely errors caused by the actual process of solving the integral 
equation and errors due to apparatus-function errors. ‘The errors due to the 
process of solution arise out of a magnification of the random errors present 
in the original observed profiles. ‘The factor in the present work is estimated 
to be 1:8; magnification of random error does not affect the far wings of the 
profiles. 

Apparatus-function errors exert their greatest effect at the centre of the 
profile where the difference between observed and true intensity is greatest. 
It is rather difficult to perform exact calculations, but by considering the core 
and the wings of the apparatus-function separately it is possible to obtain an 
approximate value for the effect of such errors, at any rate on the central intensities. 


R.m.s. errors for the apparatus-function core have been given elsewhere (4). 


If we take 3 per cent (in units of the peak intensity) as being a representative 
figure for the error in the region —0-075A to 0-075A then we find the resulting 
error in the central intensities is 0-6 per cent. Similarly, the wings of the 
apparatus-function, which may be in error by 8-5 per cent of the stated intensity 
for points beyond 1-8A, lead to an error in the central intensities of 0-25 per cent. 
The effect of errors in intermediate regions of the apparatus-function is difficult 
to assess. However, if we take a figure intermediate to those obtained for the core 
and wings, we obtain the result that the total error in the central intensities caused 
by apparatus-function errors is 0°6+0°:25+0°4=1°3 per cent. ‘The effect of 
these errors on points in the profile other than the central intensity is probably 
very small. 

Collecting the various sources of error and assuming Gaussian distributions 
we obtain the following figures for the possible errors of the true profiles: central 
intensity, 1°g per cent; flanks, I°5 per cent, wings, 0-4 per cent. 

TABLE II 
Modern observations of the D lines 


Size of Exposure 

Observer Place Date Instrument d image time cos @ 

Houtgast Utrecht Sept., Grating ‘ia 12 cM 308-g0s oO, 
1935 60, 
May, 31, 
Aug., 14, 
1936 

Houtgast Potsdam Aug., Grating 
1938 

Allen Met. Stromlo Late Prism 
1937 

Shane Berkeley ics F.-P. and A trem 


13cm 


grating 
Priester (,éttingen Sept., Lummer 
1950 _—ipsilate and 
grating 
S.-Samoilova Moscow Aug., Grating 14cm 
1951 
Bray Oxford Sept., Prism 112A 18cm 


1953 


apparatus-function half-width.) 
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2. Comparison with previous observations.—TVable 11 gives a summary of 
modern investigations of the D lines. It will be seen that of the six observers, 
only Houtgast, Priester and the present author observed at points other than the 
centre and limb. A general idea of the agreement among the various workers 
can be obtained from Fig. 2, which shows the limb profiles of D, and D,. lt 





T 


2 4s 























Fic. 2 Limb profiles of D line 


Abscissae Angstroms 


Ordinates intensity (continuum 100 per cent) 


The cos @ values for the various observers are as follows Houtgast 
author, 0°199 Shane, 0°262 ; Allen, 0°259 shicherinna-Samonlowa, 07173 The positions 
and Rowland intensities of disturbing lines are also shown The upper diagram refers to D,, thi 
lower to Dg 


Priester and the present 
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should be noted that the values of cos 6 are not the same for all observers, the values 
being o-199 for Houtgast, Priester and the present author, 0-262 for Shane, 
0259 for Allen, and 0173 for Shcherbina-Samoilova. ‘The positions and 
Rowland intensities of disturbing lines are also shown. ‘The most noteworthy 
discrepancy, which is immediately apparent from the figure, is that Priester’s 
limb observations lie considerably higher than those of the other workers. ‘The 
discrepancies are discussed in detail below. 

2.1. Central intensities. Central intensities at centre and limb only found by 
the different observers are given in Table III. At first sight it may seem that there 
are rather large discrepancies, but if we calculate the r.m.s. deviation 
s= VY d(residual)*/5 we obtain the following values: D, centre, 1°6 per cent; 
limb, 1°6 per cent; D,-centre, 1°7 per cent; limb, o-g per cent. Since the 
estimated error in the present work is 1-9 per cent (and the values for the other 
workers are likely to be similar) we can conclude that the observations of the central 
intensities agree within the experimental error. ‘There is some evidence that the 
central intensity of dD, is greater than that of D,: the difference is comparable 


with the rom error of the mean, however, and may be due to the presence ot 


the heavy blend on the violet side of Dg. 


il 
Limb 


3 


+ 
Samoilova 1) 4 
s 


2° 63 
Mean 54 6-4 
0 values tor the limb are as follow Houtgast, Priester and Bray, o Shane 

Allen, « 9; Shcherbina-Samoilova, 0:173. Continuum intensit 0") 

Flanks and wings.—‘Vhe most striking discrepancies among the various 
workers lie in the centre-limb variation of the wings. Allen and the present 
author agree in finding little or no variation, Priester finds a very large variation, 
and the other observers find a variation lying between these extremes. At points 
where the present observations show a difference between limb and centre 
profiles of 2 per cent, Priester nevertheless finds a difference of 1 4 per cent, while 
at points where the present observations show no difference at all Priester finds 
a difference of 5 per cent to 8 per cent. Shane and Shcherbina-Samoilova also 
find a large centre-limb variation, though not as large as Priester’s. ‘These 


discrepancies are shown in ‘lable IV, which gives the difference between the 


intensity at the limb and the intensity at the centre for a number of selected 
points in the violet side of D, and the red side of D,. ‘The last line gives the 


r.m.s. deviation from the mean s= \ L(residuals)? /5, which varies from about 


1} per cent to 4) per cent. ‘The corresponding r.m.s. deviations for the actual 
intensities in the profiles, as opposed to the centre-limb differences, are similar, 


being about 1 per cent to 5 per cent. ‘The latter figures are rather large in view 


of the error estimate in the present work (1) per cent for the line flank) but this 
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estimate takes no account of possible systematic error, such as systematic errot 
in the wedge calibration It is difficult to account for the discrepancies 1n the 
centre-limb differences on this basis, however, since such systematic errors 
would presumably attect centre and limb profiles equally, leaving their difference 
unaltered 


\A (A) 


Shan 
Alle nh 
Houtgast 8-8 
Priester 14 
S.-Samoilova 6°5 


Bray 


R.m.s 


yuantity tabulated 1s limb intensit\ mnt ty i rious points mn the olet si 


D, and the red side of D The f cs quan V 2 (residuals)? /5) 


While on the subject of systematic errors, it may not be out of place to emphasize 
that in line profile work an accurate determination of the nearby continuum ts 


as Important as an accurate determination of the actual profile, since ultimately 


the profile must be expressed in terms of the continuum intensity lhe precau 


tions taken in the present work to achieve this end have been described in 
Section 1.2, but it is not certain that the other observers exercised the same care 
Shane, for example, seems to have measured the continuum at a distance of only 
1A from the line, which according to the present observations lies well within 
the wing. ‘This could lead to an error of 2 per cent or 3 per cent in the wings. 
‘The use by Houtgast and Priester of Minnaert’s method of locating the continuum 
is also open to criticism In this method it is assumed that the profile wings 
follow an inverse square law, but a test of the method applied to the present 
observations of the violet wing of dD, (which is relatively free from ble nds) shows 
that the wing is incapable of being represented by such a law. It 1s estimated 
that the use of the method may lead to errors of the order of 2 per cent or 3 per cent 
in the profile wings \ further objection to the use of Minnaert’s method 1 
that application before correction for spectroscopic distortion 1s invalid, because 
the portion of the wings and flanks used may be altered by the correction process ; 
on the other hand the true level of the continuum must be known before correction 
for spectroscopic distortion 1s carried out ‘here seems no way out of thi 
vicious circle. ‘lhe use of the method is in any case unnecessary when we are 
dealing with relatively narrow lines in a region fairly free from disturbing 
lines. 

However, it seems unlikely that this and other sources of systematic erro 
can wholly account for the discrepancies in the centre-limb variation Spitzer 


(33) mentions the possibility that the solar line profiles undergo an actual variation 
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with time: Minnaert (13) has also discussed this question and he and van der Meer 
(14, 12) have made attempts to detect such a variation. In the case of the D lines 
we can assert that if the discrepancies are due to a real variation with time, then the 
variation is uncorrelated with the solar cycle. ‘This is shown by the dates at 
which the various investigations were carried out (cf. Table II). Allen’s 
observations, which agree with the present observations in showing a very small 
centre-limb variation, were made at the solar maximum of 1937, while the present 
observations were made in 1953 near a minimum of solar activity. Priester’s 
observations, which show a large centre-limb variation, were made three years 
past the maximum of 1947. Clearly there is no correlation with the solar cycle. 
‘The possibility of a connection with active regions of the Sun cannot be ruled out, 
however. 

The possibility of a latitude effect must also be considered, though here we 
have insufficient data to come to any definite conclusions. In the present work 


the limb exposures were made in the neighbourhood of the equator. ‘Those of 


Priester and Shcherbina-Samoilova were made in the region of the poles of 
the Sun. Shane and Allen do not state the latitude of their limb observations, 
and Houtgast, after finding the deviation with latitude negligible, paid no further 
attention to the question of latitude. As far as the data go, the observations 
are consistent with a bigger centre-limb variation for the polar limb than for 
the equatorial limb—1.e. a weaker polar limb profile than equatorial limb 
profile ; 

3. A theoretical discussion of the line cores.—Vhe physical mechanism 
responsible for the formation of resonance lines 1s scattering rather than absorption. 
\n atom in the ground state absorbs a quantum which raises it to a higher state. 
\fter a time of the order of 10~* seconds the atom returns to the ground state, 
emitting a quantum of the same frequency. ‘The whole process is therefore 
equivalent to scattering. A small fraction ¢« of the excited atoms are de-excited 
by superelastic collisions with electrons, photo-ionization, or induced transitions 
to higher excited states rather than by spontaneous transitions to the ground state. 
he value of « can be determined from the measured central intensity, as we 
shall show below. 

lor the moment we assume that the scattering is coherent; the question of 
non-coherent scattering 1s considered later. ‘The discussion is confined to the 
line cores and it is assumed that », (the ratio of the line to the continuous absorption 
coefficient) is independent of 7, the optical depth in the continuous spectrum, 
and that A(r), the continuous emission per unit mass, ts a linear function of 7, 
so that 


(1) 


In the line core these are likely to be reasonably good approximations, since the 
core is formed in a thin layer at the top of the atmosphere throughout which the 
physical conditions remain fairly constant. 8 is proportional to the gradient of 
B(r) in this layer, and does not necessarily have the same value as the value given 
by limb darkening, since the continuous spectrum (observed between cos @= 1 
and cos@=0-2) originates much deeper down in the atmosphere. For this 
reason, it is probably preferable to regard 8 as an unknown parameter, at any rate 
until we have a more certain knowledge of the temperature distribution in the 
higher layers, On the other hand, should our analysis of the D lines lead to a valuc 
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very different from the limb darkening value, we should rightly regard the result 
with suspicion. In actual practice we find that when the analysis is made on the 
basis of coherent scattering the two values are very simular. 

3.1. Comparison of theory and observation.—Vhe equation of transfer for 
coherent scattering in the Eddington—Strémgren notation is 


dl {r, L) 
bv — (1 4 ny 1 At, 1) (1 ey, , (1 + en,) B(r), 


where »=cos#. When», is independent of 7, and B(r) is given by (1) the solution 
(Eddington approximation) for the emergent radiation is 


I (0, 2) 6 t ;p HAI €) 


" . . 4.1 4 I + €7),,). 
B(o) j j i % j ", + pag t \ 3( nv lw) 


(3) 
A similar formula has been used by Plaskett (19); Woolley and Stibbs (40) 
have compared (3) with Chandrasekhar’s exact solution. ‘lhe difference is very 
small except at the extreme limb (4 = 0), particularly for large »,,. 

In order to substitute for /,(0,4) the actually observed quantity 
r,(0, 2) = 1,(0, »)/1,(0, 2), we require the limb darkening in the continuous spectrum 
/,(0,). It is preferable to use the observed limb darkening for this purpose, 
rather than the limb darkening given by (3) with 7,=0, thus confining the 
assumption of linearity of B(7) to the high-lying layers in which the line cores are 
formed. Observations of limb darkening have been made by Abbot (1), Moll, 
Burger and van der Bilt (16), Raudenbusch (26), Canavaggia and Chalonge (7), 
Pierce, McMath, Goldberg and Mohler (18), and Peyturaux ( 17). By interpola 
tion from graphs the limb darkening at A 5893 was obtained for each of the observers 
and a mean was taken. ‘lhe mean limb darkening is represented within 0:5 per cent 
by the formula 

I,(o, 1) = B(o)(1 + BP’ t+y'p?), 


7. 


(4) 
with B(o) = 0-340, fp’ = 2-60, y’ 0676. ‘The quantity /,(0, 1) appearing in (3) 
is obtained by combining (4) with the measured r,(o, ,.). 

In comparing theory and observation we make use of the observations of 
Houtgast, Priester and the present author, since each of these observed at several 
intermediate points of the disk in addition to the centre and limb. 
shall consider D,, which is less influenced by blends than D,,. 
an upper limit to « 


Initially we 
We first determine 


At the line centre 7, = 4, say, and since np 


i 


Jo= V\3)"0 J(: ) 


1 we have 


Hence (3) becomes 


IL ,o(O, 2) 
Bo) 


“hed, 


+ 7,9 By. It is clear from (5) that the central intensity 
depends substantially on y/(1/7 


approximately, since 1 


lo TE), and therefore does not allow us to determine 
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, 


Ne and ‘ eparately 


letting 7, s { 


However, we can find an upper limit to €, say €max, by 
ing the mean central intensity of D, at ~=1 obtained by 


Houtgast, Prieste: and the pre sent author (r $°4 per cent) we obtain € max 0°0024. 


(When we use the mean of all six observers, 7 5°4+0°7 per cent, we obtain 
€ma 0°0035 +0°0012). ‘The actual value of € is probably fairly close to this 
upper limit, owing to the high value of no (~10* according to Harris (8)). 

We can now proceed to the solution of (3), letting « =o and € = 0-0024 in turn. 
\t each point in the profile we have two unknowns, namely 7, and f, so by using 
observed values of r (the mean of the values obtained by Houtgast, Priester and 
the present author) for two values of (1-0 and 0-2) we can obtain a pair of values 
(7,, 2) which exactly reproduce the measured r at both centre and limb at this 
pointin the profile. Ideally, values of 8 obtained in this way from different points 
in the profile should be identical. However, the profile is very insensitive to the 
value of f, particularly in the inner core; this has the effect that small variations 
in the profile give rise to quite large fluctuations in 8. ‘The best fit with observation 
occurs for B=3°7+1°0. ‘The values of ,, which, with this value of B, produce 
exact agreement with observation at «= 1 are given (both for «=o and € =0-0024) 


in ‘l'able \ (It will be noticed that for AA >o-125A the value of 7, is not affected 
by «) 


"TABLE \V 
Values of n, 


Using the values of », given in ‘Table V we can now calculate the profiles at 


other points of the disl lheory and observation are compared in Fig. 3 ‘The 


4. 


agreement is good, the theoretical curves lying well within the observational 


catter, oO we can conclude the ob ervations are well repre sented by coherent 
cattering 


he value of found above is in fair agreement with the value given by 


extrapolating the limb darkening formula (4) to the extreme limb, which gives 


I the value given by radiative equilibrium is also 2-6 


Vke calculation (35). Had we decided to select the latter value 


the start, instead of following the procedure described, our results 


| 


carcely have been different, owing to the insensitivity of the core to the 
precise value of 


both be deter 
of the t 


difference 1 


) 
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‘The values of », given in ‘Table V, when doubled, ought to reproduce correctly 
the measured profiles of D,, since H, IS proportional to the oscillator strength, 
which is twice as big for D, as for D,. When the computed and observed profiles 
of D, are compared we find quite good agreement over the inner core, but beyond 
\A = 0:20A the computed profiles fall below the observed profiles, the difference 


imounting to about 10 per cent at AA=o-40A. ‘This disagreement may be a 


consequence of our assumption that 7, is independent of 7. ‘This assumption 


- AO 4A 


O35 





O30 


O25 





0:20 


O15 


0:10 





Observed and computed centre-limb variation of D,. 
Lhscissae cos 0 


Ordinates intensity (continuum = 100 per cent) 


——Hougast Priester Bray. 
Filled circles represent values computed from (4) with «=o, B , using values of ny given in the 
second column of Table | Similar results are obtained with « 


9024, using values of Ny gwen 
in the third column of Table | 

becomes progressively le 3S and less valid as we proc eed outward from the core to 
the wings, owing to the fact that the change in effective optical depth #(=p/(1 + 7,)) 
in going from centre to limb becomes greater. For example, at AA=0-30A, 
n, = 2°55, and 7 changes from 1/3-55 = 0:28 at the centre to 0°2/3°55 =0°06 at the 
limb, a change which may be sufficient to invalidate our assumption of constant 
»,. At AA=o-10A, where », rises to 22, 7 only changes from 0-04 to o-o1 and the 


error involved in assuming 7, to be constant is corre pondingly less lhe 


2QY 
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disagreement of D, with theory indicates that the assumption is only valid for 
AA < 0°20A. 

‘The above results show that it is unnecessary to invoke non-coherent scattering 
to account for the cores of the D lines. How far is this conclusion consistent with 
our knowledge of the physical mechanisms responsible for non-coherent 
scattering? With regard to Zanstra’s mechanism (41) we may remark that in 
the higher layers of the atmosphere the collisional broadening of the ground 
tate is small owing to the low gas pressure, so that non-coherency due to this 
cause probably does not appreciably affect the line core. Deeper down collisional 
broadening becomes important and the wings may be formed non-coherently, 
as Houtgast (10) has attempted to show. (The mode of formation of the wings 
is still rather obscure, however, for Harris (8) finds that the wings can be repre- 
sented by coherent scattering, a result confirmed by Hitotuyanagi and Inaba 
(9). ‘here seems no justification for assuming pure absorption, as Priester (25) 
and other Géttingen workers have done). An additional source of non-coherency 
is re-distribution due to Doppler effect (see Unséld (36)). In so far as this 
proc affects the Doppler core, our conclusion that the cores are formed by 
coherent scattering is perhaps surprising. 

Confirmation for our conclusion can be obtained by carrying out an analysis 
of the observations on the basis of the non-coherent equation of transfer similar 
to that described above for the coherent equation. When the scattering is 
completely non-coherent/, in (2) must be replaced by./ where 


{Jil di 
[lod 


Owing to the highly peaked nature of /,, it is a good approximation to replace 
J byJ,. ‘The solution of the resulting equation of transfer is readily performed 
by the method used in deriving (3). When an attempt was made to fit this solution 


to the observations a value of £ of the order of 30 was required when 7, was put 


equal to 10° (Harris (8)). Unless independent evidence should favour such a 
high value, this can be taken as confirmation for our conclusion that the D line 
cores are formed by coherent scattering. 

3.2. The numerical value of «.—-Vhe mean of the values obtained for the 
central intensity of D, at ~=1 by all six observers is 5-4 +0°7 per cent; in the 
absence of any mechanisms for producing non-zero central intensities other 
than de-excitation of the upper state of the transition, this leads to a value of 
« of 00038 + 00012 (0°0024 when we use the results of Houtgast, Priester and the 
present author only). One such mechanism is chromospheric emission: the 
presence of chromospheric detail on D line spectroheliograms (3) shows that 
some such emission undoubtedly occurs. An additional mechanism for raising 
the central intensity is the superposition (by poor seeing or other causes) of 
yranular elements having different line-of-sight velocities; the addition of 
several slightly displaced profiles raises the central intensity (27, 28, 29). 
Such smearing is a complex effect depending on seeing, vertical apparatus- 
function half-width and microphotometer-slit height. ‘The dominating factor 
in the present work is the seeing, the seeing disk at Oxford being about 3” under 
the best conditions. ‘This figure may be compared with a vertical apparatus- 
function half-width of 1°-7 and a microphotometer-slit height of 1”. Under 
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these conditions granular elements 1” in diameter are completely smeared, and 
the central intensities correspondingly raised. 

Since part of the central intensity must be due to these two processes, 
it is clear that the above value of «, derived from the measured central intensity, 
is only an upper limit, and a direct comparison between the “ observed”’ valuc 
and a value computed from theory is not possible. Nevertheless, it is interesting 
to calculate from theory the contribution made to « by various de-exciting processes 
[tis only possible to perform this calculation for the photo-ionization and collisionai 
mechanisms; induced transitions have been qualitatively considered by Unsold 
(37), but with present theory it does not seem possible to give a quantitative result 
‘The contribution to « made by photo-ionization can be calculated with the aid 
of the formula quoted by Strémgren (34) for the photo-ionization transition 
probability of a hydrogen-like atom. If we substitute a colour temperatur 
of 6500 deg. K and a dilution factor of 0-3, following Woolley (38), we obtain 
a contribution to € of o-o00g—1.e. about one-quarter the observed value. 

he contribution made by superclastic collisions with electrons, say « 
depends on the electron pressure p, and the kinetic temperature of the electron 
lhe latter we can take to be 4650 deg. K, the boundary t mpe rature of the Sun 
according to Plaskett (24). Making use of formulae given by Woolley and 
Stibbs (39), we obtain « ie 10 -.. The latter calculation ts of som« 
interest because it gives us an upper limit to the electron pressure of 
0°26 + 0°11 dynes/cm*. Whatever the contribution to the central intensity 
made by chromospheric emission or granular motion, p, cannot exceed thi 


value or the central intensity would be higher than the measured value. 
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SPEC TROPHOTOMETRY OF ‘THE OUTER CORONA 


C. W. Allen 


(Received 1956 April 26) 


Summary 


A low dispersion spectrograph with an occulting bar was used at the solar 
eclipse of 1954 June 30 to obtain three coronal spectra extending to six solar 
radu. Photometric measurements show a slight increase of reddening with 


idial distance ‘The Fraunhofer lines of the outer corona are faint and an 


} 


nterpretation of their measured intensities leads to a rather high coronal 


ctron densit\ 


1. Observations. Ihe prime purpose of — the University of London 
Observatory Eclipse Expedition to Syd Koster, Sweden in June 1954 was to make 
photometric measurements of the corona using an occulting principle to reduce 
the great difference of brightness between the inner and outer corona. For the 
direct coronal photographs this required the use of an occulting disk already 
described in (1), and for spectra the same principle calls for an occulting bar. 
It is important to notice that the bar does not affect the colour characteristics of 


the radiation. 


ihe optica ystem for the spectrograph | hown in big. 1. \n image of the 
Sun refi ted from the coelostat mirror M is formed by the tel cope lens L on 
ot the it S hich Wa in the equatorial d ction Ihe hood H 


aperture A focussed by L onto the prism P lhe image of A was 


the aperture of the collimator-prism-camera system and 
the effective aperture of the pectrograph Lhe occultins 
ded from stray light by the hood H B is oriented at right 
lit S, and the width of B wa uch that the light from the Sun 
excluded from A; i.e. width B= width A+ BA x Sun’s diametet 
(in radians). Hence during the eclipse the extreme inner corona was just excluded 
but the rest of the corona was admitted with aperture variation partially com 
pensating for brightness variation. For the corona beyond 2°3 solar radu the 


effective aperture was complete 
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‘The camera had an f/1-0 Maksutov optical system of g:2.cm focal length and 
aperture. Unfortunately there was no opportunity of testing the camera before 
departure of the expedition and considerable spherical aberration was found to 
exist (The Maksutov corrector has since been replaced by a Schmidt plate). 
With the wide angle camera it was possible to obtain spectra of the eclipse sky 


in less than a minute. In an effort to reduce the aberrations the camera aperture 
in the dispersion direction was reduced to 3-7cm by means of A. ‘The speed of 
the camera was reduced by this aperture together with the film holder obstruction 
to 37 per cent of the full (//1) power. ‘The focal length of the telescope lens 
I, was 47°0cm, the collimator focal length 124cm, the slit width 0-25 mm, th« 
pectrographic dispersion 240 A/mm at Hy, and the whole spectrum length 7 mm 
Ilford HP film was used cut to 2cm diameter disks. 

lor the standardization of brightness and colour, test spectra were obtained 
of Sun-illuminated white blotting paper (see 1), using a greatly reduced camera 
aperture. Photographic density calibration was obtained by means of a ste} 
wedge. Standardized Fraunhofer line intensities spectra were obtained by 
making exposures during twilight. A wave-length comparison was provided 
from a mercury-helium discharge tube. ‘The various test films and the eclips 
films were developed separately but with similar conditions. 

Three coronal spectra were obtained, the exposure times being 8 sec, 24 sec, 
and 72sec. In each case the spectra could be detected to the limit of the field, 
i.c. 64 solar radi in both equatorial directions. ‘The eclipse was observed 
through variable cloud which was thickest for the shorter exposures. ‘The 
resulting images revealed a continuous spectrum with Fraunhofer absorption 
lines in the outer corona and some of the strongest Flash lines at the limb. ‘The 
coronal emission line A 5303 could be detected with difficulty close to the limb 
on the short exposure and not at all on the longer exposures. Without doubt this 
was due to the smallness of the dispersion and the faintness of the line in the 
outer corona. 

The occulting bar was useful in reducing the brightness near the limb, but 
the scale of the image (solar radius =0-:16mm) was too small in relation to the 
aberrations to allow an investigation of any steep radial photometric changes. 
‘The only measurements that could be attempted were (a) the coronal spectral 
intensity distribution in relation to radial distance, and (6) Fraunhofer line equiva- 
lent widths in relation to radial distance. ‘lhe justification for presenting the 
results is that they extend further into the outer corona than those previously 
available. 

2. Spectral intensity distribution.-{n order to determine the intensity 
distribution as a function of r, the projected distance from the Sun’s centre 
in terms of the solar radius, the spectra were examined with a Hilger Recording 
Microphotometer Runs were made both along and across the spectra. 

Density calibration curves were prepared in the usual way from the step 
wedge spectra. It was found that the (Ilford HP3) film contrast was greater 
in the green-vellow than the violet- blue regions but that the band of red sensitivity 
near H« had less contrast than all the rest of the spectrum. ‘Uhis made an apparent 
change of distribution with film density which required correction. 

lhe most critical correction to be made was that for sky illumination. At 
the limiting distance of r = 6-4 the radiation was predominantly, but not entirely, 


from the sky. ‘The only way to determine the proportion of sky light at the 
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limiting distance was to use the information from the direct photographs (1), 


(which went to r= 12 where the radiation was almost entirely sky light). A sky 
correction was deduced which left the remaining coronal spectrum a little redder 
than at the Sun’s limb. For the smaller radii the influence of the sky correction 
is much reduced but we find the resulting spectral distribution still varies rather 
slowly and consistently with the radial distance. ‘This suggests that the correction 
has not introduced any considerable errors although of course it would have been 
preferable to have made the correction entirely independently of (1). 

\fter correction the various runs were smoothed in regard to wave-length 
andr. No steep colour effects were noted, nor was there any evidence that the 
radial intensity gradient was steeper in the red and violet than in the green as 
reported by Sokolova (2). In ‘Table I the results are expressed as coronal 
brightness /. at 4000 A and 6000A relative to sunlit blotting paper brightness 
/, when both corona and paper are reduced to the same value at 5000 A ‘here 


is a slight redness of the corona relative to the paper throughout Probably 


rable | 


Log (1e/1,) where I,,= coronal brightness, I,-sunlit blotting paper brightness (reduced 
to equality at 5000 A) 





5000 A 6000 A | Index of colour 
0°03 
o'o!l 


r 
) 


0°04 
0°04 


0°06 


some of this can be explained by colouring of the paper and general reddening 
of the corona on passing through thin cloud. Of more interest 1s the apparent 
slight increase of reddening as r increases. 

In order to compare the observations with theoretical values we define an 
index of colour as the logarithmic excess of brightness at 6000 A over that at 
4000 A when compared with the distribution of the mean Sun (for theory) or of 
the sunlit blotting paper (for observations). ‘The measured values taken from 
Table I are shown in Fig. 2. ‘The calculated index of colour is zero for the K 
corona (component due to electron scatter). ‘lhe index is 0-176 for the F corona 
(diffracted component) in the outer regions (where scattering coefficient 1s 
proportional to wave-length), but in the inner corona it is rather dependent 
on the size of the diffracting particles (see Blackwell 3). In Fig. 2 we differentiate 
between (a) large diffracting particles (> 2 10 *cm radius) and (4) particles of 
10 %cm radius Ihe final colour depends on K/F the brightness ratio of the 
K to the F corona lor K/F we use (1) the van de Hulst equatorial unspot 
minimum model (4) and (i) the higher A/F values derived in Section 3 below 
lhe change of observed colour with r tends to favour a large particl iZe 
(>2x10%cm radius). Either proposed value of A, F gives agreement within 
the experimental error which might be as much as 0-05, but the results are lightly 
in favour of the higher A’ fF 
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Relation between coronal colour and radial distance 
Observed (relative to sunlit blotting paper) 
Calculated (relative to mean Sun) 
a large particle D particle radius 1o°-"* cm 
an de Hulst model (11) . fion data 


raunhofer line intensities. —lbraunhofer lines were visible in the outer 
regions of the 24sec and 72 sec exposure spectra, but they appeared to be fainte1 


and more diffuse than in the test spectra taken at twilight. ‘There does not appear 


to be any noticeable breadth difference in the emission lines of the corresponding 


comparison spectra, and the natural conclusion would be that the Fraunhofer 
lines in the outer corona had suffered an extra broadening ettect (of about 5 or 10 A). 
However this conclusion should be viewed with caution since it disagrees with 
that reached earlier by Ludendortt (5) and Allen (6) who used a much higher 
pectroscopic resolution lhe weakness of the lines and the effect of film grain 
prevent any reliable photometric test of this point. 

Line weakness and grain also prevent accurate measures of the Fraunhofer 


line intensities. Only the G-band near 4300 A was strong enough to be detected 


oronal line dept) or equivalent breadth relative 
Fraunhofer Lu 
olar radial distane 
2 per cent 
cal in de Hulst data 
/ / } trian (G) and Allen (A) 
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consistently on the microphotograms. In Fig. 3 we have plotted measurements of 
the ratio of the depth of the G-band in the corona to the depth in the twilight test 
spectra. It will be remembered that this ratio should equal (F + S8)/(K + F +S) 
where S represents the sky brightness. Some allowance has been made for the 
ditfusion of bright inner coronal light onto the outer regions. It is noticed that 
even at r= 6, where A 1s negligible, the G-band depths are less in the corona than 
in the test spectra. Probably this is related to the quality of the eclipse spectra 
(associated with the extra broadening mentioned above) and the measured ratios 
have therefore been increased by 26 per cent to give the corrected curve of Fig 3 
for which the outer corona G-band has its full depth. It is found that even when 
corrected in this way the coronal line intensities are less than expected. ‘Th 
calculated curve in Fig. 3 1s from the van de Hulst model (4) for the equator 
it spot-minimum when superimposed on the sky intensity found at this eclipse 
lig. 3 also shows earlier eclipse measurements for relatively small r made by 
Grotrian (7) and Allen (6). 

he ditference between the corrected and calculated curves of Fig. 3 is too 
large to ascribe to experimental error and the measurements suggest that the ratio 
K/F used by van de Hulst is too low. From the corrected curve of Fig. 3, and 
the coronal brightness measurements of (1), one may derive values for A and / 
for the coronal equator at the 1954 eclipse. Values obtained are compared with 
those from the van de Hulst model and from Blackwell's polarization measure 
ments (8) in ‘lable I] We tind that the new value of A, and therefore the electron 
density, in the outer corona is about } times ygreate! than the model, while 
Blackwell’s A is about 2 times greate: ‘The new estimate of /' in the inner corona 


represents a significant reduction 


mMecan 


here have been many indications from radio astronomy and solar-terrestrial 
phy ics which suggest that the electron density im of Baumbach (9) 
should not have been reduced so substantially b { . and f sceyrevation 
originally thought (10). It was from absorption line intensities and polarization 
that the original segregation was made and we no find that both of these lead 


in reasonable avreement, to higher electron det the 
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4. 1955 June 20 eclipse.-An attempt was made to repeat the observations 
at the eclipse of 1955 June 20 observed from Ceylon. In this case the spectrograph 
dispersion and sky acceptance area were both considerably greater. However 
the corona was obscured by thick cloud. 


Grants for the two eclipse expeditions mentioned were made by the Joint 
Permanent Eclipse Committee, Mr J. Miller, and the University College London. 
In Sweden the spectrograph was operated by Mr A. S. Asaad, and in Ceylon by 
Vir k. W. Foster. 


lor preliminary arrangements relating to the selection of 


both eclipse sites | have much to thank Dr H. von Kliber Much help is 
icknowledged from local authorities and residents in both Sweden and Ceylon. 


Uniwersity of London Observators 
Vill Hill Park, London, N.W.7 
19gs6 April 
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EXCITATION OF CHROMOSPHERIC EMISSION 


H. H. Plaskett 


(Received 1956 \pr 


Summary, 


The presence of the Lyman-alpha line in emission across the whole disk 
of the Sun suggests that the chromosphere is excited by a mechanism addi 
tional to that g.ven by the absorption of photospheric radiation. ‘To find 
some of the prope rties of this additional excitation the observed intensity of 
a chromospheric line, as observed at the limb of the Sun outside an eclipse, 1s 
derived as an integral containing the source function \ numerical method of 
solving this integral equation is outlined, and from the solution the excitation 
temperature and the density can be found as functions of radial height in th 
chromosphere Other applications of the integral equation are briefh 


discussed 


Introduction. —'Vhe construction of biaxial Sun-followers, notably at the 
University of Colorado, has allowed long-exposure spectra of the Sun to be 
taken from rockets. On a flight to an altitude of 81 km made in an Aerobee 
Rocket on 1952 December 12, an exposure of 25 sec showed the hydrogen 


l.yman-alpha line (La) in emission (3). ‘This result has been confirmed by 


R. ‘Tousey, and in a spectrum exhibited at the Dublin meeting of the 1.A.| 
a strong emission L« could be seen stretched more or less uniformly across the 
whole disk of the Sun. It seems probable that this is a characteristic of the 
normal chromosphere rather than the result of a fortuitous flare which would, 
presumably, have given a sharply localized emission. 

Photospheric radiation by itself appears to be incapable of producing 
I.x emission in the integrated spectrum of chromosphere and photosphere. 
Excitation of chromospheric hydrogen by photospheric radiation can occur 
in two ways—on the one hand by the absorption of La quanta from the photo 
pheric continuum giving excited atoms, and on the other by the absorption 
of quanta from the continuum to the violet of the Lyman limit giving 1onized 
atoms. In so far as the first process gives detailed balancing it will never result 
in an emission line appearing in the integrated spectrum of the chromosphere 
and photosphere. For consider a chromosphere surrounding a spherical black 
body and containing atoms capable of absorbing and emitting a single line. 
‘Then from the conservation of energy, as Professor Zanstra has pointed out,* 
ince the emission in the ingle line 1s obtained solely at the expense of the ene rvy 
from that frequency in the continuum of the radiator, the total emission in the 
line integrated over the hole chromosph« re is exactly equal to the total amount 
or energy integ! ited over the whole surface of the radiator whi h has been lost 
ir that frequency from the continuum he chromospheric line-emission 


however, 18 very né¢ arly | otropK ally distributed and that fraction, approximate ly 


* Private communication in 1927. 
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one half, which lies within the solid angle subtended by the radiator will be 
effectively lost to the frequency concerned, its only effect being slightly to warm 
the radiator. Hence if we consider radiation in the frequency of the line which 
escapes outside a sphere drawn about the chromosphere there will be a net loss 
of energy, and the spectrum will show a continuum crossed by an absorption 
line at this frequency. 

‘The second process, that of photo-ionization and recombination, at first sight 
looks more promising. In particular the photospheric radiation which is absorbed 
occurs beyond the limit of the Lyman series so that the chromospheric Lz 
resulting from recombination will not, in the integrated spectrum of photosphere 
and chromosphere, be superposed on a continuum crossed by a chromospherically 
produced absorption line at the same frequency. On the other hand the 
geometrical dilution of photospheric radiation is only one half so that the cycli 
process of ionization and recombination is but little favoured in the chromosphere 
(4). Moreover it is difficult to see how the photosphere beyond the Lyman limit 
can have sufficient energy to produce by this process an La-line of the observed 
intensity ‘The radiation which escapes from the photosphere in the far 
ultra-violet, as the observations of ‘Tousey and others show (5), 1s effectively 
black-body radiation corresponding to the boundary temperature of the photo- 
sphere. If we take this to be 4650 deg. K, then the total flux available beyond 
the Lyman limit, received outside the Earth’s atmosphere in the absence of a 
chromosphere, amounts only to 7 x 10-7 microwatts cm * a negligible fraction 
of the 5x 10 * microwatts cm * received, according to Rense (3), from Lz 
itself 


P(y, z) 
4. 


~- 


y f TO OBSERVER 


Since the absorption neither of La nor of Lyman-continuum quanta from the 
photo pli re can account for the observed La emission, we conclud that there 
must be some additional mechanism of chromospheric excitation. ‘The nature 
of this excitation is the fundamental problem posed by the chromosphere. It 1 
the purpose of the pre ent paper to show how observations of chromospheri 
lines without an eclipse, such as those of Unséld (7), Keenan (1) and notably 
the interferometer observations of ‘Treanor (6), can be used to find the excitation 
temperature of the chromosphere as a function of radial height. A knowledge 
of this temperature is clearly the first step towards finding the nature of the 
additional mode of excitation. 

1. The integral equation.—-lVhe geometrical circumstances of an observation 
of the chromosphere are shown in the figure. Consider a rectangular 
y2-coordinate system with origin on the limb of the Sun and the positive y-ax! 


y for boundary temperature of 650¢ de Ix ould the two amounts of en be comp rabl 
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pointing to the observer. ‘Then at a point P(y,z) in the chromosphere, radially 
distant r km from the boundary of the photosphere and on a sight-line which 
passes z km above the origin, we shall have 


y?=2R(r — 2). (1.1) 


‘This follows since z-_r< R, where R is the radius of the photosphere. ‘Then 
if J(v) is the specific intensity at P of a chromospheric line, the equation of 
radiative transfer in the direction of y increasing will be given by 


dI(v) 
dy «,I(v) +), 


where «, is the linear coefficient of absorption at the frequency v in the line, 
and 7, is the volume coefficient of emission. Defining the tangential optical 
depth by 


x, dy, '3.2) 


the equation of transfer takes the form 
dI(v) 
dt 


The observed intensity, /,(v), in the frequency v of the chromospheric line ; 
a height z km above the limb is then given by 


I(v) 


/ 


I(v)= | | P¥ exp t) dt, 


where 


t= | «, dy. 
is the total tangential optical depth in the frequency v of th 


: line at a height 
z km above the limb. 


If the broadened state, originating from radiative and collisional damping, 
is “naturally excited” (as may be reasonably assumed), and the atom has a 
Maxwellian velocity distribution, the linear absorption coefficient is given by 


: 
hvys 


ni(r) Bas 5 ~ H(x,v). (1.4) 


Here n,(r) is the number of atoms in the lower of the two states of the line per 
unit volume at a radial height r in the chromosphere, B,, is the Einstein transition 


coefficient, and 
h Vie" 2kT, 
< m 


where 7’, is the kinetic temperature and m is the mass of the atom. 


Unsold (8) 


Following 
H(a,v) = exp(—v?*) : (7 (1.5) 


where «=(6,+6,)/b, 6, and 6, being the half-widths of the two 
states, 


broadened 


Av /— Vy" 


b b 
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and where F(v) is a tabulated function. Similarly it may be shown*, neglecting 
radiatively stimulated downward transitions, that the volume coefficient of 
emission is given by 
hv,.° I 
= nlr)As H(a, v) (1.6) 
( = by 7 
where n,{r) is the number of atoms in the upper of the two states for the line 
per unit volume at height r, and A,, is the Einstein transition coefficient for 
spontaneous emission. Hence the source function which appears in equation (1.3) 
is given by 
O\S 
- ce nr) A, 2h(v,,°)* i 
2 12 ) exp (—hv,,°/kT,) (1.7) 
kK 477 nr) By, - - 
ince A,,/ By =g,.87h(1,,")*/g,.c°. T, 1s the required excitation temperature 
at the height r in the chromosphere. It should be noted that over the width of 
the chromospheric line the source function, j,/«,, is independent of frequency 
Now re-writ equation (1 3) in the form 
t 


L (v4.9 + bv) hs exp(—?t)dt + | hs exp(— 1?) dt, (1.8) 
k 


0 1 “f K 


first integral applies to a range of y given by o * and the second 


y y-_o, and where 
dy=tht,, (1.9) 


ince the chromosphere may be assumed to be spherically symmetrical. Now 
change the variable of integration from ¢ tor. In both integrals of equation (1.8) 
ve have from (1.2) and (1.4) 


o] 
hi 12 


By j i ss H(x, v) n,(r) dy, (1.10) 


cf \ 7 


Jy 
where we have assumed that 6, that 1s the kinetic temperature, is independent 
of radial height in the chromosphere. We have also assumed that « is independent 
of height, but since under chromospheric conditions H(z, v) = exp(—v*) this 
latter assumption is not so serious. For simplicity we further assume that 

n,(r) = n,(0) exp ( Br), (1.11) 
where f, the reciprocal of the scale height, is determined by the kinetic tempera- 
ture, the turbulent velocity, and the changing excitation and ionization. A more 
general density distribution is considered in Section 3.3. With the value for 
n(r) given in (1.11) and from (1.1), equation (1.10) takes the form 


via n,(0) exp ( a) /2R Bij H (a, v) | exp ( — &*) dé. 
< 8 IA 7 J yy 


kK 


* This derivation of the emission coefficient, as well as of the absorption coefficient, is based on 
transitions between broadened energy level It is assumed, and this 1 tatistically correct for an 
assemblage of atoms, that the probability of an atom finding itself at a distance between s and s+ ds 
(measured in trequen units) from the centre of the level is given b (8 2) d. (s 5*) where 8 is the 
half half-width of the broadened energy level (cf. L Spitzer, VN , 96, 794, 1939; H. H. Plaskett, 
V.N., 115, 256, 1955) Another way of looking at this problem has been suggested by Professor 
Zanstra to whom | am indebted for reading and commenting on this paper He points out that the 
expression (1.4) for the absorption coefficient is always correct and that the expression (1.6) for the 
volume coefficient of emission can be derived on the assumption that the emission from level 2 takes 


| 


place with complete red tribution in frequency with relative weight proportional to H(z, v) 
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For the range of values of y covered by the first integral of (1.8) we have 
+ V 2R(r—z), and consequently for this integral 


27R 
n,(o) exp ( es),/ az B H(a, v)[1 —erf V P(r —2)}. 
\Y 7 


8 2) 


< 
1 Avy, 


Following Pannekoek and Unsold (g) we put 


, hv,.° . 2n 
C(2) ~s n,(o) exp ( — Pz) P - (1.12) 


where C(z) is the total tangential optical depth at height z km above the limb 
for the centre of the line, that is when v=o and H(a,v)=1. With this notation 
the variable ¢ in the first integral of equation (1.8) is given by 


t= }C(2z)H(a,v)[1—ert V BO )| 


and where consequently 
dt 
dr 


C(2)H(a, v) P_{B(r z)}-He2 B )}. 
\/7 


Similarly for the second integral of (1.8) where ) we find 


t= }C(2z)H(a,v)|1 + ert V BU )]. (1.14) 


J 


Changing the variable of integration in both integrals of (1.8) from ¢ to s, 
and noting that the limits of integration in the two integrals are respectively ¢ 
and z, and z and ™, we find after a little algebra that 
1 (vy2° + bv) = C(z)H (a, v) exp { — 4C(2z)H(a, v)} 


— | I cosh {4C(2)H(a, v) erf V B(r—2 )} 2 exp {— ptr 
V7), % 


fh? 
2) 


This is the integral equation required in order to find the source function, / 
as a function of r from the observed intensity of the chromospheric lin 
takes an apparently simpler form if we change the variable of integration to 


SC(z)H(a, v)erf V A(r— 2), (1.16) 
when 


L(v,." + bv) =2exp{—}C(z)H(a, v)} | IY cosh {df (1.17) 


hk 
; 


If, following Pannekoek and Unsdld (9g), we assume j,/«, independent of r, 
that is of ¢, the integral in (1.17) can be explicitly evaluated and we immediately 
reproduce their result 


), 


v 


1 (149° + bv) {1 —exp | — C(z)H(a, v)}]. 

2. Numerical solution of the integral equation.—Our integral equation (1.15) 
or (1.17), contains as unknowns the constants 6, n,(o) and £8, as well as j,/«, as a 
function of r. A somewhat inelegant numerical solution can be found in the 
following way. 
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Divide the chromosphere up into four spherical shells 


}, B, 


ry, in which 
kK 


» + Bs 


where the B,; are unknown constants. ‘Then our integral equation (1.15) takes 


the form 
1 (vy.° + bv) = C(2)H(a, v) exp | — $C(z)H(a, v)} 


5 ’ 
| ae 


Ihe first three integrals (with the finite limits) in this expression can be explicitly 


valuated by changing to the variable ¢ defined in equation (1.16) 
fourth integral needs more consideration. In full it is 


Ihe 
r,)| cosh {1 C(z)H(a, v) erf V BO )! 


[B(r —=z)} '*exp {| — Bl 


BBQ 


Consider first the error function in the argument of the hyperbolic cosin 
choose rz, which is at our disposal, so that V f(r, — z,,)= 2°24, wher 


limb at which observations are made, 


\¢ 
maximum distance in km above the 


of integration 


tre 
then over the range 
11-0000 


0:9985 <erf V B(? z) 


© that with an error of less than 0-2 per cent we can write 

cosh {4C(2)H(a, v)erf V B(r z)'=cosh {}C(z)H(z, v)}, 
1 constant. For example if z,,=7250km, corresponding to 10 seconds of ar¢ 
ibove the limb, and if we take B-' ~ 150 km, a value of r,=8000km will satisfy 
these conditions. Finally changing the variable of integration from r to 


our integral (2.2) takes the form 


t 


' cosh {4C(z)H(a, v)} 
|B, B, (7rg—32)} | | ” '? exp ( n) dn t B, | 


In terms of the incomplete I’-function defined as 


J 
I*(x, p) | 7! exp ( n) dn / nH” exp ( 


0 


0 


ind tabulated by Karl Pearson (2), this may be written 


\ hs cosh {4C(2z)H(«, v)} 
x [{B, — B, P(r, Br,—2z, —})}+4B,{1 —1*(Br, 
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Hence our integral equation (2.1) may be written as a linear equation in the 
unknowns, B,, namely 
1,(v,2° + bv) = G,B, + G,B, + GpB, + GsBs. (2.4) 


From the preceding discussion the dimensionless coefhicients, G,, are readily 


found to be 
G, exp {— 4C(z)H(a, v)} sinh {4C(2)H(a, v) erf V A(r, — 2)} 
G, exp { — }C(z)H(a, v)}[sinh {4C(z)H(a, v) erf V B(r, — 2)} 
sinh {4C(2z)H(a, v) erf V P(r, — z)}} 
exp {— }C(2z)H(a, v)}[sinh {}C(2)H(a, v) erf V (rs 
sinh {}C(z)H(a, v) erf V A(r,— z)}] 
LC(z)H(a, v)[1 + exp {— C(z)H(a, v)}|[1 — 1*( 
G, = }C(z)H(a, v)[1 + exp { — C(z)H(a, v)}] 
x [41 I* (Br. z, + 4) 
If H(x,v)=exp(—v*), as is normally the case for chromospheric lines, then 
each of these coefhicients, G., may be evaluated for any and v, provided that B 
and the ratio n,(o)/b are given. In effect, therefore, the equation (2.4), linear 
in the unknowns, B., also contains in the coefhicients, G the unknown } and 
n,(o)/b. 

The solution of (2.4) to find these various unknowns depends essentially 
on the fact that over the width of the emission line the source function is 
independent of frequency, and so consequently are its specific values, 
Consider for example the observed profile at two heights such that r, 

As a result of this height restriction equation (2.4) reduces to 

1 (142° + bv) = G,B, + G4By. 

\t v=o solve the resulting pair of simultaneous equations at these two heights 
for B, and B, for each of a number of pairs of arbitrary values of 8 and n,(o)/b 
From these solutions construct double-entry tables B,-'Vable O and B,-'Vable O 
giving respectively B, and B, each as a function of the arguments f and n,(o)/d. 
Now for the same two heights and a frequency Av, = (bv), similarly solve the pair 
of simultaneous equations and construct a number of double-entry tables 
B,-'Vable 1(6) and B,-'Vable I(+) for the same arguments £ and n,(o)/b, but for 
a number of assumed values of 6. Note that for constant Av, any assumed valu 
of 6b corresponds to a new value of v. Repeat this procedure for a third frequency 
\v,, obtaining tables B,-'Table L1(+) and B.-‘lable [1(b) for the same ilues of 
the thre arguments lor each value ot b Compare ty B, lable (ph) and 
B,-'Vable 11(6), and choose that pair of arguments £ and n,(o)/b which give 
identical values of 4, in the two table gnate the resulting values of B, 
by B,{b, B,n,(o)/b inally go back to B,-Table O giving values of B, for the 


arguments f§ and n,(o)/b, and choose that pair of arguments which gives a valu 


of 2B, in tl table identical with one of the B,)b, 8,n We thus have 
found the value of 6, 8 and n,(O) b which give the re quil identical values 
of B, for the three different frequencies Av=o, Av,, Av, 

‘The same three values of 6, 8 and n,(o)/b should also give identical values 
of B, n the table iB} | ible (), B,- | able I(h) ind B. lable [I(b) (Any marked 
failure to do so would be attributable either to inadequacy of the assumed density 
distribution (1.11), or to departure from the assumed constancy of kinetic 


3° 
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temperature, or finally to the crude representation of the chromosphere by 
pherical shells. Assuming that no such failure occurs we now have B,, Bs, 
b, 8, and n,(o)/b. ‘Vo find the remaining unknowns B, and B, select two observed 
profiles /,(v,,°+6v) in the height range o< z<r,, and immediatelv solve the 


resulting two simultaneous equations for these unknowns. 


3. Applications and another form of integral equation 
3.1. Optically thin chromosphere.—\f C(z), the total tangential optical depth 
of the chromosphere at the centre of the line, is so small that its square and higher 
powers may be neglected, our integral equation (1.17) takes the form 
$C (2)H(a, 2 


I (vy2° + bv) = 2[1 — }C(z)H(a, v)| Is di 


i] K, 


‘ ‘ I ‘ v , 
C(2z)H(a, v)[1 — 4C(z)H(a, v)| J 7 '*exp(—7) dy (3.1) 
\ 7 ty 
The second step follows as a consequence of changing the variable of integration, 
C, defined in (1 .16), to7 


) defined in (2.3). If now we represent the source function 


by a polynomial in r, namely 


we find 


1 (142° + bv) = C(z)H(a, v)[1 


3) 
4) 4 
pb 


Since the g, are independent of v, a procedure similar to that outlined in the 
previous section will give the values of the six unknowns. It is possible that 
equation (3.3) may be applicable to chromospheric helium which shows little 
evidence of self-reversal. 

3.2. Interpretation of spectroheliograms.—-Once the integral equation (1.15) 
has been solved for the unknowns b, n,(o) and §, and for j,/«, as a function of r, 
we can readily predict the profile of a line at the centre of the disk. ‘The radial 


optical depth, 7, in the chromosphere at the frequency v,," + bv will be given by 


hv,.° 
K, dy = 


By —— H(a, v) | n,(r) ds 
( “by " ‘ 


hv,." I . 
2 Bij _ H(a, wv) : xp ( ir), (3.4) 
tf DON “ hy 


n,(O) 
e 


the last step following from equation (1.11). ‘The total tangential optical depth 
at the same frequency and at height o above the limb will from (1.12) be 


I 27rR 
_ H(2,0)m(0) ./ 3° 


given by 


hi 12 


C'(o)H(a, v) By. 


by 
or from (3.4) the radial optical depth will be 
C(o)H(a, leapt Br). to.23 
V 2rRB ; 


Knowing 6, n,(o) and £8 this equation immediately gives us the radial optical 
depth for any Yr. When 7 O it gives 7), the total radial optical depth of the 
chromosphere, 
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With this definition the equation of radiative transfer in the direction of + is, 


Jv 
K : 
¥ 


Hence the observed intensity of a line at the centre of the disk is 


3 


I (obs) [,(phot) exp ( T1)+ | 2 exp(—r)dr. (2.0) 
0 Ky 

Knowing /,(obs) and j,/«, as a function of r for the undisturbed chromosphere, 
we can immediately find /,(phot), that is the profile of the photospheric 
absorption line as it would appear in the absence of a chromosphere. Spectro- 
heliograms give evidence of disturbed chromospheric conditions, and having 
found /,(phot) we can use equation (3.6) to make an intelligent guess as to the 
origin and nature of these disturbances. 

3.3. Integral equation for another density law.—Our integral equation has 
been derived on the assumption that the density in the chromosphere can be 
represented by 

n,(r) =n,(0) exp( — Br). 
If this should prove inadequate, consider a law of the form 
n,(r) =C, exp (— for) + ¢, exp(— fr). 
Proceeding as in Section 1 and defining 


C hv,.° 
ol) Cy EXP ( 


c 


hv,,° 


C,(2) - c, exp (— P,2) 


where C,(z) + C,(z) is the total tangential optical depth in the centre of the line 
at height z above the limb, we readily find that 


1,(v49° + bv) 2exp | LiCo(z) + Cy (2) H(a, v)] | W cosh t di » 8) 


0 K, 


As in Section 2 this may be written as a linear equation in the unknowns, B,, 
and once again the coefficients, though more complex than those in (2.5), may 
be written down and evaluated. ‘lhe solution of these new linear equations is 
likewise possible, though clearly it will present more difficulty than the simpl 
case treated in Section 2. 

To those colleagues who have patiently listened to lectures on this topic, 
and especially to Dr P J. ‘Treanor whose pertinent questions have greatly 
clarified its exposition, | wish to express my thanks 

University Observatory 

Oxford : 
1956 April 20. 
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Summary 
The week of 1954 August 


20-2 


20-27 provided an opportunity to study the 
solar region under circumstances of unusual 
‘The calcium plage formed near the central meridian of the solar 
disk between 15400™ and 15®;0™ U.'T. on August 20. In the week following 
its formation, the calcium plage grew to at least 3500 millionths of the solar 
hemisphere and became 2°4 times as bright as the dark portions of the mottled 
background of the Ky, 3, spectroheliograms. 


development of an active 
simplicity 


Isophotometric studies of the calcium plage provide quantitative data for 


comparison with sunspot, magnetic, geomagnetic, radio flux, and coronal data 


made available by astronomers in all parts of the world. At least five flares 
or subflares were observed in the region before it traversed the west limb. 
‘The study presents a comprehensive and partially quantitative account of 
of the formation, growth and decline of an isolated centre of solar activity 
The problems encountered are illustrative of those that will be met on a 
grander scale in the solar studies of the International Geophysical Year 


It is possible to study the development of active solar regions under 
circumstances of unusual simplicity during the time when the Sun is just emerging 
from minimum in the 11-year solar cycle. 


The late summer of 1954, and especially 
the week of 


\ugust 20-27 provided exactly this opportunity. 
Between 15"00™ and 15" 10™ on August 20 a small calcium plage formed very 
close to the central meridian, at about 30° south latitude. 


By 18" 00" U.T. a spot 
was Visible within the plage. 


‘There were no other disturbed regions on the entire 
visible solar hemisphere, and this fortunate circumstance continued through 
August 25 (Plate 3). In the week following its formation, the active centre was 
rotated from central meridian to solar limb. 


During this period it grew rapidly 
in size and intensity. 


It was the site of at least five flares and of enhanced radiation 
During its passage across the west limb, distinct maxima 
in both the red and green coronal lines were observed. 
at L=175°, B 30. 


at radio frequencies. 


The region was centred 


Study of this region is based on photometric measures of calcium 
spectroheliograms taken at the McMath-Hulbert Observatory and on data 
relating to the region, most kindly provided by the directors or members of the 


staff of the following observatories and establishments in many parts of the world: 
Spec tal Data 

Cavendish Laboratory, Cambridge, England 

Cornell University, U.S.A. 

Hale Laboratory, U.S.A. 

Mount Wilson Observatory, U.S.A. 

Nagoya University, Japan 
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National Research Council of Canada 
Nederhorst Den Berg (NERA), Netherlands 
Observatoire du Pic-du-Midi, France 
Radiophysics Laboratory, Sydney, Australia 
Royal Greenwich Observatory, Herstmonceux, England 
‘Tokyo Observatory, Japan 
Published Data 
High Altitude Observatory, Climax, U.S.A. 
Journal of Geophysical Research 
Quarterly Bulletin of Solar Activity 
U.S. Naval Observatory 
Calcium plage.—'Vhe photometric study of the calcium plage consists of a 
series of isophotes made for ten-minute intervals during the early development 
of the region, then at half-hourly, hourly, and longer intervals as the region changed 
more slowly (‘Table I and Fig. 1). Intensities are expressed in units of the dark 
portions of the mottled background of the K,,. spectroheliogram Viaximum 
intensity attained by the plage was 2-4 or an increase of about one magnitude. 
‘The maximum area of the plage as outlined by the 1-2 isophote was of the order of 
3500 millionths of the solar hemisphere. 
LABLE | 
Development of calcium plage, 1954 Au 27 
= Product of Excess 
Maximum ‘Total Area* Intensity and Area, for 
Intensity Disk Hemisphere [sophotes 
2 Ip XJ x Au 


3) 


4 


1540 
1550 
1600 
1632 
1700 
1300 
Igoo 
2000 
2104 


2157 


1340 


186s 
1650 oO (1110) 
Plage l] ) \ ose to limb for stometric measurements. 


* Area is that outlined by 1.2 isophote It is reported in milli of solar disk and hemispher« 
respectively 
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AUGUST 20, 1954 


AUGUST 22, 1954 


AUGUST 24, 1954 


be I phot { calcium plage, 1 lugust 2¢ 


1.1 wophot Full lines give 1.2 
differing by o.1 from the last-dra 


», 22, 24 In general, dotted line 
isophotes and the that differ by 


n e are shown by dotted line 


5 indicate 


increments of 0.2 Viaxima 


Radiation from the calcium plage can be expressed by the summation of the 
products of excess intensity and area. 


‘This product shows a steady increase 
through August 25 if areas have been corrected for projection effects and are in 
units of the hemisphere, or a levelling off after August 
in units of the disk (‘Table I and Fig. 2). 


el 


22 if the areas are ¢ xpressed 


‘The plage was still visible on August 26 
and 27 but too close to the limb for accurate photometric measurements. 


I (Areo x Excess Intensity) for Calcium Plage 
, ———————— _ 
Ll Ayala tO* 





I Apxixio® 


2000 


Hemisphere 














——— oo _— 
204 2 , 2% 264 


August 1954 


calcium plage area and excess intensity as shown by the isophotes 
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/ 


‘The plage underwent two distinct periods of growth. ‘lhe initial development, 
lasting from August 20 to 21, was followed by a four-day period, August 22 to 25, 
during which the maximum intensity was relatively constant and distinctly 
greater than on the two preceding days. ‘lhe plage declined in both area and 
intensity after August 25. For both periods of growth, the maximum intensity 
was attained relatively quickly, and the subsequent change in the plage was 
primarily an increase in area, ‘The development of the plage is shown clearly in 
the plot in Fig. 3, where for each day intensity is shown as a function of the distance 
across the plage. For this diagram,the isophotes of irregular shape were replaced 
by circles equal to them in area. In general, the intensity profile of the mature 


plage was bell-shaped with a broad maximum. 


Intensity Variation across Plage 
(Circular Approximation) 





intensity 
2.5} 


R«Solar Radius 











050 100 R 
Distonce from Point of Maximum intensity 


ity profile for calcium pla circular approximation In these measure 


corrections have been made for projection 


‘The same type of general development is apparent in the intensity profiles 
for E-W and N-S cuts across the actually measured isophot In addition, 
the detailed profiles show that there was continuous growth of the plage in the 
E-—W direction, but that a spread in the N-S direction took place only when the 


maximum intensity was increasing (big 4). 

Sunspots.—On August 20 at about 18"00"WU.'T., at the MceMath- Hulbert 
Observatory, a small sunspot was observed to have formed within the new plage 
\ telephone call was made to the Mount Wilson Observatory tor confirmation 
and cooperation in observing the newly formed plage and spot. Dr Seth B. 
Nicholson was good enough to make special records of the spot and of the hydrogen 
and calcium plages on August 20 and to put this material at our disposal 

The Mount Wilson photograph for 13" 50"U.T. and our detailed visual 
examinations of the solar disk at 13" 20" and 14” 5 mT, all indicate that the 
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pot was not present at these times. A single spot was observed by us at Lake 
\ngelus at 1800" U.'l. A Mount Wilson photograph taken at 21" 28" U.'T. shows 
that a second member had developed At that time Dr Nicholson measured the 
magnetic field and found the spots to be of new-cycle polarity, with leader V5 
ind follower R72 he position was S30 Eo1’. By 21402"02™ the bipolar 
group was well formed and the combined area of the two components, as measured 


it Mount Wilson, was 20 millionths of the solar hemisphere. 


Intensity Variation across Plage 





ntensity 
25 


North South 


* 
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\ 
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4 
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4 ‘ 4 -_" — 
000 c ‘) izorR 000 


Distance from Point of Meximum intensity 


Ihe initial outbreak of the plage and the formation of the first spot took plac« 
in the heliographic longitude corresponding to the leading portion of the matur 
plage In general, the spots were located near the east-west extremities of the 
plage 

(reas of the spot group from August 20 to 27 as measured at the Royal 


(Greenwich Observatory, at the U.S. Naval Observatory (both Naval Obs rvatory 
and Mount Wilson plates are used), and at Lake Angelus are given in ‘lable II. 
he differences in the reported areas for this relatively simple spot group illustrate 
the difficultie inherent in making such measurement and the problems 
encountered in trying to put together data from different observatori Although 
there is some ambiguity, maximum area probably took place on August 23. If one 
interprets the observations literally, the spot group disappeared on August 26 
between 13" 56™ and 14" 50" U.'T. 

Vagnetic field.—'Vhe spot group was classified at Mount Wilson as § 20 with 
new cycle polarity Detailed magnetic measures, made by Drs Nicholson and 
tabcock at Mount Wilson and the Hale Laboratories r« spective ly, are summarized 
in ‘Table III. 

It is of interest to note that the two periods otf growth shown by the plage are 
reflected also in measurements of spot area and magnetic field In the first pe riod 
the leading spot was the large! and had the greate1 magnetic ti ld. In the second 
period the reverse was true. Although a number of the plage isophotes show 
two distinct maxima, they do not reflect clearly the ystematic change ; and 


differences suggested by spot areas and magnetic fields 
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PasLe Ii 
; don voli f é phere 


(sreenwicl ’ i Viount Wilse H Obs." 


(04-64 


) (04-34 


+O 
+o 
re) 
Vleasurements made ;¢ il Observatory 
Vieasurements of daily |: drawing 
Spots are p1 it in another part of solar disk on 1 
(d) The Servi ol | cs and Radio Astronom | Royal Observatory of 


Belg im reports that on A ; P le 


pot in thi revion W ible at ¢ 9! 3 mUT 
but had d appe ired b\ 5 g aus "| pot Wi pre ywraphed at Vilount Wilson 
Observator = 43° 56 ‘ oO 3b ‘ } | a Hiulbert Observatory at 


14" sol [ 


"TABLe II] 


f 


Summary of magnetic observations of acti 
Mount Wilson Observator 
Date Leading Spot Following Spot iborator 
Aug 204 rete) 5OOo aus 300 ius MA ! new pol irit 
21°62 1500 00 
22 1000 + 2000 
Following spot grows and develop 
field stronger than that of leader 


if 


J ~400 


2100 


200 1 50¢ 


P 

{ 
J 101 JS or Posit ele recorded as region 
| \ ipproact vest limb; record 1s 


impert 
Positive field recorded; 


does not show Wh 


Coronal emission.—'Vhere was no indication of unusual coronal radiation in 
the position of this active region when it traversed the east limb about August 13, 
a week before the appearance of the plage near the central meridian. At west 
limb passage, the excess radiation associated with the region was detected from 
\ugust 25-30. All coronal stations with observations during these days agree 
in the existence of maxima in both the red and green coronal emission lines for 
the general region of the Sun in which this active area was located. 

Coronal isophot s have been prepare d« spe ially for this study by Drs J. Rosch 
and M. ‘l'rellis from the calibrated data at the Pic-du- Midi ‘The maximum green 
line intensity was 54 in units of 10°* times the intensity at the same wave-length 
for one angstrom of the photospheric spectrum lhe red line maximum was 
22 in similar units \ccording to Dr Rosch, the ratio of the intensity of the green 
coronal emission to that of the red 1s much smaller during periods of low solar 


activity than near maximum in the solar cycle, and the ratio 54/22 should be 


considered as normal. ‘lhe Pic-du-Midi isophotes shown in Fig. 5 are projected 
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on a cone tangent to the Sun at latitude S30. ‘The isophote for the calcium 
plage for August 24 is shown in a similar projection. Although there is marked 
similarity between the plage and the region of green coronal emission, the latter 
is larger and is centred about 7 to the east, and about 2 to the south of the plage. 
These differences in position may reflect the relatively greater importance of the 
following member of the spot group as the active region approached the west limb. 
‘The size of coronal isophotes of this type is necessarily influenced by the height 
as well as by the lateral extent of the emitting region. According to Dr Résch, 
isophote “30” probably circumscribes the region of excess green coronal emission 
that should be associated with the active area, since the ‘‘ background”’ intensity 
of the green line during the minimum years falls between 15 and 25 x 10°, 
With this criterion, the calcium plage and green coronal region are strikingly 
similar in size and shape. ‘They both have an east-west extent of about 2-7’. 

According to isophotes prepared at the High Altitude Observatory, Boulder, 
Colorado, and published in Report No. HAO-25, the green coronal emission 
over the region was “ moderately bright’’ and the red emission “ very bright.” 
These statements should be interpreted as referring to the intensity range for the 
separate lines and not to an absolute scale. ‘The heliographic longitude of the 
HAO coronal isophote for the region is similar to that of the calcium plage but 
the latitude is 5° less. 

Flares.—The first flare observed in the region took place on August 22, only 
48 hours after the initial formation of the calcium plage Ihere is evidence for 
the occurrence of at least five flares and subflares before the region traversed 
the west limb. Of the four flares that we observed, that of August 23 at 13" 25™ 


was the greatest. Information relating to the flares is summarized in ‘Table LV. 


‘Tasie [V 


Flares observed in active centre, 1954 August 


x Imp Position Observe Remarks 
1650 I S30 W117 f Capri, Meudor Accompanied by high 
1 McMath velocity dark flocculus 
530 We2o MeMath 
S30} =6W3!1 MecMath Accompanied by high 
velocity dark flocculus 


W6s Kodaikanal 
Woo Me Math On calcium — spectro- 
helogram 


Prominences.—'Vhe Solar Report No. HAO-25 states for the region here studied 
that “the most definite active-region-prominences of the period (July to 
September 30, 1954) were seen at the east limb passage preceding the appearance 
of the spot group.”’ Small filaments were visible on the disk near the region. 
The report of prominences at the east limb is the earliest indication of activity 
in the area in which the plage and spot later developed 

Radio- Frequency emission: (a) decimetre wave-lengths At Nagoya, Ottawa, 
and Sydney the high resolution instruments recording at approximately 7, 10, 
and 60cm showed a region of excess flux traversing the solar disk from 


close to the central meridian on August 20-21 to August 26, 27, and 25 for the 


respective frequencies. ‘lhe maximum excess flux at these wave-lengths was 
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about 6 per cent of the flux from the quiet Sun. ‘There is some indication from 
the Sydney record that excess 60 cm radiation was being emitted from the location 
of the active area one to three days before the spot and plage became visible. 

‘The three records at decimetre wave-lengths essentially agree with each 
other as to the location of the source of the enhanced radio radiation. Further- 
more, these positions correspond to the successive daily positions of the spot and 
plage within the accuracy of the data. 

Mr Covington, at Ottawa, reports that the drift curves taken between 1954 
August 22 and 25 indicate that the 10 cm emissive region ‘‘ was very unlikely to 
have been as wide as 3 minutes of arc and can be regarded as a point source. ‘The 
antenna half power width was 8 minutes of arc”. For the spot group, the 
maximum east-west extent was about 1-5’. ‘The calcium plage was larger. ‘he 
greatest east-west chord of the 1-1 isophote of the plage corresponded to about 


‘or the more certain 1-2 isophote, the maximum east-west extent was 2-3’. 


/ 
‘The great similarity between the extent of the calcium plage and the region of 


excess green coronal emission has already been pointed out 

Although the active region here studied is not well suited for clarification of 
the detailed relationships between plage area, spot area, and excess decimetre 
radiation, it does provide further evidence that the 2800 Mc/s radiation follows 
very closely the development of the calcium plage (1) (Fig. 6). When the plage 
areas are measured in units of the disk, the agreement between plage area and 
2800 Mc/s flux is better than when the areas are measured in units of the hemisphere 
and corrected for projection. ‘The use of both area and intensity (LApy * Lpgeoes) 
as a measure of the calcium plages gives only slightly better agreement, in this 
case, than does the use of area alone. Sunspot areas diminished more rapidly 
after August 23 than did either the excess 2800 Mc/s radiation or the plage area. 
‘The high resolution 2800 Mc/s records indicate that excess emission from the 
region under consideration was observed up to and including August 27. On the 

7th the calcium plage was still visible though greatly reduced in size and intensity. 
The sunspots, however, had disappeared or ceased to be visible before 14" 50™ 
on August 26. 

It is of interest to report that, in addition to the region studied, the 
high resolution instruments at Nagoya, Ottawa, and Sydney also recorded the 
presence of another newly formed, very small, plage and spot, visible on August 26 
and 27 in the eastern part of the solar disk (Plate 3). Radiation from the small 
second region was detected and reported by the three radio-frequency observers 
even though the excess flux associated with it was only one to two per cent of the 
quiet Sun radiation at the respective frequencies. 

(b) Metre wave-lengths.—-At metre wave-lengths the excess radiation was 
associated more with flare activity than with the general development of spot 
and plage. For the time interval under consideration, there was no significant 
increase in flux until after flares had started to occur in the region. A 200 Mc/s 
noise storm with increase in flux and many superposed bursts began on August 23. 
‘The flux increased slowly from 12" 15™ U.'T tothe outbreak of a flare at 13" 25™ U.T., 
after which the flux increased rapidly for several hours and then slowly declined. 
‘The noise storm was still in progress throughout August 24. On these two 
dates the active centre was between 50 and 65 trom the centre of the solar disk 
(Plate 3 and Fig. 7). ‘The 81 Mc/s records did not show a noise storm during 
this period. 
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‘The foregoing variations in the excess flux at radio frequencies are summarized 
in ‘lable V and presented graphically in Fig. 7. ‘The 6 per cent increase in flux 
at 3750 and 2800 Mc/s associated with the August 20-27 active centre is clearly 
evident in the plot of the average daily values, but on the same scale, the comparable 
increase recorded at 500 Mc/s is scarcely detectable. ‘The ten-fold greater 
increase at 200 Mc/s illustrates the greater sensitivity of this frequency to solar 
activity than is the case for its neighbouring frequencies. 


TABLE V 
Summary of daily values of flux at radio frequencies, 1954 August 18-28 Quarterly Bulletin 
of Solar Activity 
Me/se« $1 200 460 500 545 2800 3750 
Observatory Cav. Syd Tok Ned Cor Sun NBS Syd Ned Ott Nag 


Average U.T' 12)-5 o2! off yaks 6h-s oles 18h-s o2h 12h-s 18h o4h 


Date Average Daily Flux in Umts of 10~** Watts Meter~* (C/S)! 

(a) (b) 

Aug ; 8 sa 2 +16 69(0) 

‘03 ‘ 67(0) 

27 69(0) 

14 : cea “ i 

15 d 2 22 2 2 2 7 , $4 
16 7 2 384 
20 e , 2 32 
13 ‘ 5 2 ° 2° 7 & So 
15 7 } , , 25 Y 7 81 
19 7 : ; ° 25 . SO 


(a) Data from special observations of excess flux made by Swarup and Parthasarathy with 


Christiansen’s interferometer operating at 60 cm. Quiet Sun is 16 ~ 107** 

(b) Values in parentheses give flux assigned by Covington specifically to region here 
tudied, on basis of data from 150-ft. array. Reference quiet days were August 28 and 29 
Peak value of 4°1 units 1s difference in total flux on August 23 and 29 as measured by 4-ft 


reflector. Only one region of excess 10-cm emission on disk on August 23 


(c) Additional active regions are present on the solar disk on these dat 


(c) Distinctive events.—During the period August 20-27 there were fou 
burst-type events and one noise storm onset. ‘lhe two principal bursts occurred 
at both metre and decimetre wave-lengths and were associated in time with 
known flares. Data relating to these outstanding occurrences are summarized 
in ‘Table VI. 

\t metre wave-lengths, the relationship between flares and certain well 
defined increases in flux is often masked by the use of average daily flux values 
rather than those made at closer time intervals. For this study an attempt was 
made to provide as complete a coverage at 200 Mc/s as possible. ‘lable V shows 
that there is a considerable scale difference between the several observatories 
reporting at this frequency. In order to use consecutively the hourly values 
provided by Cornell University, Nederhorst den Berg Radio (NERA), and 
‘Tokyo Observatory, the values from the latter two stations were multiplied by 
a scale factor of 1-4. ‘The hourly values in Fig. 7 show clearly the distinct rise 
in flux on August 23 that followed the largest and brightest flare observed in the 


region, 
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‘TABLE VI 
Outstanding Occurrences or Distinctive Events at Radio Freque ncives 1954 August 20-26 


(Summarized from Quarterly Bulletin of Solar Activity and Special Reports) 


Max. Inten 
Frequene Start Duration (Inst.) General Description 
Dat Vc/s (U.T.) (Min.) 10°7!'wm 


Bursts 
Aug Minor burst of None reported 
23 short duration 


23 
Burst followed by 
temporary increase one 


in flux W3I 


Minor burst No observations 


0434 weak 
0440 strong 
| 


> 


0439 


0441 3" In progress 
} 


Imp. one 


0437 22 | oo 0443-0500 
= = 4 

145 | 430 W6s 
0441 

0440 

O44i 


Noise storm 


\uy 20¢ ( (Continuum 4 quiet storm reac he 3 maxi- 


23-24 Sun, 100—soo burst/hr mum 2 hours after 


with amplitude about flare Aug 


1/2 the continuum 2391 3h 25m 


Terrestrial effects.—TVhe region here studied was the first of the new cycle to 
attain a moderately high degree of activity. Did the development of this active 
centre have any measurable effect upon the karth? In seeking an answer to 
this qui stion, we have access only to ionosphe ric and geomagnetic data. ‘To the 
best of our knowledge, ionospheric disturbances did not accompany the observed 
flares. Geomagnetically the story 1s more interesting. 

The four months, 1954 May-August, marked the minimum in the cycle of 
geomagnetic disturbances. During these months there was no disturbance 
great enough to be reported as a storm. However, during September, when the 
Sun again exhibited only the most limited activity, three geomagnetic storms took 
place (Fig. 8). Itis possible that these storms reflect only the well known increase 
in geomagnetic disturbances in March and September (2). On the other hand, 
if one seeks a solar cause, one is faced with the recrudescence of geomagnetic 
disturbances in September after the definite resurgence of solar activity in 
August. More specifically, it can be pointed out that the solar flare at CMD 65° W 
on August 26, accompanied by an outburst at radio frequencies, is the last known 


major solar event prior to the geomagnetic storm of September 1, the first such 
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storm in more than four months. ‘The attempt to explain the September 


geomagnetic disturbances in terms of preceding solar activity for the very simple 


Earth-Sun case here considered involves time relationships far from those 


generally considered acceptable in such matters. 
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Conclusions. —\t has been possible to assemble a comprehensive and partially 
quantitative account of the formation, growth, and decline of an tsolated centre 
of solar activity (‘he material has been gathered together from many sources 
‘The problems encountered are illustrative of those that will be met on a grander 
scale in the solar studies of the Inte rnational Cx ophy ical Year in 1957 So. VW here 
there are several sources of comparable data (spot areas, coronal measurements, 
200 Me/s flux) the apparent differences between reporting stations can be absorbed 
in daily ave rage provide d one does not wish to take advantage of the time pre id 
in observations for 24-hour coverage. If the latter is sought, acceptable methods 
of reconciliation for data from different instruments are needed. Where there 
is only a single measure of a solar feature its ‘‘ uncertainty’ is comfortably masked 
by lack ot information. 

Study of the active region of 1954 August 20-27, has provided additional 
evidence for the concomitant growth of plage, spot, and magnetic fields. ‘Th 
two distinct periods of development were clearly marked in these three types of 
observations. ‘The flares and excess radiation at metre wave lengths took place 


31 
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entirely within the second period in which the “ following”’ spot and magnetic 
field had become dominant and the spots reached maximum and then diminished 
in area. ‘lhe principal outburst at radio frequencies occurred just eleven hours 
before the spots are known to have disappeared completely. ‘The great simplicity 
of the solar disk during this period permitted the one dimensional sweeps of the 


high resolution radio-frequency instruments to provide data unusually free from 


ambiguities. ‘The coronal observations, though excellent, all refer to the 
declining phase of the region and should, therefore, not be compared directly 
with peak values of plage, spot, and radio flux obtained between August 22-24. 

Finally, although the sequence of events here reported is believed to be 
illustrative for active regions in general, there is evidence that the development of 
the region was unusually swift. ‘lhe time interval between the formation of the 
plage and spot was only three hours, whereas Dr Nicholson writes, and the 
authors concur, that ““ when most spots form, the K plage develops a day or more 
ahead of the spots.” ‘The complete lifetime of the spot group was slightly less 
than six day Even though the coronal measures were made during the declining 
phase of the region, the Boulder group comments: “It is unusual for the coronal 
lines to reach so strong a level so soon after the outbreak of the region.’’* The 
lifetime for the region as a whole was also very short lhe plage was not visible 
during the next rotation. ‘The only trace of the region in September was a faint 
green coronal maximum observed at the east limb. ‘The fact that the time of the 
initial outbreak of the plage is known within ten minutes, and the times of both 
formation and disappearance of the spots within three hours, is another atypical 
aspect t of the re vion, but it is on t} at the Yrowing scric ot solar obs rvations may 
Ool eliminate 
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THE RELATION BETWEEN INTERSTELLAR GAS, 
DUST AND ‘THE EMISSION FROM 
NEUTRAL HYDROGEN AT 21CM 


R. D. Davie 


11956 April ¢ 


Summary 
have been made of the 21 cm emussion from neutral hydrogen 
1 of dense interstellar dust clo Iving the galactic plan 
these directions 1 1 I har irrou ing rewion and 
to den Hi clouds at a low kineti mpet ‘These 


licate i correlation between n lan or ¢ and dust in 


cloud the density of neutral hvdrog: t ‘ that of the 
Furthermors the obse ations contirm the ) har m propos 1 by 
IKahn and Seaton to explain 


itral hydroget 


Introduction 


is of great importance in the study of the interstellar medium to obtain 
more detailed observation of the relation between interstellar ra and dust 
‘The study of this problem can now be extended by observation of the emission 


lr the most abundant 


at a wave-length of 21cm from interstellar neutral hy 
interstellar gas. It 1s also of interest to investigate the effect of other interstellar 
gases and dust upon the emission from neutral hydro 

Several of the studies of the relation between the occurrence of interstellar 
dust and 21cm emission give apparently conflicting result van de Hulst, 
Miuller and Oort (1) found a small excess of 21cm emission in the direction of 
one obscuring cloud in ‘Taurus and a deficiency in the direction of another dust 
cloud Lilley (2) also examined the ‘Taurus region and found no correlation 
between 21cm emission and the dust clouds shown by McCusky’ (3) star counts 
but found a close correlation with the absorption regions obtained from count 
of external galaxies. Heeschen (4) found that the density of neutral hydrogen 
was closely related to the colour excess of globular clusters in the vicinity of the 
galactic centre. 

It is suggested in this paper that these results appear inconsistent for either 
one or both of the following reason 

(i) ‘he optical and radio results do not refer to matter at the same distance 
in the direction studied, e. g., the star counts in ‘Taurus give information about 
nearby interstellar space, whereas the absorption of the emission from external 


galaxies measures the tota/ amount of obscuring matter in this direction: the 


’ 


neutral hydrogen emission itself comes from the fudl depth of the galaxy. 


(ii) ‘The kinetic temperature of the interstellar medium varies from cloud 


to cloud within the galaxy so that the observed radio emission is a function of 
both the density and the kinetic temperature of the neutral hydrogen. The 
observations mentioned above have been made in regions away from the galactic 


3! 


4 
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plane where the neutral hydrogen is not optically thick. ‘his paper describes 
an investigation of the relation between dust and neutral hydrogen in interstellar 
clouds lying within the galactic plane and points (i) and (i1) are elucidated. In 
order to obviate the uncertainties introduced by (i), a radio and optical study was 
made of one region in Auriga which was otherwise free of large obscuration 


features. 


2. The observational data 

2.1. Equipment.—'Vhe aerial used in this experiment had a beamwidth of 
1-6, and was a focal plane paraboloid 30 ft in diameter with a dipole feed and 
reflector. ‘his was combined with a hydrogen line spectrometer of the double- 
comparison type. Spectra and drift curves across the Galaxy in Auriga were 
obtained with a bandwidth of 40kc/s and a time constant of 30 seconds. ‘The 
survey in Cygnus was made using a bandwidth of 20 kc/s. 

2.2. The Auriga region. (a) A radio investigation.—A clearly defined region 
with 21cm emission lower than its surroundings was found in the galactic plane 
in Auriga at Dec 40-5 and R.A. =05" 06™(/= 134 °5,b=3 ). Its angular extent 
was determined from drift curves across the region with the spectrometer set 
at the frequency which showed the maximum fall in temperature at the centre 
of the region. Mean drift curves for the range Dec. = 39 to45 areshownin Fig. 1. 





— om | ; = — 
ee) h.. — 
440 OS 20 os"40" 


RIGHT ASCENSION 


1 set of drift curves across the cloud in Auriga between Dec 39 and 45 


‘The temperature scale used in this investigation is derived from a comparison 
with the peak temperature of the spectrum at /= 50° and b=o” which is taken 
as 100dey.K. After the removal of the broadening effect due to the beamwidth 
of the aerial, the region which is elongated in declination, was found to extend 
from 39 °5 to 44. Its extent to half power points is 2°-5. Its width in right 


ascension, which resulted in only a slight broadening of the beam pattern, was 


estimated to be less than 1 

Further information about the low temperature region can be obtained from 
the difference spectrum, 7',(v), in the direction of the cloud. 7',(v) is defined 
as the difference between the observed spectrum and the spectrum expected if 
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the cloud were not present, 7',(v).  7,(v) is obtained by extrapolating across 
the cloud from spectra taken on either side of it. ‘lhe mean of three sets of 
spectra taken at Dec. = 40 «5, R.A. =04" 50™, 05" 06™ and 05" 22™ are shown in 
Fig.2(a). Amean difference spectrum was obtained by combining the results for 
sets of spectra at Dec, = 40 -5 and 42 with a set taken diagonally across the region 
It was then corrected for receiver bandwidth and added to two sets of observations 
obtained at Dec 40 5 using a 20kc/s bandwidth to give the resultant difference 
spectrum shown in Fig. 2(4). ‘The accuracy of any point on this difference 
spectrum due to observational errors was calculated to be less than 1 deg. K. 
‘The centre of 7',(v) is displaced 8 km/sec relative to the frequency zero, compared 
with 15 km/sec for the centre of the background hydrogen line emission spectrum, 
T,{v), in this direction. T7)(v) has a half-width of 6km/sec compared with 


20km/sec for 7,,(v). 
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bic. 2(b).—The mean difference spectrum, 1 pv) of the Auriga cloud 
(6) Optical data.—Since there is no detailed analysis of the distribution of 
stars or obscuring matter in this part of the Milky Way, an analysis was made of 


the available information in the Henry Draper and Henry Draper Extension 
atalogues his data was extended by using photographic plates of different 
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exposure times kindly taken by Dr J. Cisar with the 18” Schmidt telescope at the 
St. Andrews University Observatory. 

An inspection of a high contrast print from a long exposure (60 minutes) 
plate revealed the existence of a clearly defined absorbing region 2° long in 
declination, and o” 7" in R.A., which coincided in position and shape with the 
radio feature described above. 

Star counts were taken across the plates at a number of limiting magnitudes 
whose values were obtained by comparing the number of stars observed per 
square degree with the number for each magnitude on the galactic equator 
given by ‘Trumpler and Weaver (5). ‘The absorption produced by the dust 
cloud was o™-25 and was calculated from a comparison of the observed number 
of stars near the cloud and the number in the direction of the cloud. Its distance 
was deduced from the magnitude at which absorption attained half its maximum 


value, and was found to be 200 parsecs, thus making its dimensions 8 parsecs by 


hom hom hom h ” 
2220 22 00 2140 2! 20 
. —- 7 = 








i = 4 = - — 
ia) 
22°20" = 22"00"_—-2i"40" ~—.21"20” 20"40" 
RIGHT ASCENSION 


Fic. 3.—A map she wt position and extent of six dust and neutral hydrogen clouds in and near 


the constellation of Cyoenu 


2parsecs. ‘his regionof the sky was found to contain no further dense obscuring 
matter out to a distance of 1 kiloparsec (the limit of the longest exposure plate), 
thus confirming the original supposition that this part of the Milky Way is relatively 
free of obscuration. It is concluded that the Auriga cloud lies between the Sun 
and the bulk of the Milky Way in this direction. 

2.3. The Cygnus region. (a) A radio investigation.—A region including the 
constellation of Cygnus was first studied by Williams (6) using a receiver bandwidth 
of 40 ke/s and taking drift curves for each 2 in declination at only one frequency 
Several low temperature regions were detected. A more complete investigation 
was therefore considered desirable. 


The area of the sky surveyed extended from 19" 00™ to 23"00™ in right 


ascension and from 25 to 5s in declination. Drift curves were obtained across 
the area at intervals of two degrees in declination with the receiver tuned to the 
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frequency of maximum emission in the galactic plane. ‘The absolute intensities 
of the drift curves were obtained by taking spectra at 2 Interv als along the galactic 
plane and calibrating them against the peak intensity of the spectrum at /= 50 
and b=o which is assumed to be toodeg. KK. Six separate low temperature 


emission regions were discovered and their positions are plotted in Fig. 3. 
secause of the sharply changing background temperature in the vicinity of these 
regions, the frequency dependence of the emission was obtained by taking transits 
across the region at 30kc/s intervals instead of taking spectra as in the Auriga 
region. ‘The difference spectra could be reconstructed from this data and 
were checked with the difference spectra obtained for regions Nos. 5 and 6 where 
the background is changing less sharply. 

(6) Optical data.—Pannekoek (7) has made a photometric survey of the 
distribution of visual surface brightness of galactic star-light in the Northern 
Milky Way. ‘This distribution shows a number of discrete absorbing clouds 
which coincide in position with the regions of low radio emission observed above. 
\n approximate estimate of their absorption can be derived from his results. 
More detailed information regarding the distance and absorption of these clouds 
was obtained from star counts made by Pannekoek (8) and Miller(g). A complete 
study of the cloud at Dec. = 31°, R.A. = 21" 00™ has been made by Franklin (10). 
\ll the optical information derived for the clouds in Cygnus is given in ‘Table I]. 


Interpretation of results 

3.1. The correlation between the distribution of dust and neutral hydrogen. 
‘The 21cm results presented above differ in nature from those already published 
in that they refer to conditions near the galactic plane where the clouds appear 
as a fall in temperature relative to the surrounding neutral hydrogen emission. 
In both Cygnus and Auriga the regions of lower neutral hydrogen emission fit 
closely to the regions of dust absorption. ‘These 21 cm phenomena can only 
be adequately* explained in terms of dense neutral hydrogen clouds which 
are at a lower kinetic temperature than their surroundings. ‘lhese results suggest, 
at least qualitatively, that neutral hydrogen and dust occur together in clouds 
near the galactic plane. 

3.2. A model interstellar cloud.—A simple model of an interstellar cloud wall 
be assumed and its effects on the emission of 21 cm neutral hydrogen radiation 
will be discussed and compared with the observations reported in this paper and 
elsewher« 

The model interstellar cloud is assumed to contain dust and the various 
interstellar gases in the same proportion as they are found in the remainder of the 
interstellar medium. Such conditions are demanded by the theoretical considera 
tions of Spitzer (11), Seaton (12) and Savedoff (13). ‘he constancy of the ratio 
of neutral hydrogen to dust is suggested by the 21 cm results pres« nted in this 
paper and is also demonstrated by the hydrogen-line observations of the galactic 
centre (4) of the anticentre (2) and of the local subsystem of the Galaxy (14). 
Supporting optical evidence: for the close relationship between gas and dust 
comes from studies of the external galaxies (15), of the Andromeda Nebula (16), 
of the Carina section of the Milky Way (17), of individual emission nebulae (18) 


* The dust particle hemselves, with dimension i to ft cm, W ot abso or scatter 
21 cm radiation Furthermore, the results cannot be explained in tert fad of neutral 
hydrogen in these regions for even if the volume of the dust cloud were comp ly devoid of neutral 
hydrogen the emission would not be sufficiently reduced. 
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and of the observed correlation between the equivalent width of interstellar lines 
and stellar reddening (19) \lthough departures from this correlation between 


the occurrence of gas and dust are found, the model assumed can be regarded as 


representative of conditions in the interstellar medium. 
The effect of an increase in the density of the model cloud upon its kinetic 
temperature will now be discussed, Spitzer (20) has shown that molecular 


hydrogen is an effective cooling agent for interstellar neutral hydrogen. ‘The 
considerations used by Kahn (21) to explain the observed kinetic temperature of 
interstellar neutral hydrogen in terms of heating due to cloud collisions and 
subsequent cooling by molecular hydrogen are now applied to determine the 
concentration of molecular hydrogen required to give interstellar clouds of 
various harmonic mean temperatures Ihe results are given in ‘Table I. ‘The 
cloud concentration factor, which has been included in the ‘Table, is the density 
of dust and gas in a cloud relative to the average cloud which has a kinetic tempera- 
ture of 125 deg. K (the observed value for interstellar neutral hydrogen). Since 
the density of molecular hydrogen, N(H 11), varies as the square of the density of 
neutral hydrogen atoms the concentration factor will be given by the square 
rootof N(H 11) relative to the average cloud. It is clear that a large increase in the 


concentration factor is required to reduce the temperature significantly below 


o deg. K 
['s 
Cloud 


Harmonic mean No. of hydrogen 


temperature molecules ‘cm concentration 
deg. K V (Hi) ag factor 


5 10 I 0°44 


\nother cooling mechanism which ts more effective at these low temperatures 
has been suggested by Seaton (22) it involves the electron excitation of 
the gaseous metal constituents of the clouds, principally singly ionized carbon. 
Again, the temperature falls with increasing concentration factor. ‘Uhus if one 
or both of these cooling mechanisms are operative the kinetic temperature of the 
model interstellar cloud will be lower, the higher the density. 

3.3. The 21 cm results compared with the cloud model.—Most of the information 
about the physical conditions in interstellar neutral hydrogen clouds can be 
derived from the difference spectrum as defined in Section 2.2(a). ‘The effect 
of a cloud lying between the Sun and the Milky Way will now be considered. 
If the brightness temperature of the neutral hydrogen emission from the Milky 
Way is 7',(v) at a frequency v, then it will be reduced on traversal of the cloud 
to 7,(v)e where 7(v) is the 21 cm optical depth of the cloud at frequency v. 


The observed emission in the precise direction of the cloud will then be 


T(v) = Ty(vyer + Tf — 7 (1) 
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where 7',’ is the kinetic temperature of the cloud. ‘lhe ditterence spectrum is 
then 


T (v) = T,(v)— Tv) 
(T,(¥)- TU (2) 


(a) The Auriga cloud. (i) The kinetic temperature._Vhe kinetic temperature 
of an interstellar cloud can be obtained from its difference spectrum. From 
equation (2), 7)(v) is zero at frequencies where 7',(v) equals the kinetic 
temperature of the cloud. Inspection of the difference spectrum for the Auriga 
cloud yields two estimates of the kinetic temperature, the average of which is 
60 deg. K 

The width of the emission spectrum from a neutral hydrogen cloud is also 
related to its kinetic temperature. If 7',(v) does not vary by more than 20 per cent, 
say, over the frequency range of emission fromthe cloud, the shape of the spectrum 
T,(v) will be similar to that for the cloud. ‘The half-width of the emission from 
the Aurigacloud is accordingly 30kc/s. ‘Temperature broadening alone would be 

+ 8ke/s for hydrogen at 6odeg.K. ‘The larger experimental value may be 
attributed to the expansion of the cloud into the surrounding space which may 
be 3 times as great as the r.m.s. velocity of the cloud atoms (23). 

The observed value of the kinetic temperature may be compared with that 
obtained using the optical results and the cloud model of section 2.2(a). ‘The 
dust absorption of o™-25 takes place in a cloud depth of 2 parsec (assuming the 
cloud to be cylindrical in shape). ‘he concentration factor can be obtained by 
a comparison with an average cloud, 4 parsec in diameter (24), and producing 
o™-1 absorption, (assuming 10 clouds per kiloparsec which on the average produce 
1™-0 of absorption). ‘he corresponding concentration factor of 5 would produce 
a kinetic temperature of 55 deg. K. ‘This is in good agreement with the value 
obtained from the 21 cm results and consequently adds weight to the assumed 
cloud model. 

(ii) The mass of the neutral hydrogen cloud.._trom the known kinetic 
temperature of the neutral hydrogen cloud and the value of 7',{v) at the centre 
frequency, the optical depth of the cloud at its centre frequency (7,) can be 
derived from equation (2). lor the Auriga cloud 7, = 1 lhe number of neutral 
hydrogen atoms per cm? in the line of sight in a cloud is given by Heeschen (4) a 


T n Je’ 
V(H) al \ 


v1 hydrogen atoms/cm? 


) 


where 7 is the halfwidth of the emission spectrum (6 x 10° cm/sec) of the cloud 
and 7'.’, its kinetic temperature, is 6bodeg. K. ‘Vhus N(H) for the Auriga cloud 
is 1°6 x 107! hydrogen atoms per cm? (250 per cm*) compared with 2-4 « 10” em 

for an average cloud 4 parsec In de pth and containing 20 atoms per cm V(H) 
is 1:2 x 10%*cm @ for all he hydrogen on the line of sight in the anticentre region 
(1). ‘The mass of the Auriga cloud, assuming the optical dimensions, is 200 solar 
masses or 4:0 x 10" grm. ‘lhe number of dust particles in a line of sight through 


the centre of a cloud producing an absorption of m magnitudes ts (2) 
V(D) = 3°25 x 10°m grains cm~* 


= %-1 x 10’ grains cm for the Auriga cloud. 
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‘The mass of the dust cloud, assuming grains 3 x 10°°cm in radius and unit 
density, is 1°3x 10grm. The ratio of the density of neutral hydrogen gas 
to dust within the cloud is therefore 300. ‘This value may be compared with 
a ratio of 100 obtained by Lilley (1) at intermediate galactic latitudes near the 
yalactic anticentré 

(iit) The distance of the neutral hydrogen cloud. ‘Vhe Doppler shift of a neutral 
hydrogen cloud relative to the Sun can give the distance of the cloud in the same 
way as the distance of spiral arms of the galaxy are obtained. ‘The velocity 
displacement for the Auriga cloud relative to the Sun is 8km/sec compared with 
15 km/sec for the spiral arm in this direction. Since the random cloud velocity 
at this distance from the centre of the galaxy is 6km/sec (25) the cloud is signi- 
heantly nearer than the bulk of the hydrogen in this direction, which is at a distance 
of 1°5 kiloparsec. A distance of the cloud based on the Oort rotation model is 
07 kiloparsec but since the random motions are appreciable this quantitative 
result will be unreliable. It is clear that both the neutral hydrogen cloud and 
the dust cloud lie between the Sun and the spiral arm in this direction. — In 
principle, the observation of the difference spectrum of a neutral hydrogen cloud 
leads to an estimate of its distance. 

(+b) The Cygnus clouds It has not been possible to study each of the six 
clouds in the Cygnus region as thoroughly as the Auriga cloud because of the 
complex spectrum in this region which has three narrow peaks; in addition, 
the clouds lie slightly away from the galactic plane where the brightness tempera- 
ture is changing rapidly \ll the information which could be obtained for each 
cloud is given by ‘lable II The position tabulated is that of the centre of the 
cloud as shown in Fig. 3. ‘The kinetic temperature is a maximum value only, 
being the actual temperature recorded at the centre frequency of the cloud 
spectrum. ‘This value would be accurate if the optical depth of the cloud were 
greater than unity. 7 is the half-width of the emission spectrum of each cloud. 


The “radio” distance was calculated in the manner described in the previous 


section. ‘The opt al data has been discussed in Section 2.3 (4). 


Tasce II 
Cloud Position Radio Results Optical Results 
Dec. R.A kinetic ” oy ‘ dust distance 
temp max absorption 


deg. K) km/s (deg, K) kp m kp 


16 
emote 
remote 
remote 
remote 

0'3 


The distances obtained from the two methods, both of which are liable to 
large errors, agree as well as may be expected. In these dense dust clouds it is 
likely that the kinetic temperature of the neutral hydrogen should be very low, 
as in fact the 21 cm results show. ‘The actual temperatures are lower than 
anticipated on the Kahn model and could be more readily explained by the presence 
of ionized carbon (Seaton’s model). 
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The cloud at De 2°, R.A. =20" 50™, may be associated with the Cygnus 


rift at about 2 kiloparsec. ‘The two distant dust clouds at Dec. = 50° exhibit only 


a small fall in temperature on the 21 cm records. ‘This is probably because the 
main body of the hydrogen emitting at their characteristic frequencies lies in 
front of and not behind them. 


4. Conclusion 

The study of the relation between dust and neutral hydrogen emission in 
\uriga and Cygnus shows clearly that regions containing excess dust are cooler 
than surrounding regions. ‘his result adds weight to the mechanisms proposed 
by Kahn and Seaton to explain the high kinetic temperature of interstellar neutral 
hydrogen. 

Equation (2) indicates that a dense neutral hydrogen cloud will appear either 
hot or cold relative to its surroundings, depending upon the values of 7',(v) and 
T,’. In the galactic plane it will nearly always appear cool as in the observations 
presented above, while away from the plane it may appear bright or cold, as in the 
results obtained by van de Hulst, Muller and Oort (1) for several regions in 
‘Taurus. 

It has been suggested that the kinetic temperature of neutral hydrogen 1s 
likely to vary throughout the Galaxy (see for example (21)). ‘The observations 
by Williams and Davies (26) of absorption of the radio sources at a wave-length 
of 21cm gave evidence for a variation from 60 deg. K to 250 deg. K in the kinetic 
temperature of the spiral arms. ‘lhe present observations show that the kinetic 
temperature of individual interstellar clouds can be definitely lower than the mean 
value. 
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Summary 


‘The radial velocities of cepheid anable stars, recently augmented by 
observations ot fifty-five southern Cepheids (10), are analysed for solar motion 
galactic rotation near the Sun, and a A-term due to a general expansion of the 
system parallel to the galactic plan Solutions are carnmed out with an 
approxim ite form for the radial velocit equation for two values of the longi 
tude /, of the galactic centre, 325 and 328 , suggested by optical and radio 
data, using a correction 4M 1*4 magnitudes to Shapley’s 1940 zero-point 
of the period-luminosity curve, an Lace distance Fy to the galactic centre 
or so kiloparsecs It is sho that there 1s no evidence in the cepheid 
motions for a A-term if J, 325 , and only slight evidence in the solution 
tor /, 325 tor a negative A f form kr cos* 6, corresponding to a 
uniform contraction of the sy e alue of ky being 
kim, sec/kpe ‘The solution for /, 32 with no A-term gives to 
0 b2°t + 1°4, We 7°3 (assumed), for t components in km 
solar motion with re pect to tl circula y of the cepheid 


Dein a 


irected towards the alactic centr nthe direction of rotation, and 

ig perpendi ilar to the galact plane Ihe alue of the Oort constant A 

tound to be 19°5 + 1°9 km/sec/kp in agreement with a determination trom a 

group of nearby Cepheids given in an earlier paper (20) on the intrinsic colours 

ot cepheid iriables "The e two determination ot ] refer to a cephe id 

distance scale in which the photographic absorption averages approximately 
magnitudes per kp¢ 


‘he rotation l ystel i unctio tf distance Ft tr the 


ilactic ¢ f t} Lime ! ) umm method (4, 11) 
Weighted mea aluc ot I ul ion (Pt, Ro Rj a(R) wl Ry } al 


obtained from the radial velocit data fe 156g cephe d variable arranged 1 
I 


i 
15 groups according to distance R, the weighted mean distance of the group 
trom the alactic centre ranging trom t} arest at 5°4 kpe to the farthest at 
11°9 kp« Ihe angular velocity a /,) of th tem at the solar distance # 
is then found from the f(2, Ry) curve and t | otions of the objec n 
valactic longitude Reliable proper motions are availal for only 158 variable 


‘hese are found to yield a value of o( Ry) which d tt differ significant! 


trom the value of A— B- 26-4 km/sec/kpe obtained by Morgan and Oort (34) 


from the proper motions of early type irs and faint stars in yweneral It 


«emphas.zed that the angular velocity of ‘ d system should | 
from epheid data alone as the On f rotation is not nece 
ame as the average for a heterogeneous collection of nearby object 


Ihe paper com lude vith a comparison between the rotation {(i i) 


from cepheid variables, and the rotation 4 b rem radio dat 
chemati models of the mass distribution in the vali tern, the t 
comprehensive being a model recently constru I hmuidt (38 I hic 
cepheid group large distances from the galactic centre appear to rotate 
more rapidly than the model, and the Zroups nearest the centre more lowls 
Ihe latter discrepancy represents a striking disagreement between radio and 
optical methods both of which purport to give an observational determination 
of the same quantity. 
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Introduction.—In 1937 Joy published his monumental work on the Radial 
Velocities of Cepheid Variable Stars (1) followed by a discussion of galactic 
rotation (2) using all the data available. Joy’s list of velocities has been further 
analysed for galactic rotation in several investigations (3-9) in which there are 
important differences in the treatment of the effect of galactic rotation on the 
radial velocities, in the method of allowance for interstellar absorption, and in 
the adopted zero-point of the period-luminosity curve. Valuable though these 


analyses have been, they have all suffered to some extent from the inadequate 


distribution of the observations in galactic longitude due to the absence of radial 
velocity data within a large sector inaccessible from the northern hemisphere. 
‘This deficiency has been remedied by the recent observations of 55 variables 
presented in the communication Radial Velocities of Cepheid Variable Stars in 
the Southern Hemisphere (10) which completes an extension of Joy’s programm« 
to southern Cepheids with photographic magnitudes brighter than 10-5 at 
minimum light. ‘These observations increase the number of cepheid velociti 

to 195, including some objects belonging to Population Il. A discussion of 


tl 


¢ augmented list of velocities is given in this communication 


1. Lhe scope of the work 
Lhe object of this analysis of the data is to determine (i) the Oort constant 
of differential rotation near the Sun, (it) the solar velocity with r spect to the 
rotating system of Cepheids, (i) the level of significance of the A-term, and 
(iv) the rotational structure of the system using Camm ’s treatment of the radial 
velocity and proper motion equations (4, 11). ‘The data required for such an 
inalysis are radial velocities and proper motions, absolute magnitudes from the 
period luminosity law, and distance moduli corrected for interstellar absorption 
1.1. The period luminosity curve ‘The correlation between the periods of 
lactic Cepheids and their luminosities may be taken from the periods and 
apparent magnitudes of cepheid variables in the Small Magellanic Cloud. 
\ccording to Baade, the period-luminosity curves for type I Cepheids (classical 
Cephei variables) and type Il Cepheids (RR Lyrae and W Virginis stars) differ 
in zero-point by 1°5 photographic magnitudes. ‘lhe period-luminosity curv: 
given by Shapley (12) in 1930 was obtained by adjusting the zero-point of the 
curve for the Cloud Cepheids so as to join the RR Lyrae variables in the Galaxy 
it zero photographic absolute magnitude. It appears that Baade’s correction 
\M 1°5 magnitudes refers to the 1930 zero-point. In 1940 Shapley (13) 
lightly revised the shape of the curve and increased the luminosity of the 
classical Cepheids by o-1 magnitudes, so that Baade’s correction to the 1940 
on of the curve should be taken as AM 1-4 magnitudes. Support for 


ABLE | 


adopted pe riod—lumnostty relation 


o’b 
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such a large correction to the zero-point of the period—luminosity curve is to 
be found in some recent investigations (14, 15, 16, 17), and in the much earlier 
work of Lundmark (18) and Mineur (5, 19). 

‘The period-luminosity relation adopted in the present work is given in 
‘Table I in which P is the cepheid period in days, and M,,, is the photographic 
median absolute magnitude defined as the average of the magnitudes at maximum 
and minimum light. 

1.2. Interstellar absorption.—'Vhe most satisfactory method of correcting the 
distance moduli of the Cepheids for interstellar absorption is by means of colour 
excesses. Given a homogeneous set of magnitudes and colours, and some 
information as to the unreddened or intrinsic colours of the Cepheids as a function 
of their periods, reliable distances could be determined. ‘lhe basic datum 
required for such a determination is the relation between intrinsic colour and 
period 

[n a recent paper (20) the writer has obtained intrinsic colours for the Cepheids 
using Eggen’s photo-electric colours, and a suitable interpolation formula for the 
reddening of the B stars in surrounding fields as a function of distance and the 
gala tic coordinates. ‘The formula 


nb 
E,(r, 5) co ( I exp ( th ) 
In D \ 


for the colour exce on the system of Stebbins, Huffer and Whitford, with a 
zero-point correction of 0°04, was used to estimate the reddening in the path 
eyond the field 


1) 
ddening formula 


to a given Cepheid by interpolation or extrapolation within o1 


stars used in the determination of the parameters in the re 
In this way, estimates of the unreddened colours were obtained, and from thes« 
estimates the best linear relations for intrinsic colours at maximum and minimum 
light as a function of the logarithm of the periods were determined hese 
relations may be used to determine the interstellar absorption for the calibrating 
stars with greater accuracy for a single star than that obtained directly from the 
formula, and to obtain reliable corrections for absorption for more distant 
( epheids. We have howeve! used Pp! ictically all the availabk photo electri 
data in the calibration of the intrinsic colour In the circumstances, pending 
the completion of further photo-clectric work on cepheid colours, the only cours« 
to follow in order to keep the distance scale as uniform as possible is to determine 
the interstellar absorption for the remaining Cepheids by the method used in 
the determination of the intrinsic colour Accordingly, the cepheid distance 
moduli have been corrected for interstellar absorption using (i) A, = yE, where 


, 75 and k is the colout excess on the P V vstem, and (in) the olution by 
+") ’ 


sive approximation in the distance modulus equation witl 
Gy rsinb 
I exp 

ind ; 
using B= 100 pe, and the values of a, given by Parenago (21) for small zone 
bounded by circles of galactic latitude and longitude he distances of the 
Cepheids determined by these methods are given in ‘lable I] 

\s a check on the absorption, it may be shown (22) that the distances of the 

Cepheids in the table satisfy the test proposed by Bottlinger and Schneller (23) 


that the half-thickness of the z-distribution should not correlate with heliocentric 
distance. It has already been demonstrated by the writer (20, p. 330) that the 
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cepheid distances calculated from colour excesses, using the intrinsic colours 
determined as above, also satisfy the condition that heliocentric distances and 
periods of Cepheids are likewise not unduly correlated. 

1.3. Observational data.—Vhe data required in the analysis of the radial 
velocities are summarized in ‘Table I] which gives the galactic coordinates (/, 6) of 
the Cepheids referred to Ohlsson’s pole (a = 190°,6= + 28°); the median apparent 
magnitudes m,, from Eggen’s photo-electric observations where available, 
otherwise the revised magnitudes given by Parenago (6); the median absolute 
magnitudes M ,, from ‘lable I, using the periods given in the Russian Catalogue 
of Variable Stars; the heliocentric distances r in kpc; the normal radial 
velocities p from Joy’s work (1, 2) revised by Parenago (6) by means of improved 
ephemerides, and from the writer’s radial velocity observations of cepheid 
variables in the Southern Hemisphere (10). ‘The data on cepheid distances and 
radial velocities should be supplemented by the accurate proper motions of 
eighteen variables published by Blaauw and Morgan (17). 

[he explanation of the notes in the last column of ‘lable II is as follows: 


absorption determined from colour excess; 2, absorption determined by 


iG. 1 Distribution of cepheid variables with known radial velocities and within 3 
n projected on the galactic plane The scale of the projected heliocentric distan 
hy complete circles centred on the Sun with radii 1-0, 2-0, and 4:0 kpc, the galactoce 


cale being given by arcs of circles drawn from the galactic centre distant 8-2 kpc from the 
longitude 3425 


successive approximation in the formula for the distance modulus; 3, Cepheids 
used in the solutions for the solar motion and Oort’s constant A in Section 2 
of the paper; 4, Cepheids with accurate proper motions used in the solution 
for the velocity of rotation of the system in Section 3.2; and 5, Population type II 
Cepheids. Six of the variables in the table are assigned to this class: W Virginis, 


and five Cepheids provisionally identified as type II by Eggen. ‘The distances of 


these objects have been calculated for an assumed period-luminosity relation 
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1°5 magnitudes fainter than that given by ‘Table I for all periods. ‘l'ype | 
galactic Cepheids within 3 kpc of the Sun are shown projected on the galactic 
plane in Fig. 1. 


‘TABLE II 


Cepheid data for analysis 


Star mpg Mpg 


V336 Aql . ; 10°8 
SZ Aql . 3°5 9° 
TT Aql 
FN Aql 
BH Oph 

7) Aql 
BB Her 
FM Aql 
BL Her 
AP Her 
FF Aql 
KL Aql 

o oge 

U Vul 
C'N Lyr 

X Vul 
SV Vul 
SU Cyg 
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/ My, Notes 
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BY 
VV 
VX 
SI 
UX 
SZ 
VY 
UY 
RW 
RI am 
RX Cam 
AS Per 
SX Per 
SV Pes 
SY Au 
AN Aur 
RX Aur 
Y Aur 
AW Per 
YZ Aur 
AQ) Aut 
SZ ‘lau 
AS Au 
KT Aur 
AA Gem 
RZ Gem 
SW ‘Tau 
S'T ‘Tau 
AD Gem 
C Gem 
CR Or; 
RS Or 
W Gem 
Cs Ori 
WW Mon 
T’ Mon 
SV Mon 
TZ Mon 
SZ Mon 
‘TX Mon 
XX Mon 
AC Mon 
CMa 
TV CMa 
TW CMa 
CMa 
CMa 
Pup 
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tar : M, Notes 
X Pup 
WW Pup 
SS CMa 
WX Pup 
WY Pup 
AD Pup 
WZ Pup 
VZ Pup 
AQ Pup 
RS Pup 
AT’ Pup 
AP Pup 
AH Vel 
RZ Vel 
AX Vel 
T’ Ant 
T Vel 
SX Vel 
B Dor 
BG Vel 
V Car 
V Vel 
/ Car 
UX Car 
Y Car 
SV Vel 
VY Car 
U Car 
ER Ci 
GI Car 
yee 
V4i9gC 
AZ ¢ 
UZ Cen 
R'T’ Mus 
T Cru 
S Mus 
R Cru 
AG Cru 
R Mus 
X Cru 
Ss Cru 
V378 Cen 
CPD — 60 
XX Cen 
V381 Cen 
V 339 Cen 
V Cen 
R'TrA 
W Vir 
S TrA 
U TrA 
» Nor 
« Pav 


AL Vir 


Ww 


=) 


t 
rauwn 
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b ) ) v Notes 


In addition two parameters, the distance to the galactic centre Ry and the 
longitude of the centre 4, are required in the analysis. We shall adopt the 
distance R, = 8-2 kpe obtained by Baade (24) from the maximum in the line of 
sight distribution in the density of RR Lyrae variables in the vicinity of the 
galactic nucleus Although the value of R, determined from the distance to 
the maximum density depends upon an adopted longitude of the galactic centre, 
probably the most important factor affecting the certainty of the result is the 
absorption. No attempt will be made in this paper to determine the longitude 
of the galactic centre from the radial velocity data. Previous determinations 
of 1, from cepheid velocities have given values between 325° and 326° with an 
r.m.s. error of about 2°. More recently, higher values of 4, between 327° and 
328° have become popular, due largely to the influence of the radio determinations 
(25). However, current work (26, 27) suggests that the radio position may be 
disturbed by a point source approximately in the direction of the centre but 


located between the Sun and the galactic nucleus. Since the values of /, suggested 


by the optical and radio data lie within the range /, = 325° to 328”, we shall carry 
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out the analysis with these two assumed longitudes in order to set limits to the 


values of the parameters in the representation of the cepheid motions which are 


sensitive to the direction of the galactic centre 


2. Differential rotation near the Sun 

In a previous paper (20) a determination of Oort’s constant was given for 
a group of nearby Cepheids with distance moduli corrected for absorption by 
means of colour excesses. Using the Oort approximation to the effect of galactic 
rotation on the radial velocities, it was shown that in spite of Baade’s revision 
of the zero-point of the period-luminosity curve, Joy’s original value of Oort’s 
constant A = 20 km/sec/kpe was recovered. ‘The interstellar absorption adopted 
by Joy appears to have been unduly low by comparison with the absorption 
required by the reddening of the B stars, and this higher absorption was shown 
to compensate for the increase in the distance moduli of the Cepheids brought 
about by Baade’s revision of the absolute magnitudes 

‘The solutions for differential rotation given in this section bring the earlier 
analysis forward in two respects. In the first place, it was previously necessary 
to correct the radial velocities for local solar motion since the amount of data 
was insufficient for a realistic determination of the velocity components 
Secondly, whereas the use of the Oort approximation was a necessary part of 
the scale-factor method, it is now convenient to use higher order approximations 
which do not restrict the selection of stars to those for which the heliocentru 
distances are small. 

2.1. The radial velocity equation._-Consider a galaxy containing a highly 
flattened system of stellar objects in differential rotation about an axis of 
symmetry, and let w(R) be the angular velocity at a distance Rk trom the axis 
of rotation, w(R) being assumed to be independent of the distane« above or 
below the galactic plane for small values of If circular orbits of the objects 
about the galactic centre are assumed, the motions take place in a field of fore 
F(R) = Rw*(R) directed towards the galactic centre Ihe radial velocity of an 
object relative to a frame of reference rotating with the velocity of the system 


at the distance R, is then given by 
(R,)| sin (1— 1,) cos b (1) 


here / and 4 are the galactic co-ordinates of the object ‘This form for the radial 
elocity due to ditterential rotation was originally introduced by Bottlinger (28) 
in order to obtain, for empirical laws of force, trigonometric expansions of a 
higher order than is given by the Oort approximation 


insion of th galacti ystem on the 


low superimpose a general exp 

tation neglecting the compone nt of the elocity field of the « pansion 
init the effect of thi component on the radial velocits involve the factor 
int, which ts negligibly small for objects in a highly flattened tem Let Il 
be the 1 mmponent of the velocity of expansion at a distance FR from the axi 
of rotation, and suppose that for small distance from the galacti plane 1) i 
effectively independent of he radial velocity ic to the expansion ts then 
iven by 


ki R\r cos* bh RRR) k(R,)| co (/—1,)cosbh (2) 


where A(R)=11/R Ihe expression (2) | the A-term for a general ¢ x pansion 
parallel to the galactic plane 
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Combining the effects due to differential rotation and expansion of the 

ystem, the radial velocity equation may be written in the form 
p' = krcos* bh + f(R, R,) sin (1— 1,) cos b — g( R, Ry) cos (1— ly) cos 6, 
whe re 
[(R, Ry) = Rof[w( R) — oA R,)), o(R, Ry) = Ro[k(R) — RCRo)). 

‘he radial velocity p’ is the velocity of an object relative to a frame of referenc 
rotating with the system. Let Uy, UV, and w, be the components of the solar 
motion referred to this frame, the components wu, and Wy being in the galactic 


plane and in the direction of longitude Le and l, + gO respectively, Wy being 


perpendicular to the galactic plane and directed towards the north galactic pole. 


I‘he observed radial velocity corrected for the solar motion is then given by 


p' =p +u,cos(l—1,)cosb +7, sin(1—1,) cos b + w, sin b. (4) 
he distinction between the radial velocity p referred to the Sun, the radial 
velocity p’ defined by equation (3), and the radial velocity p, given by equation (1) 
for purely circular motion should be noted. ‘lhe latter velocity is clearly th 
observed velocity corrected for solar motion and the A-term. 
2.2. Solutions for galactic rotation.—We wish to make an analysis of the 
cepheid velocities which will lead to a determination of the Oort constant 


1 ; (a R =») . (<) 
dR)», 


the velocity components of the solar motion, and to an examination of the level 
of significance of the A-term for the cepheid system by means of equations (3) 
and (4), using a suitably approximate treatment of the equations. ‘The passage 
from the exact form (1) for the radial velocity equation for circular motion to the 
well-known Oort approximation has been shown by Camm (4 p. 166; 11, 
p. $2) to involve two steps, (i) a dynamical approximation 
f(R, R,) 2A(R—R,), (6) 
and (11) a geometrical approximation 
R R,, rco (/ l,) cos b. (7) 
Although (1) ts satisfactory for small values of R—R, so long as the curvature 
of {(R, Ro) is not large, (11) is rather less satisfactory since it involves a systematn 
error. Furthermore, if there is a change 4r in r, the corresponding change oR 
in KR is such that [4(R—R,)| ~~ [or], and the equality only applies to objects 
and /,4 180. In most cases we have 


in the galactic plane with /=/, 
| a(R Ry) | < [dr |, and the radial velocity is then by no means linearly 
proportional to the distance r. In this section we shall adopt Camm’s linear 


approximation (6) By analogy with the definition (s) of Oort’s constant, we 


do(R, R,) 
- a» 


30 that the approximation is g(R, R,) 2C(R—R,). ‘The appropriate form 
of equation (3) 1s then 
p kor co *h—2A(R R,) sin (/ 1,)cosb 
2C'(R— R,)i(r/R,) cos? b — cos (1 l,) cos b}, (g) 


where ky = k( Ry) l hus, to a consistent order of approximation we may determine 


detine 
(8) 


from « quations (4) and (9), by the ordinary method of least squares, the quantities 
ko, A, C, up, Vg, and wy, although a determination of w, would be low in weight 


for a highly flattened system of objects. A further simplification might be mad 
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by writing C=o, corresponding to an expansion parallel to the galactic plane 
such that ///R is independent of R. We shall, however, retain the term in C 
in equation (g) in order to preserve the degree of approximation, and to decide 
from the observations if there is any evidence for the function g(R, R,) in the 
cepheid motions. 

\lthough the linear approximation to f(R, R,) in equation (6) may be used 
without restriction on the heliocentric distance r, it is necessary that R— R, 
should not be too large if the higher order derivatives of the function are not 
negligible. Some idea as to the limitation of the range of values of R— R, which 


might be necessary may be gained by considering the rotation 


f(R, Ro) = $ARg{(Ro/R)? — 1} (10) 


which holds for a system in rotation under the influence of an invers« square 


law of force. A comparison between the exact function and the approximate 
form (6) shows that the error is less than ten per cent provided |R— Ry|~ 0-1 Ry. 

In ‘Table I] there are g6 Cepheids for which the galactocentric distances, 
given by 


R,?* tr*cos*h 2Ryr cos (1 1,) cos b, 


lie within the range 7:-5kpe-R~<g-okpe, comprising about half the data 
available, the velocities of 45 C'¢ pheids being taken from the recent observations 
(10). We shall restrict our solutions of the approximate equations to Cepheids 
which lie within this annular strip of width 1-5 kpc, the Sun being located 
approximately at its centre. ‘he solutions were made by the ordinary method 


of least squares from 96 equations of condition of the form 


kr cos* b — u, cos (1— 1) cos b — 2, sin (1 — 1) cos b 
w, sinb—2A(R— R,)sin(1—1/,) cosh 
2C(R— R,)}(r/R,) cos* b— cos (1— 1,) cos b} = p+ dp, (11) 


where dp, the peculiar part of the radial velocity, is a correction to the observed 
radial velocity p which includes the observational error in p and departures 
from the adopted kinematics as represented on the left-hand side of equation (11). 
We have assumed that dp is more important than the errors in r and R— Ry, 


a condition necessary for the applic ation of the ordinary method of least quares, 


Solutions for Oort constant 1, the ¢ mponent 
of the solar motion, and the K-term 
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and that op behaves like a random error in the construction of the normal 


equations. ‘The velocity component w, is essentially indeterminate, and the 


0 
tandard value uw, 7°3 km/sec has been adopted. 

».%. Discussion of the solutions.—'Vhe solutions for Oort’s constant, the solar 
motion, and the A-term are given in ‘Table III for the two assumed longitudes 
of the galactic centre, l,= 325 and 328 ‘The table also gives the r.m.s. erro! 
calculated in the usual way from the weights of the unknowns in the normal 
equations and the r.m.s. error » of an observational equation. ‘Three solutions 
were made for each longitude. In the first, all six unknowns were determined, 
in the second we have taken C=0; and finally, a solution was made without 
any A-term 

». 2.1. The K-term he solutions show that there is no evidence in th 
cepheid motions for the second part of the K-term which involves the coefficient 
( Hence o(R, Ry) 0, and k(R)=k,, a constant independent of R, is the form 
of the A-term to be examined. It is evident from the size of the r.m.s. error 
in f 


warrant for retaining the term if 4,=325 . However, the size of the negative 


» that we are dealing with a quantity of marginal significance. ‘There is no 


K-term increases with /,, and thus could become more significant if higher 


values of /, wer adopted For / 228 , we have k 


. . 2°6 + 1:1 km) sec/kpe 
\lthough it might be argued that a quantity which is two and a half times the 
r.m.s. error could be significant, the fact remains that the velocity dispersion 
is only very slightly reduced if the A-term is retained. It might be concluded 
that there is some slight evidence in the cepheid motions for a negative A-term 
corre ponding to a uniform contraction of the system if the longitude of the 
galactic centre is 328 ‘The acceptance of a negative AK-term for the cepheid 
ystem has been advocated by Parenaygo (6, 29) who has considered some ot 
the interesting consequences which arise from such a term in dynamical theory 
However, it appears that the existence of the term was inferred from solutions 
for galactic rotation in which the lower value of the longitude /,= 325° was 
assumed Lhe present work offers little support for the A-term in these 
circumstance 
Solar motion.—‘Vhe solutions for the solar motion in ‘Table ILI are the 
velocity components relative to a frame of reference moving with the circular 
velocity of the ystem at the olar distance trom the galactic centre, [Che solar 
motion relative to the circular velocity should be carefull distinguished from 
the local solar motion which is the velocity of the Sun with respect to the centroid 
of the space motions of the nearby stars. Now the well-known asymmetry in 
the frequency distribution of the space motions in the direction of galactic 
rotation gives a centroid which ditlers for different groups of stars Sut it has 
been pomted out by Vyssotsky (30) that whereas the centroids of the space 
motions for the various groups differ, the maximum of the frequency function 
for each group occurs at approximately the same value of the space velocity 
The position of the maximum has been identified by Vyssotsky with the solar 
motion with respect to the circular velocity. In order to distinguish between 
the solar motion detined in this way, and the /ocal solar motion, we shall use 
Vyssotsky's term basic solar motion. 
\ comparison between the velocity components for the standard value of 


the local solar motion (20 km/sec towards /= 23-5, b= +21 -6), the basic solar 


motion of Vyssotsky and Janssen (30), and the velocity components obtained 





No. 4, 1956 of the system of cephetd variable stars 405 


in the present work is given in ‘Table IV. In addition to the velocity components 


Uy, Up, and wp, the table gives (uy? + v,")'*, the projection of the solar velocity on 
the galactic plane Special consideration should be given to vp, that 1s the 
component in the direction of galactic rotation. It is found that the solar motion 
obtained for the Cepheids gives a rotation of the system which makes the averag« 
value of the function f(R, R,) equal to zero for Cepheids for which R= Ry. Camm 
has shown (11, p. 84) that this condition is rather sensitive to the value of the 
vy-component. If the condition is not satisfied, the @,-component of the solar 
motion is incorrect. A recent determination by Raimond (31) is deficient in 
this respect as it gives f[(R,, R,)= +7 km/sec. ‘The discrepancy in f(R,, Ry) 1s 
not so marked in the case of the local solar motion 

lhe last line in ‘Table IV shows that the components of the solar motion 
obtained in the present work are practically the same as the components of the 
basic solar motion given by Vyssotsky and Janssen. Since the motion of the 
Sun with respect to the circular velocity of rotating systems of objects 1s a 
reflection of their speed of rotation, the system of stars to which Vyssotsky and 
Jansse n's determination refers must be sharing the ime rotation as the cephe id 


svstem at the solar distance from the galactic cent: 


TABLE IV 


(‘omparison hetween velocity compo 
of the solar motion for ly 325 


(Component Ng 


Loc il 
Motior 


yer 
Solar 
2.2, Oort’s constant In an earlier paper (20) Oort’s constant was determined 
from a group of nearby Cepheids at a mean distance of 0°57 kpe from the Sun 
using distance moduli corrected for absorption by means of colour excesse 


The solution A =19°8 2 km/sec k pe was obtained by means of the Oort 


approximation Since the reddening of the B star hich was used tor the 


determination of the intrinsic colours in (20) is essentially the same as that used 
in the calculation of the distances of the Cepheids in lable [1 tor which colours 
are not available, the determinations of Oort’s constant given in lable Il] and 
in the previous paper differ only in method and weight lhe solution for 
Qort’s constant in ‘lable III] for the case /, 22 and no A-term, namely 
1g°5 + 1°9 km/sec/kpe, is practically identical with the value of 1 found in the 
earlier paper by means of the Oort approximation with / 225 and no A-term 
‘The present work theretore strongly reinforces the earlier conclusion that Baad 

revised distance moduli for the Cephe id corrected for ab orption by means of 


the reddening of B stars in surrounding field rive in Oort constant 
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1=20 km/sec/kpe approximately. ‘These determinations refer to a distance 


cale in which the average absorption 2A,,/2r is approximately 2 magnitudes 


7 


per kpe near the galactic plane. 


3. The rotation of the cepheid system 

In the preceding section, Camm’s linear approximation (6) was used for the 
functions f{(R, R,) and g(R, R,) for small values of R— Ry. ‘This step was taken 
o that the differential rotation at the solar distance from the galactic centre 
could be determined at the same time as the components of the solar motion 
vith respect to the circular velocity of the rotating system. In this way we 
have obtained a solar motion which gives the correct rotation /(R, R,) at R= R 
and a determination of the A-term which could be used to free the observed 


th 


radial velocities from the effect of non-circularity in the mean motion about the 
valactic centre. 

Having obtained the corrected velocities, the restriction on R— R, may now 
be removed, and the rotation of the entire system examined by means of a powerful 
general method introduced in 1938 by Camm (11). ‘The method starts with the 
radial velocity and proper motion equations for circular motion as written by 
Bottlinger (28) for any force function. But whereas Bottlinger only used the 
equations to obtain higher order terms in trigonometric expansions for special 
laws of force, that is for assumed forms for the rotation as a function of distance 
from the galactic centre, Camm inve rted the proble m and determined the rotation 
as a function of distance directly from the radial velocities. ‘The method was 
applied by Camm to cepheid variables (4) and planetary nebulae (11), and has 
been used extensively by Soviet astronomers with due acknowledgment. It 
appears, however, that the existence of the method has not been appreciated 
elsewhere 

3.1. Exposition of Camm’s method.—Vhe radial velocity and proper motion 
equations may be written in the form 


p, = {UR, R,) sin (1— 1,) cos b, (12) 


and 


{(R, R,)| cos (1— 1,) — (r/ Ry) cos b| — w( Ry )r cos b, (13) 


where p, and 7, are the radial and transverse velocities corrected for solar motion 
and the A-term, the proper motion in galactic longitude y,cosb being related 
to 7, by the equation 


rp, COs b u, sin (/—1,) — va cos (1 — 1,) + g(R, R,) sin (/—1,). (14) 


I. quations (12) and (13) are exact in the sense that they follow from the assumption 
of circular motion about the galactic centre without mathematical approximation. 
When applied to a system of stellar objects, a velocity dispersion is superimposed 
on the mean motion. Equations (12) and (13) should therefore contain terms 
representing the peculiar part of the radial and transverse velocities. 

Camm’s method may be divided into two steps, (1) the determination of the 
{ (R, R,) curve, (2) the determination of w(R,) from the proper motion equation, 
and hence the angular rotation of the system as a function of & from the equation 


w( R) = w(R,) + Ry f(R, Ro). (15) 


Ihe determination of the {(R, Rk,) curve for the system consists in finding mean 


values of f(R, Ry) and R tor objects conveniently arranged in distance groups. 
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The following is a suitable method of weighting equation (12) to obtain the best 
mean value of {(R, R,) from the objects within each group. If the range of vatues 
of R in the group is small, we may assume that 


f(R, Ro) =f(R,, Ro) +(R—R,Vdf/dR) pn, (16) 
where R, is some weighted mean value of R. Let w=sin(/—/,)cos6, then 
equations (12) and (16) give 

f,w + (df,/dR)w( R — R,) =p, + op,, (17) 
where /,=/(R,, Ry), and dp, arises from the velocity dispersion. If dp, can be 
treated as a random error, the ordinary method of least squares applied to 


equation (17) gives the normal equations 
f, Xw? + (df, /dR)Xw*(R — R,) = Lwp,, (15) 
f 2% 2(R R,) { (df, dR) <w*(R R, 2 — Vw(R R,)p ; (19) 
Now the weight of the solution for f, in these equations is 
Lw* — (Xw?( R — R,)]?/Xw*(R — R,)’*, 
and when Xw*(R— R,)=0, that is when R, is given by 
I Lw*R/Xw?, (20) 
the weight of {(R,, R,) takes its greatest value Lw?. ‘The mean value of /(R, Ry) 


given by equation (18) is then 


{(R,, Ry) = Lwp,/Uw* (21) 
‘This method of weighting takes two factors into account, the variation of /(R, Ry) 
within the group, and the weight of each object as a function of the longitude. 
If the velocity dispersion were small and the amount of data adequate, the method 
would give small segments of the f(R, R,) curve It is, however, sufhicient to 
have established that the appropriate mean values for any distance group are 
given by equations (20) and (21),* the /(R, R,) curve for the system being the 
continuous curve through the mean points The velocity-dispersion tor the 
objects in each group follows directly from the radial velocity equation (12), the 
value of f{(R, R,) for each object being taken from the smooth curve lhe r.m.s. 
error p, in the weighted mean value {(R,, Ro) mav then be calculated for each 
group from the relation 


jy" dU(op )° (1 1 Liz ; (22) 


Having determined the {(R, R,) curve tor the ystem, the econd step, the 


determination of w(R,) from the proper motion equation, may now be completed 


from equations (13) ind (1 }) with o(R, R.) O, have 


wk.) cosh J rit, COS b, (22) 
\M he re 
u, sin (/ in) v, cos (/ /,) {(R, R,)lco (/ /.,) () R,,) cos b}. 


Since U,, and f{(R,R,) are known from the analysis of the radial velocitic 


the contribution V to the mean motion transverse to the line of sight may b« 
found for any object. ‘he solution of equation (23) for w(R,) which then follow 
requires some consideration since the terms in the equation are subject to several 
sources Of error In the first place, for a system of objects with a velocity 


dispersion, equation (23) should contain a term representing the peculiar 


le of t tre and the 
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transverse motion which affects the term V. Secondly, the observational term 
in the equation of condition contains the observed proper motion j,cosb in 
salactic longitude and the computed distance r, both of which are subject to 
errol If the distances are known more accurately than the proper motions, 
equation (23) may be solved by the ordinary method of least squares to obtain 
the estimate 


w(R,) = Xr cos b( V — rp, cos b)/ Zr? cos* b. (24.1) 


‘his is equivalent to weighting the individual values of w(R,) in equation (23) 
by r*cos*h, Alternatively, equation (23) may be divided by r before applying 
least quares to obtain a second estimate 


w»( Ry) = X cos b(V /r — pz, cos b)/X cos? b. (24.2 


In this case, greater weight is given to the nearer objects by placing emphasis 
upon errors in V due to the velocity dispersion. Since equations (24.1) and 
(24.2) give estimates of w(R,) which fit the data from one end of the distance 
range at the expense of the other, it is advisable to make a third estimate to 
compromise between the two extremes. ‘This may be done by introducing 
weights proportional to rcoshb, so that the weighted mean value of w(R,) from 


equation (23) 1s given by 
wo Ry) = UV — ry, cos b)/&r cos b. (24.3) 
As a further refinement, the observational equations might be weighted according 
to the errors in the proper motion determinations, choosing weights proportional 
to the inverse of the errors in preference to their inverse squares. 
3.2. Application to cepheid variables. —\n Camm’s application of the method 


to cepheid variables and planetary nebulae (4, 11), preliminary forms for the 


{(R, Ry) curve were obtained for assumed values of the velocity components 


(Uy, Vp, Wy) Of the solar motion, the longitude /, and the distance R,. ‘The method 


of differential corrections was then applied to the radial velocity equation in order 
to determine corrections Au,, Avy, Awy, Al, and AR, to the assumed values by 
the method of least squares, the preliminary curve for f/(R, R,) being used to 
obtain the derivative Of(R, R,)/AR required in the equations of condition. ‘This 
refinement of the method has not been applied in the present work in which the 
{(R, R,) curve for the Cepheids is based upon an adopted longitude /,, Baade’s 
given by th 


value of the distance ms and the velocity components u and 


0 “Oo 


olutions carried out in Section 2.2. 

Cepheids listed in ‘Table I] as members of Population Il have been omitted 
in this determination of the /(R, R,) curve. ‘The remaining 18g Cepheids wer 
conveniently grouped according to distance R, and weighted mean values of R 
ind /( 2, &,) tor each distance group were calculated from « quations (20) ind (21), 


the I idial vi lox ity p for each object be ny viven by 


p p u, cos (1 i,) cos b in(l l,)cosb W, 81nd K, 


where p is the observed radial velocity in ‘Vable II lhe mean values of R and 
{(R,R,) are given in ‘Table V for / 325 and no A-term, the values of wu, 


0 0» 


and w, being given in the third column of ‘Table III. 

Having determined the f(R, Ry) values we can now proceed to the determination 
of «w(R,) for the cepheid system along the lines indicated in the preceding section 
It is to be expected that the determination cannot be made with much w ight 


on account of the small number of reliable proper motions which are availabl 
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at the present time. It is, however, essential to attempt a determination of the 


angular velocity, even from insufficient information, for two reasons. In the 
first place, this can only be done satisfactorily for distant objects—for there may 
be a local vortex—and the differential rotation effect should be carefully removed, 


TABLE \ 


Weighted mean values of R and f(R, Ro) 
for cephe id variables in distance groups 


in R | Number of | R 
objects 


as indicated in Section 3.1, and not merely thrown into another harmonic term. 
Secondly, there is no valid reason for supposing a priori that the angular velocity 
«w( Ry) for the Cepheids should be the same as that given by other groups of object 
Indeed, it has been established observationally by Stromberg (32) and envisaged 
theoretically by Lindblad (33) in the original hypothesis of galactic rotation, 
that the velocity of rotation is closely related to the velocity dispersion of each 
group or system of objects. Now the Cepheids comprise a highly flattened 
system in which the random motions of the objects, as indicated by the dispersion 
in radial velocities, are much smaller than the random motions of other group 
of objects such as the early type stars. Although the proper motion data are 
much more extensive for these groups than for the Cepheids themselves, it 1 
clearly desirable to determine w(R,) for the cepheid system from the cepheid 
data alone, instead for forcing agreement with othe groups of objects. 

‘The data used in this determination of w(R,) for the Cepheids are the prope: 
motions of 18 variables given by Blaauw and Morgan (17) in their work on the 
calibration of the cepheid absolute magnitudes. ‘The proper motion stars are 
indicated in the notes for ‘Table II which also indicate that the distance moduli 
of ten of the objects have been corrected for ab orption by means of colour 
excesses. ‘lhe distances of the remaining eight Cepheids for which photo-electric 
colours are not available are probably less accurately known. ‘The component 
of the proper motion in right ascension and declination given by Blaauw and 
\lorgan were combined to give the proper motion in galactic longitude in second 
of arc per annum. ‘The linear transverse velocity «rjy,cosh was then obtained 


for cach object in km/sec., using « = 4°735 = 10 with the distance r in kpe. ‘The 
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‘TABLE VI 


Estimates of ow Ro) for the cepheid syste m 


Estimate | Solution 1 Solution 
km 
39°24 4 33%. 


/ 


| km /sec/kp« 
| 
w (Ro) ie 
i4( Eo) 


( 4! Ro) 


26 gt 6's 2QY 


28°5 + Bo 32 





| 
| 


Vlean 85°47 7°2 32 
three estimates of w(R,) from equations (24.1), (24.2) and (24.3) are given in 
‘lable VI with the r.m.s. errors. ‘lhe second set of solutions given in the table 
vas obtained by weighting the equations inversely as the errors in the proper 
motions. ‘The solutions were all made for the case /,= 325°, using the f(R, R,) 
curve defined by the mean points in ‘Table V. 

The weighted mean value of the estimates given in ‘lable VI 1s 
u( R,) = 32 +7°5 km/sec/kpc, the large uncertainty being due entirely to the 
mall number of reliable proper motions. It is of interest to note that the same 
proper motion data which enabled Blaauw and Morgan to determine the zero- 
point of the P-L curve with considerable accuracy have given a determination 
of the rotation with low weight. ‘The uncertainty in w(R,) is such that the Oort 


constant B 12°5 + 7°5 km/sec/kpc, which follows from the relation 
B= A—w(R,) 


with A= 19°5 km/sec/kpe, does not differ significantly from the lower negative 
value B 6-9 + 1°5 km/sec/kpe obtained by Morgan and Oort (34) from the 
proper motions of early type stars and faint stars in general, and the entirely 


TABLE VII 


(Comparison of {(R, Ro) for cepheid variables with values from 21-cm radio data and schemati« 
models of the galactic system 


Cepheid Radio oI Schmid 
ariables data le] mod 


11d 11d 
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independent value of B which follows from the axial ratio of the velocity ellipsoid, 

B/(A — B)=0°24 + 0°04, obtained by Hins and Blaauw (35) from a study of 
the Radcliffe proper motions of faint stars. In the course of time, the accuracy 
of the determination tor the Cepheids will improve as more data on proper motions 
become available. It may then be possible to ascertain if the Oort constant B 
for the Cepheids differs appreciably trom these average values for heterogeneous 
collections of nearby objects. 

3.3. Discussion of the rotation curve.—We shall conclude this discussion of 


the rotation of the cepheid system with a comparison between the observed 


rotation and the rotation according to schematic models of the mass distribution 
in the galactic system. If the differential rotation of the system near the Sun 
is determined largely by the mass of the Galaxy in the inner region where R~— Ry, 


the f(R, Ry) curve may be approximated near R= R, and at greater distances by 
equation (10), namely 


{(R, Ry) =4AR,}(R,/R)? 


‘This gives the rotation for circular motion when the field of force F(R) varies 
as the inverse square. Although this representation appears to involve only 
the distance R, and the Oort constant A, it is implicit in the assumption of an 
inverse square law that A + 3B=o, so that the value of the Oort constant B 1 
implied for an assigned value of A. Now the value A = 19°5 km/sec/kpe for the 
cepheid system Is well determined from the analysis of the radial velocities, 
and a system with this value of A in an inverse square field of force would have 
B 6-5 km/sec kpe. ‘The cepheid proper motions, on the other hand, give 
an uncertain but greater negative value B 12°5 km/sec/kpe. Although the 
inverse square behaviour of the field of force near the Sun would be excluded 
by this greater negative value, the uncertainty is such that a decision cannot be 
made from the available proper motion data. Values of /(R, Ry) calculated 
from equation (10) with A= 19-5 km/sec/kpe are given in ‘Table VII for the 
values of R at which weighted mean values of {(R, Ry) were calculated for cepheid 
variables in Section 3.2. A comparison is also made in the table with values of 
{(R, Ry) for circular motion in Oort’s original model of the galactic system (36, 
37, p- 134) in which the distribution of mass is represented by the superposition 
of concentric homogeneous spheroids and a more recent model, constructed by 
Schmidt (38), in which both homogeneous and non-homogeneous spheroids 
are used. An important feature of the model adopted by Schmidt is that in the 
region R-—K, it reproduces the circular velocities obtained from the Leiden 
21-cm radio data (39). ‘These models are constructed for a Galaxy in which 
Ry =82 kpc, A=19°5 km/sec/kpe, B 6-9 km/sec/kpe, so that the circular 
velocity at the Sun is ©, =216 km/sec. Accordingly, the field of force at R= R 
is practically an inverse square as A +3B 1:2 km/sec/kpe for the models, 


whereas A+3B=0 for the inverse square law. Since most of the mass 1 
concentrated in the region R < R,, it is perhaps not surprising to find in ‘Table VII 
that for RR, the rotation f(R, R,) for the models does not differ greatly from 
that given by the inverse square law. 

A comparison Is given In Fig, 2 between the rotation for Schmidt's model 
and the mean rotation calculated for the cepheid groups. ‘The comparison show 
that the best agreement is to be found in a comparatively small range of distances 
about R= R, where the agreement is guaranteed by the appropriate choice of 
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The differential rotation {CR, R,) of the ystem of cephoid aniable star mpared 
lation curve for Schmidt's model Open circles denote mean rotations comput d for cepheid 
iniahles in distance groups by Camm’s method, assuming |, 32 and RR, 8:2 kpe The O 


constant A for the model is 19 km sec kp 


the Oort constant A for the model. ‘There are, however, three remarkable 
discrepancies between the rotation /(R, R,) for the Cepheids and the model. 
\t large distances from the galactic centre the Cepheids are found to rotat 
much more rapidly than the model. ‘There are, no doubt, appreciable uncertainties 
in the absorption at great distances in the absence of data on the colours of th 
faint Cepheids which comprise the outermost groups, but it appears unlikely 
that the cepheid rotations would be drastically modified by a revision of the 
distances of the objects in these groups. ‘The second discrepancy between the 
Cepheids and the model is the appearance of large negative values of f(R, Ro) 
near R=g:2kpe. ‘This peculiarity is also present in the rotation curve obtained 


by Camm (4) and is due, as Camm suggested, to group motions In a small sector 


of galactic longitude in the vicinity of the constellation Cassiopeia. Ihe velocity 
dispersion within the group itself is small, but the radial velocities are systemati 
cally too negative by about 20 km/sec. ‘Uhe third disc repancy 1s perhaps the 
most striking. For distances R <7 kpc, which are not shown in Fig. 2, the cepheid 
rotations given in ‘l'able VII fall well below the rotation curve from 21-cm radio 
data. 
In parti ular, the group near¢ st the galactic centre shows a striking difference 
It is rather remarkable that there should be such a large discrepancy between 
{(R, R,) from cepheid variables on the one hand and radio data on the other 
as both methods purport to give an observational determination of the same 
quantity, In the case of Camm’s method it is possible to assign an uncertainty 
to the computed rotation by means of equation (22). ‘The computed uncertainty 
for the two groups nearest the galactic centre is large, partly on account of the 
mall we ight which the method a igns to the objects near the longitude of the 
ilactic centre, and partly because the velocity dispersion appears to increas« 
towards the centre of the Galaxy. In the case of the radio method no uncertainty 
is assigned to the computed rotation, but the determination by Camm’s method 
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clearly indicates that the rotation for the group nearest the centre lies significantly 


below the rotation given by the radio determination 
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